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In  conducting  the  research  described  in  this  report,  the 
investigators  adhered  to  the  "Guide  for  Laboratory  Animal  Facilities 
and  Care"  as  promulgated  by  the  Committee  on  the  Guide  for  Laboratory 
Animal  Resources,  National  Academy  of  Sciences  -  National  Research 
Council . 
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ELECTROMAGNETIC  ENERGY  ABSORPTION  AND  ITS  DISTRIBUTION  FOR  MAN 
AND  ANIMALS  AT  DIFFERENT  FREQUENCIES  UNDER  VARIOUS  CONDITIONS 


Objectives 

The  objectives  of  the  project  were: 

1.  To  quantify  the  electromagnetic  power  absorption  and  its 
distribution  for  man  when  subjected  to  radiation  at  dif¬ 
ferent  frequencies  and  under  various  exposure  conditions. 

2.  To  verify  the  important  findings  of  part  (1)  by  exposure  of 
living  organisms. 

The  experiments  for  man  were  performed  with  reduced  proportion¬ 
ately-scaled  models  filled  with  0.9  percent  saline  solutions  and  with 
biological-phantom  materials  developed  to  simulate  electromagnetic 
properties  of  human  tissue.  The  animal  experiments  consisted  of  micro- 
wave-induced  convulsions  and  work  stoppage  in  rats. 

The  purpose  of  the  project  was  to  develop  an  understanding  that 
would  lead  to  projections  for  humans. 

Highlights  of  the  Work  Done  on  the  Project 

* 

A.  Electromagnetic  Energy  Deposition  in  Man  and  Animals 
Free-Space  Irradiation  Condition 

The  condition  that  has  been  studied  the  most  to  date  is  that  of 


A  paper  based  on  these  highlights  will  be  submitted  for  publication  in 
the  Proceedings  of  the  IEEE  (January  1980)  Special  Issue  on  Biologi¬ 
cal  and  Ecological  Effects  and  Medical  Applications  of  Electromag¬ 
netic  Energy. 
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free-space  irradiation  of  single  animals.  The  whole-body  absorption 
of  electromagnetic  waves  by  biological  bodies  is  strongly  dependent  on 
the  orientation  of  the  electric  field  (E)  relative  to  the  longest  dimen¬ 
sion  (L)  of  the  body.  The  highest  rate1  of  energy  deposition  occurs  for 
E 1 1  L  orientation  (see  Appendix  A),  for  frequencies  such  that  the  major 
length  is  approximately  0.36  to  0.4  times  the  free-space  wavelengtn  (A) 
of  radiation.  Peaks  of  whole-body  absorption  for  the  other  two  configura¬ 
tions  (major  length  oriented  along  the  direction  of  propagation  (k||  L) 
or  along  the  vector  of  the  magnetic  field  (H ] |  L)  have  also  been  reported2 
for  A  on  the  order  of  4Tib,  where  2irb  is  the  weighted  average  circumfer¬ 
ence  of  the  animals. 

Using  prolate  spheroidal  and  ellipsoidal  equivalents  of  biologi¬ 
cal  bodies,  theoretical  calculations  have  recently  been  given  in  a 
dosimetry  handbook3  for  frequencies  up  to  and  slightly  beyond  the 

_y  *  * 

resonant  region  for  the  aforementioned  polarizations  E [ |  L,  k||  L,  and 
H ||  L.  Numerical  calculations  for  a  realistic  model4  of  man  (see  Ap¬ 
pendix  B)  have  shown  a  fine  structure  to  whole-body  absorption  at  fre¬ 
quencies  higher  than  the  whole-body  resonant  frequency.  Minor  peaks  in 
the  supraresonance  region  are  ascribed  to  maxima  of  energy  deposition 
in  the  various  body  parts  such  as  the  arm  and  the  head5  (see  Appendix  C 
for  details) . 

For  the  supraresonant  region,  the  E||  L  orientation  has  been 
studied  most  extensively.  An  average  1/f  dependence  of  the  whole-body 
absorbed  dose  is  experimentally  observed5  to  frequencies  f  on  the  order 
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of  1.6  S  times  the  resonant  frequency  f  ,  where  S  is  the  relative 

res  •  r’  res 

absorption  cross  section  at  the  resonant  frequency. 

Empirical  equations  for  specific  absorption  rate  (SAR)  have  been 
derived6  for  e||  L  orientation  and  found  to  be  fairly  accurate  for  six 
animal  species5  from  25-g  mice  to  2250-g  rabbits. 

These  equations  are  given  in  Table  1. 

The  coefficients  in  these  empirical  equations  for  man  model  are 
approximatley  63  percent  of  the  live  animal  coefficients.  Further 
research  is  needed  to  resolve  this  difference. 

Electromagnetic  Absorption  in  the  Presence 

of  Nearby  Ground  and  Reflecting  Surfaces 

Only  highly  conducting  (metallic  sheet)  ground  and  reflecting 
surfaces6’7  of  infinite  extent  have  been  studied  to  date  (see  Appendix 
D) .  For  a  standing  man  model  with  feet  in  conductive  contact  with  a 
perfect  ground,  there  is  a  drastic  alteration  in  SAR  as  a  function  of 
frequency.  For  E | [  L  orientation  the  new  resonant  frequency  is  roughly 
one  half  that  given  by  Eq.  1.  At  this  lower  resonance  frequency,  the 
SAR  is  about  twice  that  at  peak  absorption  frequency  for  free-space 
irradiation. 

The  nature  of  the  ground  effects  on  SAR  (for  e||  L  orientation) 
is  such  that  even  a  small  separation7-9  from  ground  (to  break  conductive 
contact)  is  sufficient  to  eliminate  much  of  the  ground  effect  (see  Ap¬ 
pendices  D  and  E) .  For  separations  from  ground  more  than  3-4",  the 
total  energy  deposition  and  its  distribution  are  identical  to  those  for 


Table  1.  SAR  for  man  models  for  conditions  of  free-space  irradiation 


E||  L  Orientation 


Peak  absorption  or  resonant  frequency: 


f  =  11. A/L  GHz 
r  cm 

For  subresonant  region  —  0.5  f  <  f  <  f  : 


(1) 


SAR  in  mW/g  for  1  mW/cm 
incident  plane  waves 


0.522 


L 


2 


_ cm _ 

weight  in  g 


(2) 


For  supraresonant  frequency  region  —  ff  <  f  <  1.6  Sres  f r : 


SAR  in  mW/g  for  1  mW/cm2  _  5.95  Lcm 
incident  plane  wave  fields  f  ^  weight  in  g 


where  Lcm  is  the  long  dimension  of  the  body  in  centimeters  and 


S 

res 


0.481  U : — r~c - 

\  weight  in  g 


(4) 


\ 


free-space  irradiation  conditions.  Even  for  a  man  model  in  conductive 
contact  with  a  perfect  ground,  the  energy  deposition  in  the  supraresonance 
region  (f  >  2-3  f  )  is  comparable7  to  that  for  conditions  of  free-space 
irradiation. 

Other  orientations  and/or  finite  conductivity  ground  effects  on 
SAR  have  not  been  studied  to  date. 

For  highly  conducting  (metallic  sheet)  reflecting  surfaces  of 
flat  and  90°-corner  types,  enhancements6  in  SAR  by  factors  as  large  as 
27  have  been  observed  for  e||  L  orientation.  Most  of  the  work  to  date 
has  concentrated  on  frequencies  close  to  the  resonance  region.  The  ob¬ 
served  enhancement  factors  are  explained6  on  the  basis  of  antenna 
theory.10  Indeed,  for  incident  plane  waves  for  E||  L  orientation,  most 
of  the  observed  results  are  as  though  the  target  acted  like  a  pick-up 
half-wave  dipole  with  reflecting  surfaces  in  close  vicinity. 

Finite  conductivity,  finite  size  reflecting  surfaces  and  other 
orientations  have  not  been  considered  to  date.  Also  the  results  for 
frequencies  higher  than  about  8.5  times  the  resonant  frequency  (550  MHz 
for  man)  have  not  been  obtained  even  for  highly  conducting  reflectors. 

Enclosed  structures  such  as  rooms,  etc.,  may  act  as  lossy 
resonators  with  electromagnetic  fields  coupling  in  from  windows.  If 
such  structures  were  to  have  highly  reflecting  walls,  field  enhance¬ 
ments  by  two  to  three  orders  of  magnitude  may  indeed  be  possible.  How¬ 
ever,  since  typically  encountered  walls  are  not  very  reflecting,  it  is 
expected  that  field  enhancements  by  more  than  a  factor  of  5-10  may  not 
occur  realistically.  Research  into  these  structures  is  needed  in  order 
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to  describe  the  nature  of  field  enhancements. 

Head  Resonance 

We  have  recently  identified  a  frequency  region  for  the  highest 
rate  of  energy  deposition  in  the  head.  The  head  resonance5* 8 » 1 1  (see 
Appendices  C,  E,  and  F)  occurs  at  frequencies  such  that  the  head 
diameter  is  approximately  one  quarter  of  the  free-space  wavelength. 

For  the  intact  (adult)  human  head,  the  resonance  frequency  is  estimated 
to  be  on  the  order  of  350-400  MHz.  At  head  resonance,  the  absorption 
cross  section  for  the  head  region  is  approximately  3.0  times  the  physical 
cross  section  with  a  volume  averaged  SAR  that  is  about  3.3  times  the 
whole-body  averaged  SAR.  Both  values  greatly  exceed  numbers  reported 
earlier12*13  for  spherical  models  of  the  isolated  human  head.  Also 
numerical  calculations11  (Appendix  F)  using  144  cubical  cells  of  various 
sizes  to  fit  the  shape  of  the  human  head  (340  cells  for  the  total  body) 
give  local  SARs  at  hot  spots  that  are  about  5  times  the  head-average 
values. 

We  feel  that  the  phenomenon  of  head  resonance  may  be  important 
in  the  study  of  behavioral  effects,  blood-brain  barrier  permeability, 
cataractogenesis ,  and  microwave  bioeffects. 

Multianimal  Effects5 

It  has  been  shown  (see  Appendix  C  for  details)  that  for  resonant 
biological  bodies  close  to  one  another,  the  antenna  theory  may  be  used 
to  predict  the  modification  in  SAR  relative  to  free-space  values.  For 
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two  resonant  targets  separated  by  0.65  to  0.7  A,  the  highest  SAR,  150 
percent  of  the  free-space  value  has  been  experimentally  observed  for 
anesthetized  rats  and  models  of  man  for  e||  L  orientation  for  frontally 
(broadside)  incident  plane  waves.  For  tnree  animals  in  a  row  with  an 
interanimal  spacing  of  0.65  A,  the  central  animal  SAR  is  observed  to  be 
roughly  two  times,  while  the  two  end  animals  receive  an  SAR  that  is  ap¬ 
proximately  1.5  times  that  for  an  isolated  animal. 

Full  implications  of  the  multibody  effects  on  SAR  are  not 
completely  understood,  even  though  preliminary  experimental  data5  show 
that  similar  enhancements  may  also  occur  for  subresonance  and  supra- 
resonance  regions.  Furthermore,  other  orientations,  irregular  spacings, 
and  nonf ree-space  exposure  conditions  have  not  been  considered  to  date. 

B.  Other  Highlights  of  the  Work  Done  on  the  Project 

Electromagnetic  Absorption  in  a 

Multilayered  Model  of  Man 

A  multilayered  planar  model  is  used  (see  Appendix  G)  to  examine 
the  dependence  of  whole-body  power  absorption  on  the  configuration  of 
surface  layers,  e.g.,  skin,  fat,  muscle,  which  normally  occur  in  biologi¬ 
cal  bodies.  It  is  found  that  the  layering  resonance  for  three-dimensional 
bodies  (as  opposed  to  the  geometrical  resonance)  can  be  predicted  quite 
accurately  by  a  planar  model.  Calculations  for  a  multilayered  prolate 
spheroidal  model  of  man  predict  a  whole-body  layering  resonance  at  1.8 
GHz  with  a  power  absorption  34  percent  greater  than  that  predicted  by 
a  homogeneous  model. 
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Behavioral  and  Thermal  Effects  of  Microwave 
Irradiation  at  Resonant  and  Nonresonant  Wavelengths 

The  details  of  these  studies  are  given  in  Appendix  H. 

Behavioral  and  thermal  effects  of  radiating  an  animal  with  dif¬ 
fering  wavelengths  of  microwave  energy  at  the  same  power  density  were 
investigated  in  the  first  of  two  studies.  Five  Long-Evans  rats  were 
trained  to  perform  a  lever-pressing  task  and  were  rewarded  with  food  on 
a  variable  interval  schedule  of  reinforcement.  Rats  were  individually 

exposed  in  random  order  to  400-,  500-,  600-,  and  700-MHz  CW  radiation 

2 

at  a  power  density  of  20  mW/cm  with  the  long  axis  of  the  rat's  body 
parallel  to  the  vector  of  the  electric  field.  Radiation  at  all  wave¬ 
lengths  produced  rises  of  body  temperature  and  stoppage  of  lever  press¬ 
ing.  The  averaged  rise  in  body  temperature  was  greatest  and  work  stop¬ 
page  was  most  rapid  during  exposures  at  600  MHz.  In  the  second  study, 

six  rats  were  exposed  in  random  order  to  600-MHz  CW  radiation  at  power 

2 

densities  of  5,  7.5,  10,  and  20  mW/cm  while  performing  the  same 

2 

behavioral  task.  Exposures  at  10  and  20  mW/cm  resulted  in  work  stop- 

2 

page,  while  exposures  at  5  and  7.5  mW/cm  did  not.  In  addition,  three 

of  the  rats  were  subsequently  exposed  while  responding  to  600-MHz 

pulsed  radiation  (1000  pps,  3-  or  30-ps  pulse  durations  at  a  peak  power 

2  2 

density  of  170  mW/cm  (averaged  0.51  and  5.1  mW/cm  ).  No  work  stoppage 
occurred  to  pulsed  radiation.  Taken  in  sum,  the  data  show  that  the 
mature  Long-Evans  rat  is  resonant  at  a  frequency  near  600  MHz  while  work 
stoppage  during  short-term  exposures  to  600-MHz  radiation  occurs  at  a 


power  density  between  7.5  and  10  mW/cm 
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Radiators  for  Microwave  Biological  Effects  Research  — 

Waveguide  Slot  Array  with  Constant  Radiation  Intensity 

Microwave  biological  studies  with  large  targets  or  multiple 
targets  require  a  large  working  area  that  provides  near-uniform  exposures 
to  energy.  At  high  microwave  frequencies  where  high  power  sources  are 
not  readily  available,  it  is  desirable  to  use  a  radiator  that  provides 
a  constant  intensity  beam  so  that  the  working  area  may  receive  almost 
all  of  the  radiated  energy. 

Several  synthesis  procedures  have  been  used  for  design  of  a 
longitudinal  shunt  slot  array  having  a  pattern  with  nearly  constant 
radiation  intensity  over  the  beam  width  and  zero  intensity  outside.  A 
beam  width  of  15.5°  was  chosen  as  being  small  enough  for  a  radiated  beam 
to  approximate  a  plane  wave  yet  give  enough  divergence  for  convenient 
chamber  divisions.  Discontinuities  in  the  desired  pattern  force  the 
Fourier-synthesized  pattern  to  have  significant  errors  due  to  overshoot 
and  ripple.  Solutions  developed  on  the  basis  of  Fejer  sums  and  Lanczos 
weighted  least  squares  gave  patterns  with  objectionable  rounding.  A 
new  procedure  based  on  weighted  least  squares  was  developed  for  optimum 
characteristics  and  both  the  new  procedure  and  its  application  to  a 
practical  design  are  described  in  reference  14. 

A  longitudinal-shunt  slot  array  with  reflectors  has  been  designed, 
built,  and  tested.  The  antenna  is  easily  matched  near  the  design  fre¬ 
quencies  of  10.0-10.4  GHz.  The  radiation  pattern  has  an  intensity  within 
±0.5  dB  over  a  total  angle  of  15.6°  in  the  H  plane  and  15.5°  in  the  E 
plane.  Gain  is  19.7  dB  which  is  3.9  dB  above  that  of  a  horn  with  a 
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similar  aperture  area.  By  using  energy-absorbing  material,  it  is  possible 
to  produce  a  pattern  in  which  the  fields  drop  suddenly  from  the  nearly 
constant  values  to  nearly  zero. 

Biological  Phantom  Materials  for  Simulating 

Han  at  Different  Frequencies 

Biological  phantom  materials  were  developed  for  simulating  man 
over  the  frequency  range  13-230  MHz  (including  the  important  resonance 
regions  of  33  and  65  MHz).  These  materials  have  been  used  to  fill  the 
eight  proportionately-scaled  figurines  for  dosimetric  measurements  at 
300,  400,  600,  915,  and  985  MHz.  The  compositions  and  their  measured 
properties  are  given  in  Appendix  A  of  Annual  Report  No.  2. 

Procedures  for  Improving  Convergence  of 

Moment-Method  Solutions  in  Electromagnetics 

The  value  of  numerical  solutions  in  dosimetry  may  be  limited  by 
both  the  accuracy  of  the  model  as  a  representation  of  the  target  animal 
and  the  accuracy  of  the  computations  for  the  model.  Accurate  numerical 
procedures  have  been  developed  for  use  with  crude  models  such  as  the 
sphere,  infinite  cylinder,  ellipsoid,  and  spheroid.  We  have  used 
moment-method  solution  of  the  electric  field  volume  integral  equation 
so  that  an  array  of  cubes  of  different  sizes  and  compositions  may  be 
used  to  better  represent  th>-  shape  of  the  target  animal.  In  order  to 
allow  a  maximum  number  of  cubes  for  increased  detail  in  the  representa¬ 
tion  of  arms,  legs,  etc.,  we  have  used  a  pulse  function  basis  in  the 
moment-method  solution.  Others  have  also  used  pulse  function  solutions 
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for  block  models  in  geophysics  and  in  dosimetry,  but  we  have  addressed 


the  issues  of  both  estimating  and  reducing  the  errors  caused  by  the  use 
of  pulse  functions. 

A  pulse  function  basis  requires  the  assumption  that  the  cells 
are  small  enough  that  the  electric  field  may  be  represented  by  a  constant 
complex  vector  within  each  cell.  Solutions  previously  reported  by 
others  for  a  block  model  of  man  have  a  ratio  of  239:1  for  energy  deposi¬ 
tion  in  one  pair  of  adjacent  cells  at  10  MHz.  We  expect  that  large 
errors  are  present  in  solutions  exhibiting  large  cell-to-cell  variations. 
Our  first  publication  regarding  the  use  of  pulse  functions  is  Appendix  M 
which  establishes  an  upper  bound  on  cell  size.  Such  a  criterion  was  not 
previously  available  for  use  with  moment-method  solutions.  After  defin¬ 
ing  the  source  of  error,  our  next  task  was  to  reduce  the  error.  Ap¬ 
pendix  I  describes  several  novel  procedures  which  we  discovered  for  re¬ 
ducing  the  error  in  solutions  of  2-D  electromagnetics  problems.  We  were 
not  able  to  extend  the  new  procedures  to  3-D  solutions,  but  we  found  a 
different  method  which  is  useful  for  improving  convergence  in  3-D 
problems  which  is  described  in  Appendix  .1.  Armed  with  both  our  con¬ 
vergence  criterion  (Reference  15)  and  a  novel  method  of  improving  con¬ 
vergence  (Appendix  J) ,  we  have  proceeded  to  obtain  numerical  solutions 
for  a  realistic  model  of  man  both  in  free-space  and  under  other  condi¬ 
tions  such  as  in  the  presence  of  a  ground  plane  or  flat  or  90°  corner 
reflectors.  Such  applications  are  described  in  other  publications 
referenced  earlier  in  this  report. 
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Conditions  of  Strongest  Electromagnetic  Power 
Deposition  in  Man  and  Animals 
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Abstract — Strongest  power  deposition  for  biological  bodies  is 
found  for  fields  polarized  along  the  longest  dimension  for  frequencies 
such  that  the  major  length  is  about  0.4  times  the  free-space  wave¬ 
length  of  radiation.  Peak  absorption  in  the  presence  of  ground  effects 
is  observed  at  frequencies  about  one-half  as  much  as  for  bodies 
isolated  in  free  space.  At  resonance,  an  effective  absorption  area  3-4 
times  the  shadow  cross  section  is  measured.  Using  biological-phan¬ 
tom  figurines  of  the  human  body,  distribution  of  power  absorption  is 
determined.  Maximum  power  deposition  is  observed  for  the  neck  re¬ 
gion.  Initial  experiments  with  anesthetized  and  dead  rats  have  con¬ 
firmed  this  observation. 

INTRODUCTION 

KC’Al'SK  of  the  expanding  uses  of  electromagnetic 
(KM)  radiation,  it  has  become  necessary  to  estab¬ 
lish  the  most  hazardous  conditions  of  exposure  to  such 
radiation  for  humans.  In  particular,  it  is  important  to 
know  the  frequency  regions  of  most  absorption  for  various 
orientations  of  the  body  relative  to  the  polarization  of 
incident  waves  and  how  the  absorption  cross  section  varies 
with  frequency  of  KM  radiation.  It  is  also  important  to 
determine  the  distribution  of  power  absorption  in  the  body 
to  ascertain  locations  of  maximum  deposition  particularly 
close  to  resonant  frequencies.  This  knowledge  is  vital  in 
evaluating  and  establishing  radiation  safety  standards. 

Based  on  our  previous  experiments  [1}  with  rats  and 
biological-phantom  and  saline-filled  prolate  spheroidal 
bodies,  strongest  whole-body  absorption  was  found  for 
KM  fields  polarized  with  electric  field  along  the  long 
dimensions  of  the  bodies.  KM  power  absorption  was 
measured  for  bodies  of  fixed  size  at  various  frequencies 
using  a  parallel-plate  waveguide  as  a  transmission  medium 
of  plane  waves  and  compared  [2]  with  data  obtained  from 
fixed-frequency  free-space  irradiation  exposures  of  bodies 
of  varying  sizes.  The  salient  features  of  the  experimentally 
observed  results  (qualitatively  sketched  in  Kig.  I)  areas 
fi  >llows. 

1)  The  E  [|  /,  orientation  is  found  to  be  the  most 
absorbing,  and  ff||  /,  the  least  absorbing,  with  the  con¬ 
figuration  k  ||  /,  otilv  slightly  more  absorbing  than  the 
H  ||  /,  orientation.  The  vectors  E,  H.  and  k  are  along  the 


Frequency  (arbitrarv  units) 


EiK-  1  Typical  RF  absorption  curves  for  animals  and  bodies  of 
prolate-spheroidal  shape  obtained  with  a  parallel-plate  wave¬ 
guide.  Symbols  with  arrows  above  are  represented  by  boldface 
italic  in  text. 

electric  and  magnetic  fields,  and  along  the  direction  of 
propagation,  respectively,  and  /,  is  along  the  major  length 
/,  of  the  body.  These  experimental  results  are  in  quali¬ 
tative  agreement  with  recent  theoretical  results  [.'!]  on 
power  absorption  in  homogeneous  lossy  prolate  spheroids. 
Quantitative  comparison  has  not  been  possible  because  of 
the  long-wavelength  (kb  <  0.(>)  approximation  of  the 
theory.  By  contrast,  the  experimental  results  pertain 
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mostly  to  frequencies  on  the  order  of  resonant  and  higher 
frequencies  ( A/,  >  1.0). 

2 )  The  frequencies  of  peak  absorption  occur  in  the 
reverse  order.  Maximum  absorption  for  E  ||  /,  occurs  at 
the  lowest  frequency  with  A7.  -  2x7, 'X  on  the  order  of 
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TABLE  I 

Actual  and  Equivalent  Prolate  Spheroidal  Dimensions  por 
Anesthetized  Rats  op  Dipperent  Weights 


Weight  of 
the  Ret 

In  Crass 

Ret  Length 

L  from  Hose 
to  Reae  of 

Tall 

cm 

— 

Maaaured 

Circumference 
of  the  Anteal 

Nominal 

Circumference* 

2*b 
in  ema 

A# pact  Ratio 
l/2b  of  Equi¬ 
valent  Prolate 
Spheroid 

96 

15.5 

11.5 

(body) 

MR.  9.0 
(head) 

10.80 

4.50 

158 

18.0 

13.5 

(body) 

MR.  9.0 
(head) 

12.82 

4.10 

281 

21.0 

15.0 

(body) 

MX.  10. 0 

(head) 

13.12 

4.30 

190 

11.5 

18.0 

(body) 

MR.  11.0 
(haad) 

17.70 

4.33 

*  This  is  calculated  from  the  relationship  2/3  trZ.6*p  —  W. 


1.0  to  1.6  (LA  —  0.16  to  0.25).  Peak  absorption  for 
k  1 1  L  and  H\\L  orientations  occurs  at  successively 
higher  frequencies,  with  kL  for  these  configurations  on  the 
order  of  L/2b  where  2b  is  the  length  of  the  equivalent 
prolate  spheroidal  body1  along  the  minor  axes  (see  Table 
I  for  numbers  illustrative  of  90-400-g  rats).  These  values 
of  kL  for  resonance  in  the  different  orientations  are  ob¬ 
tained  approximately  by  requiring  that  for  a  sphere  [4], 
[5]  the  shortest  circumference  for  the  lossy  creeping  waves 
launched  at  the  center  of  the  shadow  plane  be  0.5  X. 

3)  In  the  frequency  region  immediately  below  reson¬ 
ance,  absorption  diminishes  rapidly  with  frequency.  An 
f*-n  dependence  similar  to  that  for  spheres  is  observed. 

4)  There  is  an  excellent  correlation  of  the  results  ob¬ 
tained  from  the  parallel-plate  waveguide  with  the  free- 
space  irradiation  results  [lj,  [2].  The  only  exception  is 
that  maximum  absorption  in  the  free-space  case  occurs  for 
body  size  of  about  0.4  X,  i.e.,  at  approximately  twice  the 
frequency  observed  in  the  parallel-plate  waveguide  situa¬ 
tion.  The  reason  for  this  may  be  that  while  a  body  isolated 
in  free  space  is  required  to  be  approximately  0.4  X  for 
resonance,  the  same  condition  is  met  in  the  presence  of  a 
ground  plane  (such  as  in  the  parallel-plate  waveguide) 
with  only  a  0.2-X  body  and  its  image  in  the  ground  plane 
acting  together. 

CORRELATION  WITH  LETHALITY  RESULTS 

A  study  [6]  has  just  been  completed  at  Walter  Reed 
Army  Institute  of  Research,  Washington,  I).  on  the 
lethality  of  100-125-  and  380-420-g  rats  and  25-35-g  mice 
when  exposed  to  microwave  radiation  at  710,  9X5,  1700, 
2450,  and  3000  MHz.  The  choice  of  frequencies  was 
dictated  by  the  anechoic  chamber  facilities  available  at 
Walter  Reed  Laboratories.  The  time  to  convulsion  was 

1  The  dimension  2b  tor  an  equivalent  prolate  spheroid  is  calcu¬ 
lated  from  the  expression  2/3  rb>Lp  *  W,  where  W  is  the  weight 
of  the  body  and  p  is  the  average  maes  density  of  nearly  1  g/cm*. 
The  nominal  circumference  2t6  so  calculated  has  been  found  to 
correspond  well  to  the  average  circumference  actually  measured  for 
rats  of  weights  varying  from  90  to  400  g  (see  Table  I). 


measured  for  E  ||  L  and  H  ||  L  configurations.  At  each 
frequency,  eleven  animals  were  studied  for  each  of  the  two 
orientations  picked  to  correspond  to  the  maximum  and 
minimum  absorbing  configurations  of  Fig.  1.  The  lethality 
data  correlate  well  with  our  results  [1],  [2]  obtained  from 
the  modeling  experiments.  In  particular,  the  following 
points  of  agreement  are  noteworthy. 

1)  The  time  to  convulsion  in  the  E  ||  L  orientation  is 
alw'ays  shorter  than  that  for  the  H  ||  L  orientation.2 

2)  On  the  low-frequency  side  of  resonance,  the  time  to 
convulsion  rises  very  rapidly,  corresponding  to  a  fast 
diminishing  absorption  in  this  region  (see  Fig.  1 ) . 

3)  The  minimum  time  to  convulsion  for  H  ||  L  orienta¬ 
tion  occurs  for  kL  ~  4.4.  From  Table  I,  the  nominal  L/2b 
for  rats  of  different  weights  is  on  the  order  of  4.4.  The 
lethality  results  are  consequently  in  agreement  with  the 
parallel-plate  waveguide  results  for  resonant  absorption 
at  kL  ^  L/2b. 

4)  For  the  E\\  L  configuration,  minimum  time  to  con¬ 
vulsion  is  obtained  at  frequencies  such  that  the  length 
from  the  animal’s  eyes  to  the  base  of  the  tail  is  approxi¬ 
mately  0.4  X.  This  is  in  agreement  with  the  free-space 
experiments  described  in  the  following  (see  Fig.  2)  and 
also  in  [1]. 

SCALING 

To  evaluate  the  validity  of  scaling,  free-space  heating 

oeriments  were  repeated*  at  710  MHz  for  saline-filled 
prolate  spheroidal  bodies  (L/2b  =  6  corresponding  to 
humans)  of  major  length  varying  from  7.6  to  25.4  cm. 
The  temperature  increase  for  a  5-min  exposure  to  100- 
mW/cm1  radiation  is  plotted  in  Fig.  2.  From  the  amount  of 
energy  absorbed,  the  relative  absorption  coefficient  [7]  S 

1  The  only  exception  to  this  observation  is  400-g  rate  at  2450  and 
3000  MHz.  Por  these  animals,  such  frequencies^  are  4-5  times 
higher  than  the  resonant  frequency  of  the  E  1 1  L  orientation.  A 
considerably  reduced  total  power  deposition  under  these  circum¬ 
stances  may  well  be  the  reason  for  this  anomaly. 

‘Experiments  were  previously  done  at  1700  MHz  with  smaller 
siae  prolate  spheroidal  bodies  with  L  of  4.8,  6,  7.2,  9.6,  and  14.4  cm 
(see  |1|). 
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Fig.  2.  Temperature  rise  in  saline-filled  prolate  spheroids  with  an  aspect  ratio  L/2b  of  6.  Frequency,  710  MHz;  ir¬ 
radiation  time,  5  min;  incident  power  density,  100  mW/cm*;  polarization,  E  ||  L. 
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Fig.  3.  Relative  absorption  coefficient  5  for  different  size  prolate  spheroids  (L/2b  -  6  in  each  case).  Polarization: 


is  calculated*  and  is  plotted  in  Fig.  3.  Maximum  absorp¬ 
tion  is  observed,  as  before  [1],  [2],  for  bodies  of  major 
length  on  the  order  of  0.38-0.48  X,  and  whole-body  deposi¬ 
tion  of  about  3.5  times  as  much  as  that  called  for  by  the 
shadow  area  is  measured.  This  verifies  experimentally  the 
concept  of  scaling.  Bodies  scaled  down  in  all  dimensions 

*  The  relative  absorption  coefficient  S  is  defined  as  the  absorption 
cross  section  divided  oy  the  physical-shadow  cross  section.  Power 
absorbed  in  milliwatts  -  (4180  X  temperature  increase  of  the 
body  X  weight  in  crams) /irradiation  time  in  seconds.  This  divided 
by  the  incident  field  intensity  in  milliwatts/square  centimeter  gives 
the  absorption  cross  section. 


result  in  peak  absorption  frequencies  scaled  up  by  the  same 
factor  and  vice  versa. 

CORE  TEMPERATURE  OF  RATS  EXPOSED 
IN  DIFFERENT  ORIENTATIONS 

In  order  to  dramatise  the  polarisation-dependent 
variability  of  the  EM  hazard,  250-g*  Wistar  rats  were 
exposed  to  100-mW/cm*  fields.  The  rectal  temperature  of 
the  animal  was  measured  under  irradiation,  using  the 

*  This  particular  weight  was  selected  to  observe  the  L  ar  0.4-X 
condition  for  maximum  absorption  in  the  free  space. 
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Fig.  4.  Core  temperature  of  250-g  rats  exposed  to  710-MHz  radiation,  (a)  Orientation  E  1 1  L\  incident  power  density, 
100  mW/cm*.  (b)  Orientation  K  ||  L;  incident  power  density,  100  mW/cm’.  (c)  Orientation  H  ||  L,  incident  power 
density,  100  mW/cm1.  (d)  Orientation  E  ||  L;  incident  power  density,  50  mW/cm'. 


relatively  nonperturbing  liquid-crystal  fiberoptic  tem¬ 
perature  probe  [8].  The  digital  voltmeter  reading  as  a 
function  of  time  was  recorded  and  is  given  in  Fig.  4(a), 
(b),  and  (c)  for  the  three  orientations.  The  core  tempera¬ 
ture  read  off  the  calibration  chart  of  the  probe  is  marked 
alongside  the  curves  in  Fig.  4.  For  an  incident  power 
density  of  50  mW/cm1  [Fig.  4(d)],  the  power  deposition 
in  the  E  ||  L  orientation  was  still  much  higher  than  that 
for  the  other  two  configurations  at  100  mW /cm2.  The  in¬ 
crease  in  animal  temperature  as  a  function  of  time  for  the 
two  power  densities  in  the  E  ||  L  orientation  is  plotted  in 
Fig.  5.  F  or  this  configuration,  the  relative  absorption 


coefficients  S  of  3.8  and  4.26  are  calculated  from  the 
initial  slopes  of  the  temperature-rise  curves  for  100-  and 
50-mW/cm2  incident  field  intensities,  respectively.  This 
may  be  compared  to  <S  on  the  order  of  3.5  for  resonant-size 
saline-filled  prolate  spheroids  ( Fig.  3). 

DISTRIBUTION  OF  POWER  DEPOSITION— 
“HOT  SPOTS” 

To  determine  the  distribution  of  power  deposition,  an 
18.4-cm-tall  figurine  of  the  female  body  has  been  used  as  a 
mold  to  form  a  cavity  in  styrofoam.  The  size  of  the  figurine 
has  been  selected  such  as  to  allow  resonance  (Fig.  3)  at 
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Irradiation  tine  In  Minutes 

Fig.  5.  Core  temperature  increase  of  250-g  Wistar  rats  exposed  to 
710-MHz  free-space  radiation  (E  ||  /.  orientation). 


Fig.  6,  Photograph  of  the  parallel-plate  radiation  chamber. 


710  MHz  for  free-space  radiation  ( /,/A  ~  0.44 )  and  ;.t  a 
frequency  about  half  as  much  in  the  parallel-plate  radia¬ 
tion  chamber.  In  order  to  simulate  ground  effects  on  KM 
power  absorption,  a  larger  parallel-plate  radiation  chamber 
(I'ig.  (i)  has  been  designed  and  fabricated.  For  the 
chamber,  the  central  working  area6  consists  of  a  copper 
plate  of  63.5-em  width  separated  bv  a  2f>.4-em  clearance 
from  the  ground  plate  of  1 17-cm  width.  The  overall  trans¬ 
mission  length  of  lil.S  cm  is  occupied  by  two  symmetrical 
tapered  sections  of  fit -cm  axial  length,  the  central  working 
area  of  fil-om  length,  and  the  7..r)-em  end  sections  con¬ 
nected  to  UG/iX/lJ  input  and  output  coaxial  connectors. 

From  KM  field  t  heory,  a  body  reduced  by  a  factor  of  0  in 
all  dimensions  may  bo  used  to  obtain  RF  absorption 


*  the  uniformity  of  fields  in  the  central  working  area  was  <wt.ab- 
lished^by  heating  a  19-cm-long  saline-filled  prolate  spheroid  in  the 
E  ||  1j  orientation.  For  the  spheroid  placed  at  the  center  of  the 
plates,  a  temperature  increase  of  1  23°C  wns  measured  for  an  inci¬ 
dent  field  intensity  of  49.lt)  mW/cm*  at  360  MHz.  By  comparison, 
a  temperature  increase  of  II. TO  at  43. SX  mW/crn'  was  observed 
for  the  same  spheroid  placed  10.2  cm  off  the  center  line 


characteristics  of  the  full-scale  body,  provided  an  ir¬ 
radiation  frequency  scaled  up  by  a  factor  of  0  is  used. 
However,  it  is  necessary  that  the  complex  permittivity 
(t,  —  ja/wto)  represented  by  the  scaled-down  model 
correspond  to  tin1  value  at  the  lower  frequency  charac¬ 
teristic  of  the  whole  body.  With  these  scaling  precautions, 
the  distribution  of  power  deposition  in  the  reduced  sc  .no 
model  is  identical  (though  tin1  magnitudes  are  higher  by  a 
factor  of  0)  to  that  of  the  full-scale  body.  Since*  the  pro¬ 
jected  frequency  of  maximum  absorption  for  adult 
humans  is  .,n  the  order  of  65-75  MHz  in  free  space  and 
about  half  as  much  in  the  presence  of  ground  effects,  the 
biological-phantom  material  that  is  used  should  have  the 
dielectric  properties  of  humans  in  that  frequency  range. 
A  composition  [9]  of  the  material  with  3.26  percent  NaCl, 
S.74  percent  Superstuff  (obtained  from  Whamo  Manu¬ 
facturing  Company,  San  Gabriel,  California),  and  87.0 
percent  water  has  a  measured  dielectric  constant  of  60  and 
a  conductivity  of  4.39  mho/m.  For  use  at  a  frequency  of 
3.r>0  MHz,  this  corresponds  to  tr  -  66  —  /225.6.  By  com¬ 
parison,  for  muscle,  skin,  and  tissues  with  high  water 
content,  t,  =  XX  —  ,/250  at  .30  MHz  [4].  Because  of  the 
relative  closeness  of  the  two  dielectric  constants,  the 
mixture  of  this  composition  was  used  to  fill  the  figurine- 
shaped  cavity  for  experiments  in  the  parallel-plate  radia¬ 
tion  chamber.7  Temperature  under  irradiation  was  re¬ 
corded  using  the  liquid-crystal  temperature  probe  [s']. 
From  the  increase  in  temperature,  the  absorbed  power 
density  in  milliwatts/gram  was  calculated  from  the 
relationship  [7]  (4180  X  temperature  increase) /irradia¬ 
tion  t  ime  in  seconds.  The  absorption  coefficient  a.8  defined 
as  milliwatts/gram  of  absorbed  power/milliwatts/square 
centimeter  of  incident  field  intensity  for  various  parts  of 
the  body,  has  been  calculated  and  is  given  in  Table  1 1 .  For 
E  ||  Ij  orientation,  the  strongest  intensity  of  power  deposi¬ 
tion  is  observed  in  the  neck  area  of  the  body.  For  the 
density  of  power  absorption  this  is  followed  by  the  shins, 
the  thighs,  the  chest,  the  eyes,  and  the  pudendal  region,  in 
that  order.  In  Fig.  7  the  thermographic  record  of  the 
temperature  before  and  after  1  min  of  free-spare  irradia¬ 
tion  with  a  field  intensity  of  100  mW/cm2  at  710  MHz  is 
shown.  1'he  pattern  of  power  deposition  bears  a  remarkable 
resemblance  to  the  distribution  obtained  in  the  parallel- 
plate  radiation  chamber. 

In  order  to  obtain  a  quantitative  information  on  the 
power  deposition  under  free-space  irradiation  conditions, 
a  12.1-cm  resonant-size  (L/X  ~0.4)  figurine  was  used  at 
9.X.3  MHz.2  A  field  intensity  of  100  mW/cm2  was  used  to 

7  For  free-space  irradiation  experiments,  the  composition  of  the 
material  used  is  9.25  percent  NaCl,  R.25  percent  Superstuff,  and 
32.5  percent  water.  This  mixture,  having  a  measured  dielectric 
constant  of  66  and  a  conductivity  of  10.09  mho/m,  corresponds  to 
a  c  *mplex  relative  permittivity  o  of  66  -  /255  when  used  at  710 
MHz. 

•  The  parameters  X  and  a  are  related  in  that  X  =  (o  X  weight  of 
the  body) /(shadow  area  of  the  body).  For  animals  and  other  bodies 
of  irregular  shape,  .S  is  a  more  difficult  parameter  to  calculate  Ire- 
cause  of  the  difficulty  in  evaluating  the  shadow  area  of  these  bodies. 

9  K.  I,.  Hunt  of  the  Department  of  Microwave  Research,  Walter 
Reed  Army  Institute  of  Research,  Washington,  D.  C.,  participated 
in  these  experiments. 
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TABLE  II 

Measured  Distributionof  Power  Deposition  in  an  18.4-cm-Tall 
Biologic  a  l-Ph a ntom  Human  Figurine  (Parallel-Plate 
Waveguide  Radiation  Chamber) 


Frequency 

a 

\~T  ~ 

Pudendal 

HHz 

l./» 

Neck 

Shin 

r„. 

{  Region 

|||  L  orientation 

230 

0.1*1 

2.895 

2*0 

0.1*2 

3.381 

0.788 

0.135 

220 

0-  16b 

1.927 

300 

0.18* 

3.097 

2.2  56 

0.8*0 

0.4*2 

320 

0.146 

3.367 

330 

0.213 

6.761 

1.86  3 

180 

0.233 

4.  16  8 

1.249 

19$ 

0.2*2 

6.133 

k ||  1.  orientation 

393 

0.2*2 

0.09* 

fi||  L  orientation 

300 

0. 18* 

0.130 

Note:  Symbols  with  arrows  above  are  represented  by  boldface 
italic  in  text. 


(b) 

Fig.  7.  The  thermographic  recording  of  the  central  section  of  an 
18.4-cm  biological-phantom  figurine  before  and  after  free-spaee 
plane-wave  irradiation  for  1  min  at  100  mW/cm*  at  710  MHr, 
(LA  ~  0.44).  (a)  Before  exposure  to  radiation,  a8.5°C  baseline; 
3°C  full  scale,  (b)  After  irradiation  in  the  E  ||  L  configuration ; 
waves  incident  on  the  front  side  of  the  body,  6°C  baseline;  10°C 
full  scale. 

sptH'd  up  the  heating  process  and  thus  reduce  the  heat 
diffusion  to  cooler  areas  of  the  body.  Temperature  under 
irradiation  was  recorded  using  the  liquid-crystal  tempera¬ 
ture  probe.  Due  to  the  extremely  large  energy  deposition 


rate  for  the  neck  region,  exposure  times  on  the  order  of  .r>  s 
were  used  when  measuring  temperature  increases  for  this 
part  of  the  body.  Irradiation  times  of  10-20  s  were  used 
for  measurements  in  other  parts  of  the  body.  The  measured 
rates  of  temperature  rise  are  used  to  calculate  the  a’s  for 
various  parts  of  the  body  in  a  1.75-m  man  and  these  are 
given  in  Table  III.  The  patterns  of  power  deposition  ob¬ 
tained  from  both  the  parallel-plate  waveguide  radiation 
chamber  and  free-spaee  irradiation  are  very  similar.  For 
the  E  ||  L  resonance  condition,  the  very  maximum  power 
deposition  is  observed  for  the  neck  region  of  the  body. 
Under  both  free-spaee  and  parallel-plate  waveguide 
irradiation  conditions,  an  a  parameter  of  approximately 
20  times  the  whole-body  average  is  observed  for  the  neck 
region  of  the  body. 

A  point  of  concern  is  the  much  lower  values  observed  in 
the  parallel-plate  radiation  chamber  as  compared  to  the 
values  measured  under  free-spaee  irradiation.  This  may  be 
seen  by  comparing  the  values  in  Table  II  to  those  of 
column  .3  of  Table  III.  The  fact  that  a  larger  scaling  factor 
0  is  used  for  the  12.1 -cm  figurine  than  for  the  l<S.4-em 
figurine  should  account  only  for  a  factor  of  l.o2  in  going 
from  Table  II  to  Table  III.  Somewhat  reduced  value's  in 
Table  II  will  also  be  expected  because  of  lower  field  in¬ 
tensities  and  the  resultant  heat  diffusion,  particularly 
from  hot  spots.  In  order  to  create  a  more  satisfactory 
representation  of  the  real-life  ground  effects,  a  monopole- 
above-ground  radiation  has  been  designed  and  is  currently 
under  construction  at  our  laboratory.  The  new  radiation 
chamber  ii>  sketched  in  Fig.  8.  The  radiator  consists  of  a 
45°  corner  reflector  in  conjunction  with  a  quarter-wave 
monopole  above  ground.  The  gain  of  such  an  antenna  has 
been  calculated  [10]  and  is  tabulated  in  Table  IV.  For  the 
antenna-to-corner  spacing  of  0.6.r>  X,  the  calculated  radia¬ 
tion  pattern  of  the  antenna  system  is  plotted  in  Fig.  9.  The 
field  intensity  at  a  far-field  distance  of  fi  ft  for  3.50  W  of 
radiated  power  may  be  as  high  as  60  mW/cm2.  Such  field 
intensities  are  adequate  for  the  planned  experiments.  The 
radiation  chamber  is  being  constructed  with  the  provision 
for  alterations  for  use  at  different  frequencies  in  the  300- 
1000-MHz  band.  Such  a  radiation  chamber  is  not  easily 
amenable  to  wide-band  whole-body  absorption  measure¬ 
ments  such  as  have  been  possible  in  t  hi'  parallel-plate 
t  ransmi ssion  -line  «si ng  a  network  analyzer.  The  monopole- 
above-ground  radiation  chamber  will,  however,  offer  a 
real-life  simulation  of  ground  effects  on  EM  power  absorp¬ 
tion  in  man  and  animals. 

A  significant  result  of  the  experiments  with  biological- 
phantom  figurines  is  that  a  high  power  deposition  in  the 
neck  region  is  observed.  Power  absorption  coefficient  a  = 
6.76  for  the  neck  region  for  E  ||  /,  orientation  is  over  17 
times  larger  than  the  average  value  for  saline-filled  pro¬ 
late  spheroids  of  the  same  overall  dimension  (see  Table  V 
for  whole  body  or  average  a  for  a  19-cm  saline-filled  pro¬ 
late  spheroid  at  different  frequencies).  For  k  ||  /,  and 
H  ||  orientations,  the  rate  of  heating  is  minimal  (Table 
II)  and  this,  for  the  neck  region,  is  a  factor  of  40  to  ,r>0 
times  smaller  than  the  rate  for  the  E  ||  L  configuration. 
For  these  orientations  the  results  once  again  are  in  agrec- 
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TABLE  III 

Rates  of  Energy  Deposition  in  Various  Parts  of  the  Body 
Measured  Data  and  Projected  Rates  of  Energy  Deposition  in  a  1.75-m  Man 


Region  of 
tha  tody 

985  MHz  Resonant  Size  (12.1  c»;  L  *  0.4  ») 
Bio  logical- Phan to«  Figurine 

1.75 

Kan,  67.9  KHz 

Rata  of  Teapereture  , 
Incraaaa  at  100  mV/ at1 
•C/Mln. 

(W/kgal/faU/ca2^ 

(V/kga)/{  aW/ca2) 

Rate  of  Energy  Deposit  Ion 

For  Incident 
Fields  of 

10  aW/ca2 
W/kga 

For  Incident 
Field#  of 

50  aW/ca2 
W/kga 

Heck: 

] 

Upper  part 

128.3 

89.38 

6.16 

61.6? 

308. 10 

Lcwer  part* 

65.2 

45.42 

3.13 

31.31 

156.57 

Shin  a 

9.5 

6.62 

0.66 

6.56 

22.80 

Thighs 

10.8 

7.52 

0.52 

5.19 

25.95 

Cheat 

8.64 

6.02 

0.41 

6.15 

20.75 

Pudenda  1 

Ration 

3.2 

2.23 

0.15 

1.54 

7.7 

lyae 

0.9 

0.63 

0.043 

0.43 

2.15 

Ana 

? 

* 

? 

_ _ _ 

Note:  E  ||  L  resonance  condition  (L/\  =  0.4). 

*  The  rate  of  temperature  rise  for  the  lower  part  of  the  neck  is  considerably  lower  because  of  its  vicinity  close  to 
the  torso  which  provides  a  large  heat  drain. 


-10  *SB 


Fig.  9.  Calculated  radiation  pattern  at  the  ground  plane  of  a  mono¬ 
pole  in  a  45°  corner  reflector. 


TABLE  IV 

Calculated  Field  Intensities  Due  to  a  Quarter-Wave 
Monopole  in  a  45°  Corner  Reflector 


Dlatance  to  the 
Corner  of  the 
Reflector 

Driving  Fewer 
lapedance 

Calculated  Cain 
of  the  Antenna 

Field  Intensity 
at  R  -  6  feet 

a 

I 

Ohae 

aW/ca2 

0.55 

16  *  J56.5 

75 

62.4 

0.6 

23.5  ♦  J64.6 

79.3 

66.0 

0.65 

32  ♦  165.0 

82.7 

66.9 

0.7 

47.5  ♦  j61.4 

73. 1 

60.9 

0.75 

59.5  ♦  J59.7 

70.87 

59.0 

Note:  Transmitter  Power  *  350  W. 
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TABLE  V 

Absorption  Parameter  a  or  a  1!)-cm  Saline-Filled  Prolate 
Spheroid  as  a  Fcnction  or  Kreqcency 


Frequency 

MHz 

Number  of 

L/\ 

(average ) 

240 

1 

0.132 

.2943 

260 

3 

0.163 

.  3097 

\m 

\ 

0.V90 

320 

2 

0.203 

.3214 

360 

3 

0.229 

.4089 

380 

4 

0.241 

.4908 

393 

4 

0.231 

.  3786 

(b) 


Fig.  10.  Thermographic  recording  of  the  central  section  of  an  18.4- 
cm  biological-phantom  figurine;  free-space  irradiation  of  710 
MHz  (temperature  of  the  body  before  irradiation  ■—  7  J°C).  (a) 
After  irradiation  in  the  k  ||  L  configuration  for  10  min  at  100 
mW/cms;  baseline  8°C;  3°C  full  scale,  (b)  After  irradiation  in  the 
H  ||  /-configuration  for  10  min  at  100  mW/cm';  baseline  8°C; 
3°C  full  scale. 

merit  with  these  obtained  from  free-space  irardiation.  For 
up  to  10  min  of  exposure  to  710-MHz  free-space  fields  of 
intensity  100  mW/ein*,  negligible  modifications  of  the 
temperature  pattern  were  observed,  and  this  is  illustrated 
in  Fig.  10. 

The  thermographic  record  of  the  temperature  before 
and  after  1  min  of  free-space  irradiation  (£  ||  /-  orienta¬ 
tion)  at  2450  Mllz  is  shown  in  Fig.  1 1.  Maximum  heating 
at  2450  MHz  (7,/X  1.5),  approximately  4.75  times  the 

resonant  frequency,  is  still  observed  in  the  neck,  the  arms, 
and  the  legs.  As  expected,  the  temperature  rise  is  con¬ 
siderably  lower  than  that  fur  710-MHz  radiation.  Also,  the 
torso  is  heated  hardly  at  all,  which  is  a  feature  distinct  from 
the  power  distribution  at  the  resonant  frequency  of  710 
M  Hz  (  Fig.  7 ) .  The  quantitative  values  of  a  measured  with 


(b) 


Fig.  11.  The  thermographic  recordings  of  the  central  section  of  an 
18.4-cm  biological-phantom  figurine  before  anil  after  free-space 
plane-wave  irradiation  for  1  min  at  100  m\V/rm!  at  2450  MHz 
(/,/A  =  1.50).  (a)  Before  exposure  to  radiation,  H.0°C  baseline, 
1°C  full  scale,  (b)  After  irradiation  in  the  E  ||  )■  configuration; 
waves  incident  on  the  front  side  of  the  body;  1 0°C  black  level; 
1°C  full  scale. 

the  liquid-crystal  temperature  probe  for  a  somewhat 
smaller  size  figurine  are  given  in  Table  VI.  The  hottest 
parts  of  the  body  are  tin1  neck,  the  legs,  and  the  arms. 

The  large  power  deposition  observed  in  the  neck  region 
of  the  figurines  has  prompted  us  to  look  for  similar  effects 
in  animals.  Initial  measurements10  conducted  with  the 
liquid-crystal  probe  implanted  subcutaneously  in  the  neck 
region  of  anesthetized  5XX-g  long  Evans  rats  have  shown 
a  =  0.479  at  500  MHz  as  against  a  of  0.5X  for  the  deep 
rectal  region.  For  dead  animals,  on  the  other  hand,  the 
corresponding  numbers  are  a  =  0.X55  for  the  neck  region 
and  a  =  0.251  for  the  core  of  the  animal.  Further  experi¬ 
ments  arc  currently  in  progress  to  ascertain  the  nature  of 
neck  heating  and  its  biological  implications.  Since  reports 
in  the  literature  suggest  that  continued  low-level  exposure 
to  microwave  irradiation  in  both  humans  [11]  and  ani¬ 
mals  [12]  produces  fatigue  and  hypoaetivity,  the  possi¬ 
bility  exists  that  this  may  be  linked  to  the  heating  in  the 
neck  region.  The  heating  in  the  neck  area  may  produce 
changes  in  circulating  blood  hormones  such  as  thyroxine. 
In  future  experiments,  blood  thyroxine  levels  will  be 
examined  by  serum  thyroxine  and  protein-bound  iodine- 
level  assay  techniques. 

SUMMARY 

The  UF  power  deposition  is  found  to  vary  significantly 
with  orientation  and  with  frequency.  The  strongest  ab¬ 
sorption  is  found  for  waves  polarized  along  the  long  di- 

10  Threw  experiments  arc  being  performed  by  J.  D’Andrea  of  our 
laboratory. 
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.  TABLE  VI 

E  ]|  /.,  Off-Resonance  Condition  (LA  =  0.0X5) 

Measi  red  Data  and  Projected  Rates  of  Energy  Deposition  in  a  1.75-m 


2450  HHi 

12.1  <•»>  Blol-gi.  ji-Ptuntar  Flgut 


hr 


Region  of  Rate  of  Temperature 


’  <W/kga)/Uw/c»  ' 


1 


J.*40 
I.  20 

!  - 


I  Thighs 
:  Chest 


0.425 

0.2 


0.0f> 
0  02 
0.01 


P.i.I  in-  iJ.nl  T  t'T  lii.-ld.i.l 


J  h'i  i  11. -1 

J  .'1  i  II. ID 


0.  SH 
o.:i 


Man 


inctisiun  of  bodies  at  frequencies  such  that  the  major 
length  is  approximately  0.4  times  tin1  wavelength  of 
radiation  for  bodies  in  free  spare.  Peak  absorption,  in  the 
presence  of  ground  effects,  is  observed  at  frequencies  about 
one-half  as  much  as  bodies  isolated  in  free  space.  A  sig¬ 
nificant  result  of  the  experiments  with  biological-phantom 
figurines  is  that  a  high  rate  of  energy  deposition  is  ob¬ 
served  for  the  neck  region.  Power  absorption  coefficient  a 
for  the  neck  region  for  E\  \  L  orientation  is  about  20  times 
larger  than  the  whole-body  average.  This  may  be  bio¬ 
logically  significant  and  is  being  investigated  further. 
Initial  experiments  with  anesthetized  and  dead  rats  have 
indeed  demonstrated  a  higher  power  deposition  in  the  neck 
region  than  in  the  core  of  the  animals.  For  figurines  for 
k\\L  and  H  \\  1.  orientations,  the  rate  of  heating  is 
minimal,  and  for  the  neck  region  is  a  factor  of  40  to  50 
times  smaller  than  the  rate  for  the  E  '!  /.  resonance  con¬ 
dition.  The  frequencies  of  maximum  absorption  for  adult 
humans  is  projected  to  be  on  the  order  of  05-75  MHz  for 
free-space  conditions  and  about  half  as  much  in  the 
presence  of  ground  effects. 
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Abstract — Numerical  calculations  of  absorbed  energy  deposition  have 
been  made  for  a  Mock  model  of  man  that  is  defined  with  careful  attention 
given  to  the  biometric  and  anatomical  features  of  a  human  being.  <  alcu 
lated  post-resonant  absorption  and  distribution  of  energy  deposition 
through  the  body  have  better  agreement  with  experimental  results  than 
previous  calculations  made  using  less  realistic  models. 

Manuscript  received  September  18.  1978,  revised  April  9,  1979 
The  authors  are  with  the  Departments  of  I  leuncal  Engineering  and 
Bioengineering,  University  of  1  lab  Sail  I.ake  <  its.  I "I  841 12 


I  IsiHoin  chon 

I  HAVE  used  numerical  methods  with  a  realistic 
model  lo  calculate  the  deposition  of  electromag¬ 
netic  energy  in  man.  A  block  model  of  man  was  chosen  to 
allow  maximum  freedom  in  defining  both  the  shape  and 
content  of  the  model.  We  have  used  the  moment-method 
solution  of  the  electric-field  integral  equation  with  a  pulse 
function  basis  and  delta  functions  for  testing  |l),  1 2 ).  The 
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Fig  I  Improved  model  of  man 


following  sections  of  this  paper  contain  brief  descriptions 
of  the  improved  model  and  numerieal  method',  .i'  weil  a' 
a  comparison  of  the  calculations  with  othet  niimencal  and 
experimental  results. 

II  I)i  si  k ip i tt is.  ni  mi  Modi  i 

We  have  given  careful  attention  to  hiomeinc  and 
anatomical  diagrams  in  defining  the  shape  and  content  of 
a  model  having  improved  realism  \  total  ol  ISO  cubical 
cells  of  various  si/es  was  used  to  obtain  a  best  fn  of  the 
contour  on  diagrams  of  the  50th  percentile  standard  man 
[3|.  Sizes  and  placement  of  the  cell'  aie  shown  in  lag.  I 
Height  of  the  model  is  I  7<;  in 

Partitioning  of  bone  fat.  skin.  muscle.  lung  tissue  air. 
heart,  brain,  kidnev  liver  and  spleen  through  the  cells 
was  done  with  the  aid  of  anatomical  cross  sections  [4|.  |5|. 
Whole-hods  volume  fractions  of  each  tissue  tvpe  are 
agreement  with  published  values  I  he  volume-weighted 
complex  permittivitv  of  each  cell  was  found  using  proper¬ 
ties  reported  in  the  literature  |('|  |.'| 

We  have  forced  a  plane  of  svmmetrv  between  the  left 
and  right  bodv  halves  in  selecting  ihe  size,  pkuement.  and 
volume-weighted  complex  permittivities  of  the  cells  I  In- 
symmetry  has  been  utilized  so  th.it  the  size  of  the  matrix 
used  in  the  numerical  solutions  is  dcicimnud  In  mi  cells 
rather  than  the  total  of  lf<<>  cells  Ihcsvmniciic  mav  onlv 


be  utilized  for  incident  fields  m  vchich  both  /.  and  A  are 
con'amed  in  the  plane  of  svmmetrv. 


Ill  Stli  L  IKIN  lit  mi  M-virix  l-.lg  MIUN 

Since  one  plane  of  svmmetrv  is  used  and  there  are  three 
unknown  field  components  in  each  cell,  the  matrix  is  270 
hv  2"’0  complex,  which  corresponds  to  a  system  of  540 
simultaneous  equations  m  540  real  unknowns.  Solution  of 
the  matrix  equation  is  complicated  by  the  matrix  being 
both  full  and  asv mnielric  t  tilizing  svmmetrv.  or  cells 
iiaving  different  size'  ot  pe:  nattiv  ities.  would  each  be 
suflicient  to  destroy  svmmetrv  of  the  matrix 

We  have  used  the  successive  ovc-rreluxution  (SOR> 
method  |4]  for  iterative  solution  of  the  matrix  equation  It 
is  known  that  a  matnx  having  at  least  one  complex 
eigenvalue,  as  m  the  present  case,  will  not  allow  conver¬ 
gence  with  SOR  unless  all  eigenvalues  have  real  parts  with 
magnitudes  lews  than  unity  Convergence  ot  SOR  ha'  the 
further  restriction  tii.it  qlc  H-j.ixation  paianiete:  must  be 
less  than  a  critical  value  that  depends  upon  the  matrix  |4|. 
We  have  found  that  coiHergence  with  SOR  requires  a 
relaxation  factor  Ic-"  than  .intv  at  lit  MHz.  and  the 
critical  value  decreases  with  .tic reusing  frequenev  until  it 
reaches  zero  near  vn  MHz  <  onvercence  with  SOR  is  not 
possible  at  frequencies  exceeding  40  MIC  wilri  tile  present 
model  of  man.  so  we  nave  been  forced  to  use  Cruuss 
Jordan  and  other  noniterative  method'  j  loj  In  the  prc'ent 
problem  when  SOR  mav  be  used  in  piucc  of  noniterative 
methods,  the  computer  time  is  tvpicallv  reduced  hv  about 
a  factor  of  four. 

We  have  oidcrecl  numbering  o‘  .c-lis  m  the  mode1  in 
such  a  wav  that  the  cliffeiencc  between  the  number'  for 
two  cells  i'  .in  me rcasine  function  ot  the  di'tunce  between 
the  cells.  We  have  foimcd  the  matrix  using  an  arrav  of 
3  >  3  blocks  for  each  cell  pair  which  is  inside  out  from  the 
structure  used  previously  1 1  1 1.  Ihe  two  changes  just  de¬ 
scribed  have  caused  the  matrix  to  have  a  rapid  decrease  in 
magnitude  of  the  element'  with  increasing  distance  from 
the  diagonal.  I -or  all  f icqucnc ics  used  so  far  with  the 
realistic  model  of  man  (10  not)  Milzi.  the  onlv  matrix 
elements  having  appreciable  magnitude  are  contained  III  a 
band  about  the  diagonal,  which  amounts  to  approxi¬ 
mately  10  percent  of  the  entire  matrix.  Ihe  approximation 
of  forced  banding  mav  be  applied  to  such  a  matrix  to 
allow  either  a  large  reduction  in  cost  of  computations  or  a 
significant  increase  in  number  of  cells  [12].  We  have  not 
yet  implemented  fotccd  banding,  but  we  estimate  that  the 
method  will  allow  enough  cells  so  that  reliable  calcula¬ 
tions  may  be  t"  icle  foi  a  realistic  model  of  man  to 
frequencies  as  high  a'  24l'i)  MHz 

l\  (iiNVIKi.IMI  (  llNSIDl  K  MIONS 

Numerical  solutions  using  a  pulse  function  basis  can 
or.lv  be  ex  iv  l  ;!  there  is  no  variation  ol  /  within  each  cell 
ol  t he  mode.  Solutions  home  huge  diflcrenccs  in  calcu¬ 
lated  cruiL-v  deposition  toi  .ni|. icent  cell'  must  have  signi- 
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ficam  variation  of  l  within  the  cells,  which  suggests  lack 
of  convergence.  Previously  reported  solutions  for  block 
models  of  man  have  a  ratio  of  234  :  1  for  energv  deposi¬ 
tion  :n  one  pair  of  adjacent  cells  at  10  MHz.  |1|.  The 
arrangement  and  different  sizes  of  cells  in  the  present 
model  cause  the  maximal  ratio  of  energv  deposition  for  a 
pair  of  adjacent  cells  to  he  8  .  1  at  the  same  frequency. 

If  A  is  the  side  of  a  cubical  sell  and  A  represents  the 
magnitude  of  the  complex  propagation  vector  within  the 
cell,  then  substantial  sanation  of  the  electric  field  must 
occur  wninn  the  cell  if  AA  •  \  h  [13]  I  hits  the  pulse 
function  approximation  mas  not  he  justified  at  frequen 
oes  above  200  MHz  with  the  new  model,  suggesting  an 
increasing  ertor  at  higher  frequencies.  Volume-weighting 
of  the  complex  permittivity  within  each  cell  mas  also  he 
unjustified  at  frequencies  above  200  MHz  though  the 
results  app  at  reasonable  to  at  least  500  MHz. 

A  numerical  solution  using  a  pulse  function  basis  re¬ 
sults  in  a  single  value  representing  /  within  each  cell.  It  is 
possible  to  use  the  /.  values  to  calculate  l/2n/5£*  for 
each  cell,  and  to  use  a  volume  average  to  estimate  the 
specific  absorption  rate  (.SAR  I  l.arge  number s  of  cells 
must  be  used  in  order  to  find  accurate  values  of  SAK  bv 
such  a  procedure.  1  or  example,  our  calculations  of  the 
S\R  a  12-cm  muscle  cube  at  I  MHz  show  an  error  of 
3"  percent  with  X  cells.  2h  percent  with  23  cells,  and  20 
percent  with  M  cells.  I  he  delta  functions  used  for  testing 
enforce  the  integral  equation  at  the  center  of  each  cell  so 
that  the  calculated  t.  '  allies  are  most  representative  of  the 
cell  centers  Inspection  of  the  solutions  suggests  th.it  the 
I.  cat  value-  of  /  have  appreciably  less  error  than  occurs 
in  the  SAK 

l!  there  is  much  sanation  ot  /  within  a  scattering  body, 
then  csen  if  we  had  exact  values  of  /.  at  the  cell  ccnteis. 
appreciable  error  would  be  expected  in  the  calculated 
SAK  Ihe  3.  suppo-edl’.  exact,  local  values  of  energv 
deposition  constitute  samples  of  an  otherwise  unknown 
population  so  that  if  large  eell-lo-eell  variation  is  present, 
statistical  procedures  may  he  used  to  calculate  confidence 
limits  for  values  of  whole-buds  or  part-body  SAR  ( 1 4 j 
I  or  man  at  8b  MHz  with  /.  /  and  A  front-to-back.  the 
HD-percent  confidence  limit-  for  whole-hods  SAR  are 
■  12  2  percent  >'t  the  mean  for  the  realistic  model  of  man. 
and  i  143)  percent  of  the  mean  for  ealulalions  reported 
earlier  having  somewhat  create;  eell-to-eell  variation  [1], 

We  have  found  that  accuracy  is  improved  bv  using  a 
three-dimensional  mterpol.mt  with  the  /.  values  initially 
calculated  for  each  -ell  to  account  for  some  of  the  varia¬ 
tion  of  f.  within  each  cell  1 1  5 1  1  rthnear  and  triquadratic 

mterpi ilant-  have  both  been  used  to  estimate  the  variation 
i  f  /■  between  'he  ccii  center-  Ihe  mterpol.mt  is  integrated 
m  calculating  the  SAR  I  or  the  12  c  in  muscle  cube  at  I 
MHz  the  SAR  found  with  2  :  cells  and  the  biquadratic 
tnterpol.ini  have  ciior  comparable  with  the  calculation 
ba-ed  on  M  ceil-  without  the  mtetpolant  Ihe  inciease  in 
cod  due  to  use  o'  ihe  interpolant  is  about  I  percent  Ml 


values  of  energy  absorption  in  tin-  paper  have  been 
calculated  using  interpolant-. 

\  O  IMIVXRIM  IN  OI  (All  I  I  Min-'  w  |  |  II  ( )  I  III  K 
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big.  2  show-  tile  values  ol  whole-bods  S.AR  calculated 
for  man  in  free  space  with  /.  /  and  an  incident  intensity 
of  1  inW  cm'  Points  on  lag  2  represent  the  results  of 
numerical  solutions  fot  the  realist, c  model  of  man.  both 
with  and  without  the  use  of  an  interpolant.  Ihe  homoge¬ 
neous  approximation  of  two-thirds  the  complex  permittiv¬ 
ity  of  muscle  u.i-  used  in  all  calc  illations  fm  the  figure  It 
is  seen  in  big.  2  that  the  fractional  correction  which  is 
made  by  the  interpolant  increases  at  high  Irequencies 
where  increasing  error  t-  expected. 

Values  of  whole-body  SAR  below  resonance  in  1  ig  2 
are  typically  wr’un  ID  percent  of  values  cal- dialed  lot 
prolate-spheroidal  and  ellipsoidal  models  using  the  -urnc 
homogeneous  approximation  of  complex  permittivity  [by1 
hxpernnenta!  results  have  shown  that  the  SAR  lot  /. 
polarization  is  inveiseiy  proportional  to  frequency  tor 
frequencies  from  ahou"  I  to  n  times  that  of  resonance  [i'l. 
[18]  bor  small  animals,  such  as  the  rat.  the  high  value  ol 
the  re  man!  frequency  allow-  a  reduced  magnitude  lot 
the  relative  permittivity  of  tissue  so  that  post-resonant 
calculation-  for  prolate-spheroidal  models  ma\  be-  nude 
using  the  extended  boundary  condition  method  (  H>|.  Such 
results  show  the  experimentally  observed  1  f  hehaviot. 
but  the  calculations  there  cannot  be  directly  applied  to 
man.  A  1  /  dependence  of  the  post-iesonant  SAR  is  a!-o 
expected  from  antenna  theory.  Ihe  anticipated  I  f  post¬ 
resonant  hehaviot  t-  evident  m  I  ig.  2  bather  ituniettcai 
solutic'ns  suggested  a  signtl icaiitiv  faster  real  oif  between 
I  1  ’  and  f  4  [  I  ]. 

big  3  also  shows  the  whole-body  SAR  of  man  in  free 
space  with  I.  1.  and  an  incident  intensity  of  I  m\V  enr. 
Ihe  solid  curve  in  big  3  was  obtained  by  splitting  the  20 
values  obtained  in  numerical  solution-  fot  the  homoge¬ 
neous  realistic  model  of  man.  Points  on  big.  3  represent 
experimental  values  obtained  for  human-shaped  figurine- 
filled  with  saline  or  biological  phantom  mixtures  [IT’j,  (ISj. 
The  numerical  solutions  suggest  that  the  resonant 
frequence  is  77  MHz  which  is  somewhat  Ingliet  than  the 
value  of  bX  71  MHz  determined  experimentally . 

An  important  contribution  that  can  be  expected  from 
numerical  solutions  for  a  realistic  model  ot  man  is  the 
distribution  of  energy  deposition  through  the  model.  Pro- 
late-spheroidal  and  ellipsoidal  models  cannot  be  expected 
to  proxide  an  accurate  description  of  the  distribution  of 
energv  deposition  since  the  distribution  is  strongly  depen¬ 
dent  upon  detailed  geometry.  When  the  homogeneous 
approximation  of  complex  permittivity  is  used,  the  distti- 
bution  of  energy  deposition  m  the  realistic  model  of  man 
is  in  good  agreement  with  that  found  experimentally  for 
homogeneous  models  j!7|.  |1X]  I  able  I  contains  numeri¬ 
cal  and  experimental  values  for  the  distribution  of  energy 
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incident  mumsiis  =  I  mW,  mi*. 
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f  ig  *  ’  .'hofe-KuK  S  \R  for  model  of  m.in  m  free  space  Ilk  froni  to-h.uk. 

im.nl**  ini«?ns:«\  =  I  mW  cm*.  Squares  represent  experimental  points  for  sal  ne-filled  figurines. 
Friar.  ;ie  repr  *.*nt  experimental  values  for  fieunnes  ion!,nmiig  biological  pfi.mtoni  mixtures 
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!.-\l!  values  arc  in  units  watts  pc  kilogram  ;V>  kci  f<*<  m-.ixlcn! 
miensi!)  of  !  mV*  cri':;  F  I  .  front -to- hack 
•Phantom  models  of  sr.au  [<  k«ndh:  o«  ui .  I'^'M 
■'(Chen  f (  j!  .  Woi, 

VaLuisiu-nx  f ■  T  a  rejlisih  model  of  m.t..  h>  the  .1  i*t>  ■  fs 
C  .>mmetYi.i!  doll-  used  .0  form-  f  •'  'he  phantoms  .tie  one.!  '  have 
exaggerated  neck  constrict/  -it  Neck  hcaliii^;  for  the  do'R  ' -  .dvnit  ter’ 
times  [hal  of  the  expected  f.tr  nun 
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Fig.  4  Local  SAR  values  (waits  ke  per  mW'  cnt‘>  ■  !0o  for  homoge¬ 
neous  model  of  man  with  verlic  I  pola  'ration  at  80  MMr 

deposition  for  man  near  resonance  in  free  space.  The 
deposition  of  energy  in  the  regions  of  the  e\e,  heart,  thigh, 
and  calf  has  been  severely  underestimated  in  previous 
calculations  made  with  less  realistic  models  ( 1 1.  ( I  1 1. 

When  the  inhomogeneous  complex  permittivities  are 
used  with  the  model,  a  change  of  less  than  2  percent 
typically  occurs  in  the  whole-body  average  SAR.  but  the 
distribution  of  energv  deposition  is  changed  I  igs.  4  and 
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f  :g  4  I.sai  SAR  values  'w..--.  sc  ;  o  VV  on'1  .  j<gi  t..:  nhom... 

geneous  model  :  man  with  vertical  -  .or.-atior.  .a  g'*  Vffi/ 


5  illustrate  the  distribution  of  ibs  ■:  :>ed  energy  ir.  mar.  .it 
SO  MU/  >n  free  -pace  fo;  the  homogeneous  and  inhomo¬ 
geneous  models,  respectively.  One  difference  is  that  the 
inhomogeneous  model  has  reduced  absorption  in  regions 
with  high  bone  content,  burlier  workers  have  not  reported 
results  ;.>r  inhoni. 'geneous  models  having  a  celi-bv-ceil 
representation  of  the  dielectric  properties 

V  I  C  om  1  i  sp  ISIS 

Numerical  calculations  sit  absorbed  energy  deposition 
have  been  made  for  a  model  of  man  that  has  been  made 
with  careful  attention  to  both  biometric  and  anatomical 
diagrams.  \  total  of  ISO  cells  was  used  in  the  model, 
which  is  a  relatively  large  number  for  numerical  solutions 
but  facilitates  accurate  modeling. 

Calculated  post-resonant  absorption  and  distribution  of 
energy  deposition  through  the  body  have  been  found  to 
have  better  agreement  with  experimental  results  than  pre¬ 
vious  calculations  made  using  less  realistic  models. 
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Fine  structure  in  (he  whole-body  resonant  curve  for  radio-frequency  energy  deposition  in  man  can  be 
attributed  to  part-body  resonances.  As  for  head  resonance,  which  occurs  near  751)  MH/  in  man,  the  ah 
sorptive  cross  section  is  nearly  three  times  the  physical  cross  section  of  the  head.  The  arm  has  a  prominent 
resonance  at  150  MH/.  Numerical  solutions,  antenna  theory,  and  experimental  results  on  animals  have 
shown  that  whole-body  energy  deposition  may  be  increased  by  50  percent  or  more  because  of  multiple 
bodies  that  are  strategically  located  in  the  field.  Empirical  equations  lor  SARs  are  also  presented  along 
with  test  data  for  several  species  of  laboratory  animals.  Barbiturate  anesthesia  is  sufficiently  disruptive  of 
thermoregulation  that  STs  of  colonic  temperature  yield  energy  dose  values  in  several  mammals  that  com¬ 
pare  quite  favorably  with  those  based  on  whole-body  calorimetry. 


I.  INTRODUCTION 

Some  of  the  experimental  results  on  absorption  of  RF 
energy  that  have  been  obtained  from  animals  and  from 
models  of  man  have  recently  been  evaluated  in  the  light  of 
numerical  calculations  on  an  improved  model  of  man 
[Hagmann  et  at.,  1977).  Two  new  effects  that  were 
predicted  by  numerical  solutions  and  were  confirmed 
experimentally  are  part-body  resonances  and  multibody  ef¬ 
fects.  Both  are  described  in  this  paper  for  the  first  time. 
New  empirical  equations  for  the  SAR  for  EIL  orientation 
are  also  given  along  with  test  data  obtained  from  several 
vertebrate  species. 

2.  PART  BODY  RESONANCES 

Figures  1  and  2  show  part-  and  whole-body  SARs 
calculated  for  the  improved  model  of  man  at  various  fre¬ 
quencies  in,  respectively,  free  space  and  in  electrical  contact 
with  ground,  when  the  vector  of  the  electric  field  is  parallel 
with  the  model's  long  axis  (Fill.)  and  the  direction  of  pro¬ 
pagation  (k)  is  from  the  ventral  to  dorsal  aspect.  Fine  struc¬ 
ture  in  the  whole-body  SARs  is  due  to  part-body 
resonances. 

In  previously  reported  experiments  (with  figurines)  on 
whole-body  SARs  tor  man,  an  anomalous  ini  nase  in  the 
rate  ol  energy  deposition  was  observed  in  the  legion  of  470 
MH/  tor  Fill,  and  kill,  orientations  \Gundhi  et  ui.  1977). 
Theoretical  studies  of  a  10-cm  isolated  spherical  model  of 
the  human  head  have  shown  that  its  first  geometrical 
resonance  occurs  near  4.50  MH/  j  J  nines  and  Spiegel,  1974). 
Our  numerical  calculations  (Figures  I  and  2)  for  the  im¬ 
proved  model  of  man  (Figure  .7)  with  Fill,  have  shown  that 
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resonance  of  the  intact  human  head  occurs  at  a  frequency 
near  350  MHz.  For  head  resonance,  numerical  solutions 
give  the  S  parameter  (ratio  of  absorptive  cross  section  to 
physical  cross  section)  of  the  intact  human  head  as  3. 00  for 
man  in  free  space  and  2.66  for  man  on  a  ground  plane.  The 
magnitude  of  resonant  absorption  of  the  intact  human  head 
is  sufficient  to  cause  the  experimentally  observed  increase 
in  whole-body  SAR.  The  numerically  obtained  head- 
resonance  frequency  of  350  MHz  is  in  agreement  with  the 
experimental  data  presented  in  this  paper,  but  is  somewhat 
lower  than  the  450-  to  470-MHz  values  obtained  from 
earlier  experimental  data  and  is  also  lower  than  the  first 
geometrical  resonance  of  the  10-cm  isolated  sphere.  Further 
work  is  in  progress  to  resolve  these  differences. 

Joines  and  Spiegel  1 1 974)  calculate  that  the  ratio  of  the 
absorptive  cross  section  to  the  physical  cross  section  of  the 
isolated  human  head  is  approximately  1.1  at  geometrical 
resonance  near  450  MHz  and  1 .4  at  a  second  resonant  peak 
that  occurs  near  2. 1  GHz,  if  appropriate  allowance  is  made 
for  the  inhomogeneous  structure  of  a  multilayered  model. 
The  present  studies  indicate  that  the  magnitude  of  the 
geometrical  resonance  of  the  intact  human  head  is  more 
than  twice  that  of  the  isolated  human  head  for  either  of  the 
previously  identified  resonances. 

Table  1  gives  experimental  values  of  SARs  of  figurines 
filled  with  biologically  simulating  materials  with  kill,  pro¬ 
pagation  from  head  to  toe.  For  these  experiments,  ac¬ 
curately  scaled  figurines,  25.4,  33.0,  and  40.6  cm  in  length, 
were  exposed  in  free  space  to  radiation  at  2450  MHz.  The 
whole-body  dose  was  obtained  by  a  Thermonctics  Model 
2401 -A  gradient-layer  calorimeter.  The  heads  of  the 
figurines  were  attached  to  torsos  with  a  saline-soaked  cloth 
(to  maintain  conductive  contact),  which  permitted  easy 
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Fig.  1 .  Part-body  SARs  for  a  homogeneous  model  of  man  in  free  space;  Ell.,  k  ventral  to  dorsal,  incident 

power  density  of  I  mW  cm  '. 

separation  during  subsequent  measurements  of  head  SARs.  loss  by  evaporation  for  the  duration  of  its  placement  in  the 
To  determine  the  energy  dose  in  the  head  after  separation,  calorimeter. 

the  part  of  the  polyurethane  mold  that  held  the  biologically  We  believe  that  the  significance  of  the  head’s  resonance, 
simulated  head  was  wrapped  in  Saran  Wrap  to  prevent  heat  as  shown  in  Table  1,  should  be  extended  to  studies  of 
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Fig  2  Part-body  SARs  for  homogeneous  model  of  man  in  contact  with  a  ground  of  infinite  conductivity, 
Elf..  k  ventral  to  dorsal;  incident  power  density  of  I  mW  cm  *. 
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I  if.  5.  An  unproved  model  ol  man 


animals  F  or  animals  with  a  head  similar  in  si/e  to  that  of 
' he  laboratory  rat,  resonance  is  near  the  commercially 
significant  frequency  ol  2450  V1H/.  The  head  s  resonance 
was  not  obtained  in  earlier  numerical  solutions  where  the 
ratio  ol  head  10  whole-body  SAR  was  reported  as  1.07  al 
FIX)  MH/  and  1.09  at  400  MH/  \Chen  et  al..  1470). 

The  Irequency  ol  maximal  ciieigy  deposition  in  the  arms 
s  approximately  150  MH/  when  the  arm's  length  is  nearly 
tall  a  tree-space  wavelength.  Al  resonance,  the  arm  has  an 
absorptive  cross  section  that  is  2.5  limes  the  phvxical  cross 
section 

f  Mt  i  nnouv  i  iik  is 

li  ’las  been  shown  dial  die  energy  absorbed  O'  man  n 
lesonance  corresponds  to  ihe  energy  received  O'  a  halt 
wave  Jipnle  of  length  approximately  equal  lo  the  neighi  of 
die  'Man  \Oarulhi  el  ai.  1976.  I977|  Antenna  theory  has 
been  ised  io  calculate  the  effective  area  (or  gain)  pet  dipole 
loi  a  pair  ol  thin  dipoles  in  broadside  (sagittal)  conligura 
don  iKJK.  k  perpendicular  to  the  line  joining  t he  dipole's 
.enters),  f  igure  4  gives  Ihe  variation  of  effective  area  as  a 


function  of  spacing  lor  half-wave  dipoles  F-iguri  5  gives  du 
variation  of  effective  area  as  a  function  ol  Irequency  to- 
dipoles  of  length  !  75  meters  at  a  fixed  spacing  of  0.65  ' 
Standard  numerical  procedures  have  been  used  10  coinpuli 
the  required  values  of  mutual  impedance  of  dipoles  | Hake’ 
ei  al.,  I962| 

figure  4  shows  iilai  if  two  thin  hall-wa>  dipoles  an 
tangent,  eac)'  w  ill  receive  approximalely  one  hah  the  energy 
ii  would  receive  il  isolated  in  free  space  t  o  separations 
greater  than  two  01  ditee  Iree-space  wavelength.-,  the  amen 
nas  have  little  coupling  so  that  each  has  at  etlective  area 
approximating  the  isolated  tree-space  value  Ai  a  spacing 
ol  0.65  a.  each  dipole  will  receive  about  50  percent  mon- 
energy  that  it  would  it  it  were  isolated,  Figure  5  shows  t hat 
al  a  fixed  spacing  ol  0.65  iwe'  dipoles  with  lengths  of  I  .'5 
meters  will  have  an  enhancement  of  aboir  50  peteent  in 
received  energy  over  a  fairly  wide  band  ol  frequencies 
ralher  lhan  only  ai  ihe  frequency  for  whict  ihey  are  half¬ 
wave  dipoles  a  series  of  oscillations  and  then  a  roll-off  are 
seen  ai  ihe  high-trequency  end  in  Figure  5  Since  antenna 
theory  has  been  shown  10  predict  accurate  SARs  for  human 
models  and  lor  animals  with  and  without  nearby  reflectors 
[Gantlhi  ei  at. .  1977).  q  is  anticipated  that  the  variation  of 
SARs  tor  the  Kill,  orientation  for  near-resonant  larget-  as  a 
function  ol  spacing  will  be  given  by  t hi  curve  of  Figure  4 
The  variation  ol  SAR  with  frequency  lo>  biological  targets 
that  arc-  separated  a  distance  ol  0.65  h  would  similaro  be 
given  by  the  curve  in  Figure  5. 

Numerical  solutions  have  been  obtained  lo-  du 
multibody  effect  using  the  improved  model  ol  mar 
| Hanmann  ei  al.,  1 977 1  The  solutions  have  been  obtain-.-' 
for  ihe  configuration  ol  two  men  standing  side  by  sidi  ha 
ing  Kill,  and  k  directed  v en trails .  I.acF  ol  symmelrv  tv.- 
ween  the  left  and  right  sides  in  each  model  has  necessitate.' 
the  use  ol  18(1  independent  cells  in  llu  lormulation  o'  -in 
matrix.  Since  there  are  three  unknown  field  component'  -ot 
each  cell,  the  matrix  i'  540  hy  540  complex,  which  .  o- 
responds  lo  a  system  ol  1080  simultaneous  equation  n 
1080  real  unknowns.  Tola'  compute-  tunc  'or  each  nui: 
body  problem  was  19  lion r s.  25  minutes  .viili  deduaic.l 
usage  ol  a  PDI’-IO  digital  compute:  Weak  diagonal 
dominance,  which  is  inherent  in  the  matrix,  .ontrihuk-  to 
the  conditioning  to  allow  reasonable  round-off  ertot'  with 
such  a  latgc  matrix.  Noniterative  main'  solutions  wen  re 
quired  lor  convergence. 

The  data  points  in  Figute  4  and  5  represent  mtme-ical 
solutions  of  calculated  enhancement  11  SAR  lot  human 
models  due  to  two-body  clu-.is  Tin  numerical  solutions 
arc  reasonably  consistent  wm  aiiicnn,-  dieorx  lot  mat.  at 
resonattvi  as  is  show  11  m  Figuu  4.  bm  signifisani  depat  tine 
Iti'tu  anieuna  theory  0  sect  al  subiesonaiH  liequeiisies 
(l  igure  5i  -Ml  calculations  bused  on  ameiiu.l  theory  have 
been  /err>  older;  that  is.  the'  have  no  collection  lot  non 
sinusoidal  current  disit  thultotis.  Siciutkam  lust  and 
second  oivlet  correction'  arv  tequtted  when  die  antennae 
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TABU:  I.  lix  peri  mental  SAR  lor  kill,  propagation  from  head  to  toe.  Phantom  models  of  man.  Incident 

power  density  =  10  mW  cm 
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arc  not  very  thin  and  the  halt  'englh  is  much  different  from  MHz.  The  whole-body  SAR  is  45.4  percent  greater  than 

X/4  [Am*,  1956],  Thus,  the  i/able  lateral  extent  of  the  that  for  man  in  free  space.  In  general  the  body  halves 

model  of  man  is  thought  :o  he  responsible  for  the  dif-  closest  to  each  other  have  somewhat  less  enhancement  than 

ference  between  numerical  solutions  and  antenna  theory  at  the  more  distant  body  halves,  but  the  asymmetry  is  reversed 

low  frequencies.  in  the  arms. 

Some  asymmetry  is  present  in  the  distribution  of  A  greater  enhancement  in  the  SAR  by  a  factor  as  large  as 
absorbed-energy  density  between  the  two  body  halves  in  all  2.5  is  anticipated  for  the  inner  targets  trom  antenna  theory 

numerical  solutions  obtained  for  the  two-body  effect.  for  three  or  more  bodies,  each  of  which  is  separated  0.65  X 

Figure  6  illustrates  the  calculated  local  SAR  values  for  one  from  the  closest  neighbors.  Table  2  gives  the  results  ot 

of  two  models  of  man  when  the  separation  between  the  preliminary  experimental  tests  of  the  multibody  effect  with 

models  is  0.65  X  =  2.53  m  at  the  resonant  frequency  of  77  anesthetized  Long-Evans  rats.  The  configurations  used  in 
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Fig.  4.  Variation  of  gain  enhancement  as  a  function  of  spacing  of  a  Fig.  5.  Variation  of  gain  enhancement  with  frequency  for  a  pair  of 

pair  of  half-wave  dipoles  in  broadside  (sagittal)  configuration.  1.75-me»er  dipoles  in  broadside  (sagittal)  configuration  at  a  spac- 

Data  points  represer.  calculated  enhancements  of  SARs  in  human  iug  of  0.65  X.  Data  points  represent  calculated  enhancements  of 
models  due  to  two-body  effect.  SARs  in  human  models  due  to  two-bodv  effect. 
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Fig.  6.  Local  SAR  values  (mW  g  ’  per  mW  cm  ')  x  100  for  one  of 
l«o  homogeneous  models  of  man  having  a  cenler-to-cenler  spac¬ 
ing  of  0  65  X  at  77  MH/.  Incident  power  density  of  1  mW  cm 
distance  between  closest  arms  of  the  two  models  =  0.51  X. 


Fig.  7.  Geometry  used  in  measurement  of  enhancement  o!  SARs 
due  to  proximity  of  other  biological  bodies,  (a)  two-a-iimal  ex¬ 
posure;  (b)  three-animal  exposure. 


preciably  greater  enhancement  in  SAR  than  is  obtained  in 
either  animal  of  a  two-animal  exposure.  Experiments  ai 
several  frequencies  with  varying  interanimal  spacing  need 
to  be  done  before  the  full  implications  of  multibody  effects 
on  the  SAR  are  understood. 


the  tests  are  shown  in  Figure  7.  Note  that  for  each  of  the 
animals,  an  experimental  frequency  considerably  higher 
than  its  resonant  frequency  was  used,  confirming  the  fairly 
wide-band  nature  of  the  multibody  effect.  Note  also  that 
the  central  animal  in  a  three-animal  ensemble  has  ap- 


4.  EMPIRIC  AI.  EQUATIONS  FOR  W HOI  F-BOIA 
AVERAGE  SAR  UNDER  FREE-SPACE  IRRADIATION 

The  su>  -<onant  frequency  dependence  [ Gandhi  er  at.. 
1977)  an.  observed  1//  dependence  in  the  suprareso- 
nant  region  have  been  used  to  develop  the  empirical  equa- 


I  ADI  F  2.  ExperimcnlaOrnhancemcnt  in  SAR  due  to  proximity  toother  biological  bodies.  Anesthetized 
l  ong-F.vans  rats  with  Fit.  at  245(1  MM/.  Centcr-to-center  interanimal  separation  =  0.65  X  (  -  7  96  cm). 
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TABLE  3.  Empirical  equations  for  SARs  in  human  models  and  in  laboratory  animals.  Free*  space  irradia 
lion,  KIK  polarization,  resonant  frequency,/.  =  11.4  l .  1  GHz. 
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lions  in  Tabic  3  for  the  whole-body  average  SAR  for 
human  models  and  for  laboratory  animals  with  E  polariza¬ 
tion.  While  the  numbers  of  biologically  simulating  models 
of  man  are  the  same  as  those  reported  earlier  [Gandhi  el  a!., 
1977),  the  coefficients  for  S....  and  SARs  of  laboratory 
animals  are  59  percent  higher  and  are  derivable  from  anten¬ 
na  theory  for  the  capture  cross  section  of  0.38  \  dipoles  by  a 
procedure  similar  to  that  of  Gandhi  el  at.  [I976|.  The  re¬ 
vised  higher  values  for  SARs  are  in  good  agreement  with 
experimentally  determined  values  (Table  4)  for  several 
animal  species.  For  the  data  in  Table  4,  the  animals  were 
anesthetized  (sodium  pentobarbital,  45  mg  kg  ')  to  prevent 
movement  and  to  maintain  thereby  a  fixed  orientation  and 
also  to  limit  pharmacologically  the  normal  thermo¬ 
regulatory  functions  of  the  animal  ( Putthoff  el  at.,  1977). 
The  experimental  SAR  was  calculated  from  the  increase  in 
colonic  temperature  (cloacal  temperature  in  the  case  of 
lizards  and  birds)  after  three  minutes  of  tree-space  irradia¬ 
tion  at  100  mW  cm'J.  For  the  500-g  rat  and  the  120. 3-g 
dove,  freshly  euthanized  animals  were  used  and  the 
calorimeter  was  employed  to  determine  the  whole-body 
dose;  the  procedure  was  similar  to  that  outlined  by  Phillips 
el  ai.  (1975).  The  points  to  note  in  Table  4  are:  (1)  The 
SARs  from  the  empirical  equations  for  laboratory  animals 
(for  the  supraresonant  region  shown  in  Table  3)  correlate 
well  with  experimental  measurements  on  several  species 
that  range  in  mass  from  18.8  to  2245  g,  and  in  length,  from 
8  to  44  cm.  For  these  animals,  a  whole-body  SAR  varying 
by  a  factor  of  36.1:1  is  observed  at  50  mW  cm'*  for 
2450-MFIz  irradiation.  (2)  For  ring  doves  there  is  a  signifi¬ 
cant  difference  between  experimental  values  of  SARs  and 
those  calculated  using  the  empirical  equation.  The  reason 


for  this  discrepancy  is  not  clear.  Subtracting  the  mass  of  the 
plumage  from  the  whole-body  mass  of  the  birds  increases 
the  values  calculated  via  the  empirical  equation  by  a  few 
percent  but  does  not  explain  the  observed  difference  for 
this  avian.  (3)  Reliable  correlations  between  whole-body 
calorimetric  values  on  freshly  killed  animals  and  SARs 
based  on  measurements  of  colonic  temperature  of 
anesthetized  animals  confirm  the  rapid  hemodynamic 
dispersion  of  heat  in  living  animals.  Previous  measurements 
of  field  intensities  needed  for  comparable  times  to  convul¬ 
sion  [Gandhi  and  Hunt ,  in  preparation]  at  different  fre¬ 
quencies  has  also  demonstrated  that  whole-body  integral 
dose  and  dose  rates  are  important  parameters  in  the  study 
of  living  animals. 

5  CONCLUSIONS 

Numerical  solutions  based  on  an  improved  model  of  man 
have  demonstrated  that  the  deposition  of  energy  at 
supraresonant  frequencies  has  a  fine  structure  that  is 
ascribablc  to  resonance  of  body  parts.  Failure  to  allow  for 
the  detailed  structure  of  the  body  of  man  through  use  of 
prolate  spheroidal,  ellipsoidal,  and  other  idealized  model' 
would  cause  the  part-body  resonances  to  he  overlooked. 
The  enhancement  of  electromagnetic  energy  deposition  in 
the  head  at  head-resonant  frequencies  may  be  important  in 
future  studies  of  the  blood-brain  barrier  and  of  animal 
behavior.  Numerical  solutions,  antenna  theory,  and  ex¬ 
perimental  results  with  animals  have  shown  that  energy 
deposition  is  altered  by  multibody  effects.  A  maximum  of 
about  50  percent  enhancement  in  SAR  is  found  with  two 
bodies,  but  significantly  greater  enhancement  is  possible 
with  three  or  more  bodies,  interbody  spacing  being  critical. 
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I  Mil  l  4.  Whole  body  S  \H>  -n  iii\V  g  '  for  Kill.  ..>h  ot 

'Cveru*  sptvics  o!  Iabt>iau>ry  ani’naK  irradiated  by  24*9  MU/ 
energy  n.  tree  space.  I  reshly  mihu.ii/ed  animals  and  a  liicrmo- 
neties  Model  2401 -A  gradieru-la>er  calorimeter  were  u*ed  to 
measure  ihe  500  g  rat  and  the  120.3-g  dose.  Other  measurements 
are  based  on  anesthetized  animals. 
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Antenna  theory  and  the  observed  !//  dependence  at 
supraresonant  frequencies  have  beer  u  ed  to  develop  em¬ 
pirical  equations  for  the  whole-bod,  liAR  for  models  of 


i.ian  and  for  laboratory  animals  with  the  vector  of  the  t 
field  parallel  to  ihe  body’s  long  axis.  The  empirical  SAR 
values  are  shown  to  be  in  good  agreement  with  the  ex¬ 
perimental  values  obtained  for  several  animal  species. 
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.1.  led  So  ec  polarization  currents  in  the  images  contribute 
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where  A  is  an  m  by  m  matrix  formally  defined  by  A'  = 
first  m  rows  of: 
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In  implementing  compression,  the  computation  of 
matrix  .-1  is  made  by  making  n  modified  passes  through 
the  loops  that  would  normally  be  used  to  accumulate  the 
scattering  matrix. 

We  have  used  compression  with  image  theory  to  treat  a 
body  with  one  plane  of  symmetry  above  a  ground  plane,  a- 
is  shown  in  f  igure  1 .  The  field  relationships  used  for  a  body 
with  one  plane  of  symmetry  in  front  of  a  flat  leflector  are 
shown  in  figure  2.  It  is  also  possible  to  use  compression 
with  image  theory  to  treat  a  body  in  a  corner  reflector,  it 
the  corner  reflector  has  a  dihedral  angle  equal  io  1H0  ,\ 
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where  V  is  an  integer,  field  relationships  used  for  a  body 
with  one  plane  of  symmetry  in  front  of  a  90  cornet  reflec¬ 
tor  are  shown  in  Figure  3. 
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I  here  are  two  types  of  symmetry  relationships  used  for  h 
values  m  the  ihree  figures:  (1)  Two  points  loeated  sym- 
metrieally  in  respect  to  a  plane  of  symmetry  of  the  body 
that  contains  both  k  and  K.  will  necessarily  have  K  values 
that  arc  mirror  images  so  that  the  two  halves  are  in¬ 
distinguishable  under  reflection.  (2)  Two  points  located 
sy  tnmelrieally  m  respect  to  a  plane  of  symmetry  of  the  body 
that  contains  k  and  is  perpendicular  to  E  will  necessarily 
have  1.  values  that  are  negated  mirror  images,  so  that  if  the 
values  are  exchanged,  they  will  correspond  to  the  solution 
tor  oppositely  directed  K  .  This  relationship  is  also  found 
tor  two  points  located  symmetrically  in  respect  to  a  ground 
plane  or  to  a  flat  reflector  as  required  for  symmetrical  mo¬ 
tion  ol  a  cltatgc  and  its  image. 

An  alternative  to  the  use  of  image  theory  with  compres¬ 
sion  is  to  use  special  Green's  functions  in  place  of  the  l'ree- 
space  Green's  lunction.  Special  Green's  functions  are 
available  that  could  be  used  to  account  for  scattering  ob¬ 
jects  located  near  targets  having  various  shapes  and  even 
finite  conductivities  | Tm.  I •  >7 1 ;  Banos,  1966).  We  have  not 
yet  pursued  the  use  of  such  Green’s  functions. 

fhe  results  of  calculations  for  grounding  and  reflector 
effects  are  described  in  the  next  sections. 

3.  GROUND  EFFECTS 

All  calculations  of  ground. ng  effects  have  assumed  that  a 
man  is  standing  on  or  above  a  perfectly  conducting  ground 
plane  that  is  infinite  in  extent.  The  incident  field  is  verti¬ 
cally  polari/ed  with  k  directed  from  the  vernal  to  the  dorsal 
aspect  of  the  man.  Internal  fields  in  the  model  are  (he  same 
as  those  in  one-half  of  a  double-man,  which  consists  of  the 
model  and  its  image  in  tree  space.  A  180-cell  model  of  the 
5()th-pcrccntile  standard  man  has  been  used  for  the  calcula¬ 
tions  \Huxmann  el  a!.,  1977|. 

I  he  calculated  resonant  frequency  ol  man  standing  on  a 
ground  plane  is  4'  Mil/  as  compared  to  77  MU/  for  man  in 
free  space,  f  he  specific  absorption  rate  (SAK)  of  man  on 
the  ground  plane  at  the  reduced  resonant  frequency  is  32.5 
percent  greater  than  the  SAK  for  man  in  free  space  at  the 
free-spaee  resonant  frequency.  At  lower  frequencies,  the 
SAK  of  man  standing  on  a  ground  plane  may  be  increased 
considerably  from  that  ol  the  same  model  in  free  space. 
I  his  increase  results  from  these  frequencies  being  closer  to 
resonance  for  the  grounded  man.  For  example,  at  It)  Mil/ 
the  SAK  of  mail  standing  on  a  ground  plane  is  0.0163  \V 
kg  1  per  ni\V  cm  ,  which  is  about  seven  times  that  found 
tor  the  same  model  in  Iree  space. 

Figures  4  and  5  illustrate,  respectively,  the  distribution  ol 
ahxot bed  energy  in  man  at  It)  Mil/  in  free  space  and  in  con¬ 
tact  with  the  ground  plane.  I  lie  enhancement  of  local 
etietgv  deposition  due  to  the  grounding  effect  is  approx 
imatclv  a  factor  ol  60  in  the  area  ot  the  heel 

I  ignre  6  shows  the  pat t  body  and  whole  body  SAK  lor 
the  model  ol  man  with  dillerent  spacings  from  the  ground 


tig.  4.  local  SAK  values  (\V  kg  1  per  mW  cm  ')  x  100  for 
inhomogeneous  model  ol  man  in  free  space  at  )0  MH/ 


plane  at  10  MH/.  Figure  6  is  consistent  w iih  the  experimen¬ 
tal  observation  on  models  that  a  small  separation  that 
breaks  electrical  contact  with  the  ground  plane  is  sufficient 
to  eliminate  much  of  the  grounding  effect  | Gandhi  ei  at., 
1975).  Several  calculations  made  for  the  grounded  resonant 
frequency  of  47  MH/  show  a  tall-off  in  magnitude  of 
grounding  effects  with  increasing  distance  from  the  ground 
plane,  which  is  similar  to  the  results  at  It)  MH/  (Figure  6). 
A  model  composed  of  292  cells  was  made  by  dividing  each 
of  the  six  lowest  cells  in  each  leg  of  the  180-cell  model  into 
eight  identical  cubical  cells.  Six  solutions  made  for  the 
292-cell  model  and  duplicated  with  the  180-cell  model  sug¬ 
gest  that  the  fall-off  of  the  grounding  effect  with  increasing 
distance  from  the  ground  plane  may  be  slightly  slower  than 
is  found  with  the  180-cell  model. 

Figure  7  shows  part-  and  whole-body  SARs  for  the  model 
of  man  standing  on  a  ground  plane  at  various  frequencies. 
The  part-body  resonances  shown  in  Figure  7  have  also  been 
observed  in  calculations  for  the  same  model  of  man  in  free 
space  [Gandhi  el  ul.,  1977a].  The  enhanccmci.  in  absorp¬ 
tion  of  energy  due  to  the  presence  of  a  ground  plane  is 
greatest  in  the  leg  and  least  in  the  head.  Figure  8  shows 
whole-body  SARs  for  the  model  of  man  both  in  free  space 
and  when  standing  on  a  ground  plane  at  various  frequen¬ 
cies.  Grounding  effects  are  most  prominent  at  low  frequen- 
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I  ig.  5.  l  ocal  SAR  values  (\\  kg  Per  mV.  cm  |  •  UK)  lor 
inhomogencoii-  model  ot  man  standing  on  a  ground  plane  a! 
Ill  Mil/ 

cies.  Above  a  frequency  of  aboul  200  MH/,  pari-  and 
whole-body  SAR  values  show  little  dependence  on  the 
presence  of  a  ground  plane. 
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All  calculations  of  reflector  ettccls  have  assumed  that  a 
mail  is  standing  in  from  of  a  reflector  that  is  perfectly  con 
ducting  and  infinite  in  estenl.  The  incident  field  is  wrlicalb 
polarized  with  k  directed  from  the  ventral  lo  the  dorsal 
aspect  of  the  man.  One  image  is  required  for  a  flat  reflector 
and  three  are  needed  for  a  90  corner  reflector. 

The  enhancement  of  SAR  due  to  reflector  effects  is 
approximately  the  same  as  the  ratio  of  effective  area  of  a 
half-wave  dipole  with  a  reflector  to  that  of  an  isolated 
dipole  in  free  space  [Gandhi  el  al.,  1977b],  Antenna  theory 
has  been  used  to  perform  numerical  calculations  made  for 
the  model  of  man.  Standard  procedures  have  been  used  for 
computation  of  the  required  values  of  mutual  impedance  of 
dipoles  | Baker  and  LaGrone,  1962], 

f  igures  9-12  show  the  enhancement  of  energy  absorption 
as  calculated  for  the  model  ol  man  with  reflector  effect' 
and  the  corresponding  gain-enhancement  factors  calculated 
for  thin  dipoles  from  antenna  theory.  V  alues  of  enhance 
mem  for  lire  whole-body  average  SAR  have  been  used  in 
preparing  figures  9-12.  but  the  curves  for  part -body  energy 
absorption  are  nearly  identical,  which  indicates  that  the 
distribution  of  energy  through  the  model  is  almost  indcpcii 
dent  of  the  presence  of  a  reflector.  I  spenmental  values  in 
Figures  9-12  are  for  models  of  man  and  were  reported 
earlier  [Gandhi  el  al.,  1977b]. 

f  igures  9  and  10  show  the  variation  of  enhancement  ot 
energy  absorption  as  a  function  of  spacing  from  a  reflect vm 
at  the  free- space  resonant  frequency  of'  MHz.  In  figure  '> 
the  numerical  values  for  the  model  of  matt  ate  marked!' 
below  the  curve  for  thin  dipoles  at  small  values  of  'cptra 
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I  If.  ft.  W  hole-body  SARs  in  man  lot  different  -pavings  Horn  the  ground  at  in  Mil.'.  I-  III  .  k  ventral 
dorsal,  incident  power  density  ol  I  mVV  cm 
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Fig.  7.  Part-body  SARs  for  a  homogeneous  model  of  man  on  a  ground  plane;  EIL,  It  ventral  to  dorsal, 
incident  power  density  of  1  mW  cm'1. 


tion.  Antenna  theory  requires  that  the  gain  be  small  when 
the  spacing  is  not  much  greater  than  the  thickness  of  the 
dipole,  but  such  variation  is  missed  in  the  thin  dipole 
approximation. 

Figures  1 1  and  12  show  the  variation  of  enhancement  of 
energy  absorption  as  a  function  of  frequency  for  fixed 
ratios  of  spacing  to  wavelength.  The  antenna-theory 
calculations  were  made  for  a  thin  dipole  with  a  length  of 
1.7J  meters,  which  is  identical  to  the  height  of  the  model. 
Figures  1 1  and  12  show  that  the  reflector  effects  are  found 
over  a  wide  range  of  frequencies.  Values  for  the  model  of 
man  have  a  pronounced  roll-off  at  higher  frequencies, 
which  is  not  seen  in  calculations  for  the  thin  dipoles.  In 
Figures  1 1  and  12  it  can  be  seen  that  the  rate  of  roll-off  of 
the  SAR  at  high  frequencies  is  consistent  with  experimental 
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Fig.  8.  Whole-body  SARs  for  homogeneous  model  of  man;  EIL,  k 
ventral  to  dorsal,  incident  power  density  of  1  mW  cm''. 


results  that  are  obtained  from  models  of  man.  Antenna 
theory  requires  that  the  gain  decrease  at  higher  frequencies 
when  the  effective  area  of  the  dipole  images  is  blocked  by 
the  physical  cross  section  of  the  dipole,  but  such  variation  is 
missed  in  the  thin-dipole  approximation. 

All  gain  calculations  have  been  based  on  zero-order 
antenna  theory;  that  is,  they  have  no  correction  for  non- 
sinusoidal  current  distributions.  Significant  first-  and 
second-order  corrections  are  required  when  the  antennae 
are  not  very  thin  and  the  half-length  is  much  different  from 
X/4  [King.  1956],  The  sizable  lateral  extent  of  the  model  of 
man  is  thought  to  be  responsible  for  the  difference  between 
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Fig.  9.  Variation  of  gain  enhancement  wuh  spacing  for  a  half-wave 
dipole  in  from  of  a  flat  reflector.  Squares  represent  numerical 
solutions  for  enhancement  of  SARs  in  man  at  77  MHz  (L,\  = 
0.45).  Solid  circles  represent  experimental  values  for  anatomically 
scaled  Figurines  at  L<\  -  0.417. 


-  45  - 


28 


HAGMANN  AND  GANDHI 


4 ’ 


aiswHCE  in  awLosim 


Fig.  10.  Variation  of  gain  enhancement  with  spacing  for  a  half¬ 
wave  dipole  in  front  of  a  90°  corner  reflector.  Squares  represent 
numerical  solutions  for  enhancement  of  the  SAR  in  man  at  77 
MHz  (t/X  =  0.45).  Solid  circles  represent  experimental  values  for 
anatomically  scaled  Figurines  at  L/\  =  0.417. 

numerical  solutions  and  antenna  theory  at  low  frequencies 
in  Figure  12. 

The  numerical  solutions  indicate  that  for  frequencies 
near  resonance,  the  enhancement  in  SAR  due  to  a  reflector 
is  approximately  equal  to  the  enhancement  in  gain  of  a 
half-wave  dipole  with  the  same  reflector  configuration.  Ex¬ 
perimental  studies  on  models  of  man  have  also  shown  that 
the  enhancement  of  SAR  is  equal  to  the  enhancement  in 
gain  of  a  dipole,  but  the  gain  must  be  calculated  for  the 
finite  experimental  reflector  rather  than  for  an  infinite 
reflector  having  the  same  dihedral  angle.  For  certain  com¬ 
binations  of  length  and  width  of  the  reflecting  sheets,  a 
finite  comer  reflector  may  exhibit  supergain  -  a  gain 
significantly  greater  than  that  for  infinite  reflecting  planes 
[Cottony  and  Wilson,  1958;  Schelkunoff  and  Friis,  1952J. 


Fig.  II.  Variation  of  gain  enhancement  with  frequency  for  a 
1.75-meter  dipole  0.75  X  from  a  flat  reflector.  Squares  represent 
numerical  solutions  for  enhancement  of  SARs  in  man.  Solid  circles 
represent  experimental  values  for  anatomically  scaled  Figurines. 
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Fig.  12.  Variation  of  gain  enhancement  with  frequency  for  a 
1.75-meter  dipole  1.5  X  from  a  90°  corner  reflector.  Squares  repre¬ 
sent  numerical  solutions  for  enhancement  of  SARs  in  man.  Solid 
circles  represent  experimental  values  for  anatomically  scaled 
Figurines. 

Examples  of  the  effect  of  such  supergain  combinations  are 
seen  in  the  experimental  values  in  Figures  10-12. 

Figures  13  and  14  were  copied  from  an  experimental 
study  of  finite-size  comer  reflectors  [ Cottony  and  Wilson, 
1958).  In  both  figures  the  angle  of  the  aperture  was  ad¬ 
justed  for  maximum  gain  at  a  particular  distance  of  the 
dipole  from  the  reflector.  For  reflector  lengths  greater  than 
about  0.5  X  (Figure  13),  the  aperture  angle  is  90°  and 
distance  of  the  dipole  from  the  reflector  is  1.5  X.  For  reflec¬ 
tor  lengths  greater  than  1 .0  X  (Figure  14),  the  aperture  angle 
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Fig.  13.  Contours  of  constant  gam  for  a  corner-reflector  antenna 
with  the  dipole  in  third  position  for  various  reflector  sizes.  For 
lengths  greater  than  -  0.5  X,  the  aperture  angle  is  90’  and  the 
distance  of  the  dipole  from  the  reflector  is  1.5  X  [Cotranv  and 
Wilson,  19581. 
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I  ;g  14  (  .mi. 'ins  ,.l  soust.iiu  gum  ini  ,i  s  or  ncr  -  reflector  amentia 
a  rtf.  the  dipoli  in  -es.oi;d  position  tor  vunou-  relleetor  sizes.  1  oi 
•  eng  ( hx  gustier  than  1.0  A.  ihe  aperture  angle  is  65  alls!  the 
distance  of  the  dipole  from  the  reflector  is  1.2  X  | Cottanv  and 
Hi/utn.  W5H|. 

is  65  and  di'Utnce  of  the'  dipole  from  the  relleetor  is  1 .2  a. 

Note  in  I  teiiie  13  that  lor  proper  combinations  ol  length 
and  widths  of  the  refleetmg  planes,  the  eiihaiieement  in 
S\K  o\et  the  free  space  value  may  he  3(1  to  4?  percent 
greater  Ilian  that  lot  inlmire  reflectors  I  or  such  snpergain 
combinations  the  SAK  mac  he  enhanced  (see  figure  I  3)  be 
as  much  as  14  to  14.'  dB  (a  lactor  ol  25  to  28)  rather  than 
12.9  dB  (19.5)  anticipated  for  infinite  reflectors,  or  the  fac¬ 
tor  of  16  that  is  projected  from  consideration  of  the  local  l. 
fields  in  the  standing-wave  pattern  of  a  VO  corner  reflec¬ 
tor.  As  is  shown  in  f  igure  13,  the  dimensions  2  X  x  3.5  X  of 
the  reflectors  used  for  the  experimental  results  of  f  igures  10 
and  12  correspond  to  one  of  the  supergain  combinations. 

Supetgain  combinations  of  length  and  width  are  also  pre¬ 
sent  lor  a  65  corner  reflector  (see  figure  14).  Antenna 
calculations  with  image  theory  are  restricted  to  corner 
reflectors  having  dihedral  angles  of  "  'here  K  is  an 

integer.  As  seen  in  figure  14,  gain  c.  'll  is  not 

limited  to  such  special  angles,  and  we  expect  tna,  .  .nil'ieant 
enhancements  of  the  SAR  will  occur  for  a  wide  variety  ol 
reflectors, 

5.  (  ONI  I  t  SKINS 

Image  theory  has  been  used  to  treat  grounding  and 
relleetor  effects  so  that  all  calculations  are  limited  to  a 
perfectly  conducting  ground  plane,  or  to  reflectors  with  an 
infinite  extent.  F  xperimcntal  results  suggest  that  reflectors 
ol  finite  si/c  mav  provide  significantly  more  or  less 


enhancement  in  energy  deposition  than  that  calculated  tor 
infinite  reflectors.  Ihe  effect  ol  finite  conductivity  on 
grounding  or  reflector  effects  with  models  of  man  has  not 
yet  been  evaluated  either  theoretically  or  experimentally. 

Ihe  decrease  in  resonant  I  requeue)  and  alteration  in 
energy  deposition  due  to  grounding  may  be  explained  by 
considering  a  double-length  man  in  tree  space.  If  the  model 
ot  man  had  a  cros-  section  independent  of  height  and  a 
complex  permittivity  independent  of  frequency,  the 
grounded  resonant  I  requeue)  would  he  one-half  the  tree- 
space  resonant  frequency,  hut  a  ratio  of  0.610  is  obtained 
for  ihe  model  that  has  been  used. 

The  dependence  of  reflector  el  feels  on  spacing  and  I  re¬ 
queue)  are  in  agreement  with  the  gain  enhancement 
calculated  for  dipoles  hv  antenna  ilieoiy  The  enhancement 
m  energy  deposition  due  to  icflcctor  efteets  becomes  small 
lor  small  values  ol  separation  or  for  high  frequencies, 
which  effects  agree  with  antenna  theory  only  if  finite  width 
of  the  dipole  is  considered. 
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ABSTRACT 

A  realistic  human  model  and  improved  numerical  methods  have 
been  used  for  calculation  of  deposition  of  electromagnetic  energy. 
Unlike  earlier  solutions  both  the  average  absorption  and  the  dis¬ 
tribution  of  absorbed  energy  within  the  model  are  in  good  agree¬ 
ment  with  experimental  measurements  made  using  phantom  models.  The 
distribution  of  absorbed  energy  is  frequency-dependent  and  may  be 
explained  in  terms  of  resonance  of  the  various  body  parts.  Numer¬ 
ical  solutions  for  man  near  a  ground  plane  and  near  reflectors  are 
presented  for  the  first  time.  At  10  MHz  the  specific  absorption 
rate  (SAR)  of  man  standing  on  a  ground  plane  is  about  seven  times 
that  for  man  in  free  space.  Multibody  effects  have  been  predicted 
from  antenna  theory  and  observed  with  experiments  using  anesthetiz¬ 
ed  rats.  For  two  resonant  targets  separated  by  0.65  X  an  increased 
SAR,  170  percent  of  the  free-space  value,  has  been  observed. 


Proceedings  of  the  Workshop  on  Physical  Basis  of  Electromange  tic 
Interactions  uith  Biological  Systems ,  University  of  Maryland, 
June  15-17 ,  1977  (Editors:  L.  S.  Taylor  and  A.  Y.  Cheung) 
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DESCRIPTION  OF  MODEL  AND  NUMERICAL  TECHNIQUES 

Solutions  in  electromagnetics  are  facilitated  by  choosing  a  model  of 
simple  geometry.  Then  an  exact  analytical  solution  may  be  possible,  or  a 
numerical  solution  may  be  expressed  in  terms  of  easily  constructed  whole- 
domain  basis  functions.  We  have  rejected  such  methods  since  emphasis  was 
placed  on  a  realistic  model  of  man.  When  simple  geometries  are  not  possible, 
one  can  always  attempt  a  solution  by  "brute-force"  using  moment  methods  with 
a  subsectional  basis.  We  have  used  moment  method  solution  of  the  electric 
field  integral  equation  with  a  pulse  function  basis  and  delta  functions  for 
testing1 . 

A  total  of  180  cubical  cells  of  various  sizes  were  used  to  obtain  a  best 
fit  of  the  contour  on  diagrams  of  the  50th  percentile  standard  man2.  Sizes 
and  placement  of  the  cells  is  shown  in  Fig.  1.  Anatomical  cross  sectic  3 3,4 
were  used  in  partitioning  bone,  fat,  skin,  muscle,  lung  tissue,  air,  heart, 
brain,  kidney,  liver,  and  spleen  through  the  cells.  Whole-body  volume  frac¬ 
tions  of  each  tissue  type  are  in  agreement  with  published  values.  Properti.t, 
reported  in  the  literature  for  the  tissue  types5-*7  have  been  used  to  calcu¬ 
late  the  volume-weighted  complex  permittivity  of  each  cell. 

Previously  reported  solutions  for  models  of  man  have  a  ratio  of  239.1 
for  energy  deposition  in  one  pair  of  adjacent  cells  at  10  MHz8.  The  arrange¬ 
ment  and  different  sizes  of  cells  in  the  present  model  cause  the  maximum  ratio 
of  energy  deposition  for  a  pair  of  adjacent  cells  to  be  8:1  at  the  same  fre¬ 
quency.  The  pulse  function  approximation  is  invalid  if  conditions  force  sig- 
nificant  variation  of  E  within  the  cells  and  a  large  difference  between  cal¬ 
culated  energy  deposition  in  adjacent  cells  suggests  the  existence  of  such 
variations. 

If  A  is  the  side  of  a  cubical  cell  and  k  represents  the  magnitude  of  the 
complex  propagation  vector  within  the  cell,  then  substantial  variation  of  the 
electric  field  must  occur  within  the  cell  if  kA  >  '/(T 9 .  Thus  the  pulse  func¬ 
tion  approximation  may  not  be  justified  at  frequencies  above  200  MHz  with  the 
new  model  suggesting  an  increasing  error  at  higher  frequencies.  Volume 
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weighting  of  the  complex  permittivity  within  each  cell  may  also  be  unjusti¬ 
fied  at  frequencies  above  200  MHz  though  the  results  appear  reasonable  to  at 
least  500  MHz. 

A  numerical  solution  using  a  pulse  function  basis  results  in  a  single 

value  representing  E  within  each  cell.  It  is  possible  to  use  the  E  values  to 
-*■  ->-* 

calculate  1/2  crE  *  E  for  each  cell  and  use  a  volume  average  to  estimate  the 

g 

specific  absorption  rate  (SAR)  .  Large  numbers  of  cells  must  be  used  in  order 
to  find  accurate  values  of  SAR  by  such  a  procedure.  For  example,  our  calcu¬ 
lations  of  the  SAR  of  a  12cm  muscle  cube  at  1  MHz  show  an  error  of  37%  with 
eight  cells,  26%  with  twenty  seven  cells,  and  20%  with  sixty  four  cells.  The 
delta  functions  used  for  testing  enforce  the  integral  equation  at  the  center 
of  each  cell  so  that  the  calculated  ^  values  are  most  representative  of  the 
cell  centers.  Inspection  of  the  solutions  suggests  that  the  local  values  of 

4 

E  have  appreciably  less  error  than  occurs  in  the  SAR.  If  there  is  much  vari- 

-4  -4 

ation  of  E  within  a  scattering  body  then  even  if  we  have  exact  values  of  E  at 
the  cell  centers  appreciable  error  would  be  expected  in  the  calculated  SAR. 

We  have  found  that  accuracy  is  improved  by  using  a  three-dimensional  inter- 
polant10  with  the  E  values  initially  calculated  for  each  cell  to  account  for 
some  of  the  variation  of  E  within  each  cell.  Trilinear  and  triquadratic 

“4 

interpolants  have  both  been  used  to  estimate  the  variation  of  E  between  the 
cell  centers.  The  interpolant  is  integrated  in  calculating  the  SAR.  For  the 
12cm  muscle  cube  at  1  MHz,  the  SAR  found,  using  twenty  seven  cells  with  the 
triquadratic  interpolant,  has  error  comparable  with  the  calculation  using 
sixty  four  cells  without  the  interpolant.  The  increase  in  cost  due  to  use 
of  the  interpolant  is  about  one  percent.  All  values  of  energy  absorption  in 
this  paper  have  been  calculated  using  interpolants. 


DISTRIBUTION  OF  ABSORBED  ENERGY  —  PART  BODY  RESONANCES 

The  principal  contribution  from  the  free  space  calculations  is  the  dis¬ 
tribution  in  energy  deposition.  Figure  2  shows  the  part-body  and  whole-body 
SAR  for  the  homogeneous  model  of  man  with  two-thirds  the  complex  permittivity 


of  muscle  with  E 


L  and  k  directed  front-to-back.  Calculated  values  of 


a 


FREQUENCY  (MHZ) 


Figure  2. 

Part-body  SAR  for  homogeneous  model  of  man.  Incident  intensity  = 


1  mW.cm' 


whole-body  SAR  are  typically  within  20%  of  values  found  for  prolate-cpheroidal 
and  ellipsoidal  models11.  Unlike  earlier  numerical  models  the  distribution 
of  absorbed  energy  within  the  model  is  in  good  agreement  with  that  found  ex¬ 
perimentally  for  homogeneous  phantom  models12.  For  free  space  irradiation  with 
this  polarization  near  resonance,  the  local  absorption  in  the  legs  and  neck  is 
considerably  higher  than  the  whole-body  average  while  the  torso  has  less  than 
average  absorption. 

When  the  inhomogeneous  complex  permittivities  are  used  with  the  model,  a 
change  of  less  than  two  percent  occurs  typically  in  the  whole-body  SAR,  but  a 
more  significant  change  occurs  in  the  distribution  of  energy  deposition. 

Figures  3  and  4  illustrate  the  distribution  of  absorbed  energy  in  man  at  80 
MHz  in  free  space  for  the  homogeneous  and  inhomogeneous  models,  respectively. 
One  difference  is  that  the  inhomogeneous  model  has  reduced  absorption  in 
regions  with  high  bone  content. 

In  previously  reported  experimental1 3  results  on  whole  body  SAR  for  man, 
an  anomalous  increase  in  the  rate  of  energy  deposition  was  observed  in  the 
region  of  470  MHz  for  k  |j  L  and  E  ||  L  orientations.  This  has  now  been  iden¬ 
tified  as  the  first  (geometrical)  resonance  frequency  of  the  head.  Continu¬ 
ing  experiments  have  given  the  absorption  cross  section  for  the  head  region 
as  large  as  3.1  times  the  physical  cross  section,  perhaps  on  account  of  re¬ 
flections  from  the  nearby  torso.  At  the  head  resonance  frequency,  an  SAR  4.5 
times  the  average  value  for  the  rest  of  the  body  has  been  experimentally  ob¬ 
served. 

Higher  order  resonances  based  on  a  multilayer  formulation  may  exist  for 
the  head  region  at  higher  frequencies.  Initial  numerical  results  show  one 
such  resonance  at  frequencies  of  the  order  of  2000  MHz.  It  is,  however,  an¬ 
ticipated  that  the  overall  absorption  cross  section  at  these  frequencies  may 
not  be  as  large  as  that  for  the  first  resonance,  where  an  enhancement  factor 
of  3.1  has  been  observed  for  the  head  region. 

The  frequency  for  maximum  energy  deposition  in  the  arms  is  approximately 
150  MHz  (see  Fig.  2)  with  an  absorption  cross  section  that  is  2.3  times  the 
physical  cross  section. 


GROUND  EFFECTS 


Experimental  observations  of  ground  and  reflector  effects13  have  been 
reported  previously  but  numerical  methods  have  not  been  available  to  explain 
the  observations.  We  have  used  image  theory  to  reduce  the  problem  of  man 
above  ground  or  in  front  of  reflectors  to  a  multibody  problem  in  free  space. 
Symmetries  have  been  used  to  reduce  the  matrix  size. 

All  calculations  of  ground  effects  have  assumed  that  man  is  standing  on 
or  above  a  perfectly  conducting  ground  plane  that  is  infinite  in  extent.  The 
incident  field  is  vertically  polarized  with  k  directed  front-to-back  on  the 
man.  Internal  fields  in  the  model  are  the  same  as  those  in  one  half  of  a 

double-man  consisting  of  the  model  and  its  image  in  free  space. 

The  resonance  frequency  of  man  standing  on  a  ground  plane  is  one  half 
that  for  man  in  free  space.  The  SAR  of  man  on  the  ground  plane  at  the  re¬ 
duced  resonant  frequency  is  within  2  percent  of  the  SAR  for  man  in  free 

space  at  the  free-space  resonance  frequency. 

At  10  MHz  the  SAR  of  man  standing  on  a  ground  plane  is  0.0163  watts/kg 
per  mW/ cm2 ,  which  is  about  seven  times  that  found  for  the  same  model  in  free 
space.  The  enhancement  in  SAR  due  to  the  ground  effect  is  found  since  the 
frequency  of  10  MHz  is  much  closer  to  the  grounded  man  resonance  frequency 
(“35  MHz)  than  to  the  free-space  resonance  frequency  (“70  MHz) . 

Figures  5  and  6  illustrate  the  distribution  of  absorbed  energy  in  man  at 
10  MHz  in  free  space  and  in  contact  with  the  ground  plane,  respectively.  The 
enhancement  factor  in  local  energy  deposition  due  to  the  ground  effect  is 
about  60  in  the  heel  area. 

Figure  7  shows  the  part  body  and  whole  body  SAR  for  the  inhomogeneous 
model  of  man  with  different  spacings  from  the  ground  plane  at  10  MHz.  A  small 
small  separation  that  breaks  electrical  contact  with  the  ground  plane  is  suf¬ 
ficient  to  eliminate  much  of  the  ground  effect. 
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Figure  6.  2 

Local  SAR  values  (watts/kg  per  mW/cm  )  x  100  for  inhomogeneous  model  of 
man  standing  on  a  ground  plane  at  10  MHz. 
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REFLECTOR  EFFECTS 

All  calculations  of  reflector  effects  have  assumed  that  man  is  standing 
in  front  of  a  reflector  that  is  perfectly  conducting  and  infinite  in  extent. 

•f 

The  incident  field  is  vertically  polarized  with  k  directed  f ront-to-back  on 
the  man.  One  image  is  required  for  a  flat  reflector  and  three  are  needed  for 
a  9C°  corner  reflector. 

At  65  MHz  the  computed  SAR  is  4.87  times  the  free-space  value  when  man 
is  0.1875  X  in  front  of  a  flat  reflector  and  16.6  times  the  free-space  value 
when  man  is  1.5  X  in  front  of  the  axis  of  a  90°  comer  reflector. 

Antenna  theory  may  be  used13  to  calculate  the  ratio  of  effective  area 
of  a  half-wave  dipole  with  a  reflector  to  that  of  the  dipole  in  free  space. 
Such  ratios  are  within  16%  of  the  above  calculated  factors  of  enhancement  of 
SAR  due  to  reflector  effects. 

For  certain  length-to-width  ratios,  the  experimentally1 3  observed  en¬ 
hancement  in  energy  deposition  is  30  to  40%  higher  than  that  anticipated 
from  antenna  theory  and  calculated  numerically  for  reflectors  of  infinite 
dimensions.  This  phenomenon  is,  once  again,  in  agreement  with  experimentally 
obtained  antenna  gains  for  finite  size14  comer  reflector  antennas.  Signi¬ 
ficantly  'enhanced  rates  of  energy  deposition  are  projected  for  all  kinds  of 
comer  angles  (not  just  the  values  corresponding  to  180° /n,  where  n  is  an 
integer)  and  for  reflector  lengths  and  widths  that  are  no  more  than  a  frac¬ 
tion  of  a  wavelength  at  the  resonance  frequency. 

MULTIBODY  EFFECTS 

Driving  point  inpedance  values  for  a  broadside  array  of  two  half-wave 
dipoles15  have  been  used  to  prepare  Fig.  8  which  shows  the  variation  of  ef¬ 
fective  area  per  dipole  with  spacing.  If  the  two  dipoles  are  tangent,  each 
will  receive  approximately  one  half  the.  energy  it  would  receive  if  isolated 
in  free  space.  For  separations  greater  than  two  or  three  free  space  wave¬ 
lengths,  the  antennas  have  little  coupling  so  that  each  has  an  effective  area 
approximating  the  isolated  free  space  value.  At  a  spacing  of  0.65  X,  each 
dipole  will  receive  about  56%  more  energy  than  it  would  if  it  were  isolated. 
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SPACING  IN  WAVELENGTHS 


It  has  been  shown13  that  the  energy  absorbed  by  man  at  resonance  corresponds 
to  the  energy  received  by  a  half-wave  dipole  of  length  approximately  equal  to 
the  height  of  the  man.  Mutual  impedance,  and  hence  driving  point  impedance, 
values  are  dependent  upon  the  shape  of  the  antenna  elements,  but  we  may  ex¬ 
pect  that  the  variation  of  SAR  with  spacing  for  two  men  near  resonance  is 
similar  to  Fig.  8. 

Preliminary  experiments  using  animals  have  confirmed  the  existence  of 

-f.  A 

multibody  effects.  We  have  used  anesthetized  adult  rats  with  E  | |  L  ori¬ 
entation  and  frontal  (broadside)  incidence.  In  tests  with  480  ±  50  gram 
rats  at  600  MHz  with  an  incident  intensity  of  100  mW/cm2,  isolated  animals 
had  an  average  temperature  increase  of  0.675°  C/ min.  For  two  animals  placed 
0.65  X  apart,  the  monitored  animals  had  an  average  temperature  increase  of 
1.14°  C/min.  A  70  percent  increase  in  SAR  was  caused  by  the  presence  of  the 
second  animal. 

Antenna  theory  suggests  similar  enhancements  in  SAR  for  targets  in  the 
subresonance  and  supraresonance  regions,  also  for  spacing  on  the  order  of 
0.65-0.7  X.  Furthermore,  a  greater  enhancement  in  the  SAR  by  a  factor  as 
large  as  2.5  is  anticipated  for  the  inner  targets  with  three  or  more  ele¬ 
ments. 
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Head  Resonance:  Numerical  Solutions 
and  Experimental  Results 
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Abstract — We  have  used  numerical  solutions  and  experiments  with 
phantom  models  of  man.  and  experiments  with  the  Ding  Ivans  rat  to  show 
the  existence  of  head  resonance.  Greatest  absorption  in  the  head  region  of 
man  occurs  at  a  frequency  of  about  375  MHz.  Absorption  is  stronger  for 
wave  propagation  from  head  to  toe  than  It  Is  when  the  electric  field  is 
parallel  to  the  long  axis.  The  highest  absorption  cross  section  for  the 
human  head  Is  projected  to  be  approximately  3.5  times  Its  physical  cross 
section. 

I  I  M  Root  (  I  ION 

Ii  HAVE  previously  reported  numerical  solutions 
for  the  deposition  of  electromagnetic  energy  in  a 
realistic  model  of  man  which  showed  the  existence  of 
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resonances  for  body  parts  such  as  the  head  and  arms,  as 
well  as  for  the  whole  body  [1],  I  urther  yvork  has  been 
done  regarding  head  resonance,  since  we  believe  the  phe¬ 
nomenon  may  be  important  in  the  study  of  behavioral 
effects,  blood-brain  harrier  permeability.  cataractogeneMs. 
and  other  microwave  hioeffects. 

Our  results  show  that  the  Tirst  resonance  of  the  intact 
human  head  occurs  at  a  frequency  of  about  3“'5  Mil/  and 
has  an  .S'  parameter  (ratio  of  absorption  cross  section  to 
physical  cross  section)  of  about  3.5  for  incident  plane 
waves  propagating  from  head  to  toe  Earlier  calculations 
for  a  spherical  model  of  the  isolated  human  head  showed 
a  geometrical  resonance  with  .S' =  I  I  near  450  Mil/,  and  a 
second  resonance  with  .V  =  I  4  occurring  near  2.1  ('dl/ 
when  allowance  is  made  for  the  inhomogeneous  smicture 
by  using  a  multilayered  model  |2|.  The  results  of  our  study 
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'UggcM  thill  the  absorption  in  the  head  legion  !'  ^ t n ■ : i ! » 
dependent  on  the  presence  of  the  test  ..f  the  h.  ■■!» 

In  the  tchovv.ne  sections  o|  tin-  p.tper  we  present  new 
data  regarding  resonance  of  the  nn.ui  him:. in  he  nl  . >t • 
t  une.!  u-:ni:  h . ■ ; ( ;  n  iir.crti.b  '« >i u t ;> > ri'  ..nd  p ! i . i n t ■. ■  i . :  in.  .1 
el'  fi-t  man  .i'  well  a'  e\pei  imental  leMiii  -  Hu  the  I  ,  ne 
I  v  ail'  nit 


II  \l  Ml  KIl  VI  Suit  I  i'  IS- 

I  ie  I  illustrates  a  reah'tie  model  .>1  mm  which  we  a>ed 
in  earlier  computation'  showing  he.u!  tesojiancc  |l|  I  lie 
mode!  w  as  obla ined  bv  iNiit  a  l> - 1 .i !  .  ’ t  I  M  t  cubic.: i  i elK  o| 
\  arti'U'  'i/e'  to  obtain  a  he'!  f . t  ot  the  e ■  >;i r ■  ■  u i  <  a  .It  :- 
er.ini'  of  the  Mlth  pereenlhe  'taiulati!  man  j  -|  Mi  mimeti- 
eai  m'I n tn*ii'  ■  ep» -i  leei  m  till'  paper  ate  t - ■  i  a  m.'.le!  huune 
mere  detailed  modeling  -  'I  tile  neck  anil  he.ui  a-  .It.  m  n  a, 
Fee'  -  and  '  Representation  ot  the  te't  >>!  the  bods  i'  a.' 
'hi'wti  't!  I  te.  1.  Nule  that  the  number  of  .ell-.  Ini'  heel! 
tnere.t'eil  Iren,  12  te  14-1  m  the  he.ui  region.  and  Item  l  to 
22  m  the  neek  region  let  a  whole  bodv  tetal  ei  .'-40  ceils 

We  ha\e  used  a  montent-metheil  solution  ef  the  eleittn 
field  integral  equation  with  a  pulse  function  basi'  and 
delta  funetien.s  let  testing  1 1  j  and  an  ir.tcrpniant  fellon  up 
for  improvement  of  convergence  1 4 j.  Stnee  there  are  three 
unknown  field  components  per  cell,  even  using  one  plane 
of  svmmetrv  of  the  model  the  matrix  is  '10  •  5  It)  complex, 
whivh  correspond'  to  a  'Vstem  of  1020  simultaneous  equa¬ 
tion'  in  1020  real  unknowns.  Noniterative  matrix  solu¬ 
tions  are  required  since  the  matrix  has  at  least  one  com¬ 
plex  eigenvalue  having  a  real  part  with  magnitude  exceed¬ 
ing  tinitv  j 5 j.  A  total  of  17  h  of  dedicated  usage  of  a 
PDP-11)  digital  computer  was  required  for  each  solution. 
1  ortunatelv.  weak  diagonal  dominance  which  is  inherent 
in  the  matrix  formulation  contributes  to  the  conditioning 
to  allow  reasonable  round-off  errors  with  such  a  large 
matrix  si/e. 

Pigs  4  and  5  show  the  head  and  wholc-bodv  specific 
absorption  rate  (S\R)  calculated  for  the  new  model  of 
man  m  free  space  with  an  incident  plane  wave  having  two 
different  polarizations.  Note  that  head  resonance  is  more 
pronounced  when  the  propagation  vector  is  parallel  to  the 
bodv  axis  (A  /  )  with  propagation  from  head  to  toe  than  it 
is  when  the  electric  field  is  parallel  to  the  long  axis  (  /  /  I 

L xpermient.il  data  for  the  I, one  I  vans  rat.  described  i  i 
Section  IV  of  this  paper,  also  suggest  that  head  resonance 
is  stronger  for  A  /  than  for  /  /  polarization.  Some 

experimental  data  obtained  using  biological  phantom 
models  i'f  man  are  presented  for  comparison  in  the  next 
section  of  this  paper 

Pigs  (i  and  7  illustrate  the  distribution  of  absorbed 
energv  in  the  head  region  near  head  resonance.  I  he 
presence  of  a  hot  spot  exists  near  the  center  of  the  head 
which  is  also  seen  in  the  distribution  at  the  geometrical 
resonance  for  a  sphere  (2|.  I  he  distribution  shown  in  Pigs. 
6  and  7  is  onlv  approximate  since  a  homogeneous  model 
was  used  in  all  calculations. 
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Fig  2  Detailed  modeling  of  the  neek  and  head,  front  \:ew 


Fig  3  Detailed  modeling  of  the  neck  and  head,  side  view 
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We  have  used  accurate!'  scaled  figurines  [f>  1  >>f  2<>.?-. 
25  4-.  WO-.  and  406-cni  length  - ; t h  biological  phantom 
mixture^  |7|.  described  ir.  Table  I.  to  model  the  exposure 
of  man  to  plane  waxes  in  free  space  All  xalues  of  ab¬ 
sorbed  dose  were  measured  using  a  Ihermonctics  model 
2401 -A  gradient  laser  calorimeter.  1-ach  figurine  has!  the 
head  attached  to  the  torso  h_x  a  laser  of  salmc-soaked 
cloth,  thus  maintaining  conductixe  contact  but  allowing 
ease  separabihts  for  measurement  of  energx  deposition  in 
the  head,  for  determination  of  the  absorbed  dose  in  the 
head,  after  separation,  the  part  of  the  pole  urethane  mold 
holding  the  biological  phantom  material  shaped  in  the 
form  of  the  head  was  wrapped  in  Saran  Wrap  to  limit  loss 
of  water  bs  evaporation. 

Fable  II  gives  the  experimental  values  of  SAR  for 
figurines  containing  biological  phantom  mixtures  and  ir¬ 
radiated  with  propagation  from  head  to  toe.  A  plot  of 
data  from  Fable  II  on  t  ig  5  shows  good  agreement 
between  the  numerical  solutions  and  experiments  with 
phantom  models  of  man  the  values  in  Table  II  suggest 
that  for  propagation  from  head  to  toe.  the  magnitude  of 
absorption  m  ihe  head  region  near  head  resonance  is  not 
stronglv  dependent  on  orientation  of  the  K  vector 

Both  the  numerical  solutions  in  Section  II  and  the 
phantom  figurines  have  used  homogeneous  models  The 
next  section  of  this  paper  presents  experimental  data 
obtained  with  the  Tong  1  vans  rat 

IV  liXPI  RIMI  NI  M  Rl-SI  I  IS  IIIR  till  I  o\<;  I  s  xss 
Rai 

Table  ill  gives  the  experimental  salues  of  S.AR  for 
several  tests  made  with  both  freshls  killed  and  anesthe¬ 
tized  I  one  1  vans  rats  I  he  number  of  measurements  (n) 
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uni]  c  alculated  standard  errors  in  (Ik-  measurements 
(Sl: VI )  are  also  given  in  the  I  abb  V  liing  o.ija'c'Ms  'it. it 
head  resonance  for  a  medium-sized  r.ii  -In  mid  cum  near 
the  lest  frequence  of  2450  Mil/. 

1  he  anesthetized  rats  were  gixeii  a  4s-mg.kg  dose  of 
sodium  pentohai hitol  to  facilitate  use  of  a  fixed  orienta¬ 
tion  and  limit  thermo-regulators  functions.  Calorimetric 
measurements  were  used  to  determine  the  absorbed  dose 
in  the  freshJs  killed  rats.  These  measurements  were  made 
using  a  I  hermonetic.x  model  2401-A  Seebeek  envelope 
gradient  layer  calorimeter.  I  he  dose  was  also  determined 
b\  measurement  of  the  rate  of  increase  in  rut  colonic 
temperature  and  brain  temperature,  using  iicpud  crystal 
temperature  optical  fiber  (1.(  01  I  probes,  in  the  anesthe¬ 
tized  rats  l  or  brain  temperature  the  1(01  probe  was 
implanted  in  a  triphine  hole  .1  mm  postcrioi  to  the  bregma 
cranial  suture.  2  mm  lateral  to  the  midline  cranial  suture, 
and  6  k  mm  below  the  top  surface  of  the  cortex 

The  results  in  I  able  III  show  lhat  head  resonance  is 
more  pronounced  for  A  /.  than  for  /  /  orientation 

which  was  also  noted  with  the  numerical  solutions  Note 
that  the  ratio  of  hcasl-to-whole-bodx  heat  content  is  less 
for  the  anesthetized  ral  than  for  the  dead  rat  blood 
circulation  mas  he  reducing  the  (clause  magnitude  of 
heating  in  the  head  region. 


\  St  MM.ARX  AND  COM  I  I'SIONS 

Vs  e  base  used  three  methods  to  studs  head  resonance: 
numerical  solutions,  experiments  with  phantom  models  of 
man.  and  experiments  xxith  the  Long  Is  a  n  s  rat.  All  three 
approaches  show  the  existence  of  head  resonance.  The 
phenomenon  appears  to  have  a  greater  magnitude  for 
svase  propagation  from  head  to  toe  than  for  /.  orienta¬ 
tion.  Tor  the  former  orientation  the  head  absorption  cross 
section  as  high  as  T5  times  the  physical  cross  section  is 
projected 

The  numerical  solutions,  xvhich  are  supported  bx  the 
experimental  results  for  phantom  models  of  man.  suggest 
that  the  absorption  is  much  stronger  than  would  be  pre¬ 
dicted  using  a  sphere  to  model  the  isolated  head.  The 
strong  dependence  on  polarization  would,  of  course,  also 
be  missed  using  a  spherical  model. 

We  believe  that  the  enhanced  absorption  in  the  head 
region  mas  make  head  resonance  significant  in  the  studs 
of  behavioral  effects,  blood-brain  barrier  permeability, 
cataractogenesis,  and  other  microwave  bioeffects. 

Ru  t  R!  \<  t  s 

|]|  OH  (iandhi  and  M  J  Hagmann.  "Some  recent  results  on  deposi 
tion  ot  electromagnetic  energs  in  animals  and  in  minlels  of  man.1 
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Electromagnetic  Absorption  in  a  Multilayered 

Model  of  Man 


PETER  W,  BARBER,  member,  ieek.OM  P.  GANDHI,  ee.li.ow.  ieee.  MARK  J.  HAGMANN, 
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Abstract -A  multilayered  planar  model  is  used  to  examine  the  depen¬ 
dence  of  whole-body  power  absorption  on  the  configuration  of  surface 
layers,  e.g„  skin,  fat,  and  muscle  which  normally  occur  in  biological 
bodies.  It  is  found  that  the  layering  resonance  for  three-dimensional 
bodies  (as  opposed  to  the  geometrical  resonance)  can  be  predicted  quite 
accurately  by  a  planar  model.  Calculations  for  a  multilayered  prolate 
spheroidal  model  of  man  predict  a  whole-body  layering  resonance  at 
1.8  GHz  with  a  power  absorption  34  percent  greater  than  that  pre¬ 
dicted  by  a  homogeneous  model. 

Introduction 

ECENT  interest  in  quantifying  both  the  hazardous  and 
potentially  beneficial  effects  of  nonionizing  electromag¬ 
netic  (EM)  radiation  on  man  has  been  the  impetus  for  a  great 
deal  of  experimental  and  theroretical  research.  Of  particular 
interest  are  calculations  to  determine  the  relationship  between 
incident  power  density  and  the  resulting  absorbed  power  due 
to  whole-body  irradiation  of  man.  Recent  theoretical  methods 
that  have  been  used  include  a  perturbation  approach  ( 1 1  useful 
for  analyzing  homogeneous  ellipsoidal  models  and  an  integral 
equation  method  [2)  for  analyzing  homogeneous  prolate 
spheroidal  models.  The  absorption  characteristics  of  inhomo¬ 
geneous  irregular  shaped  models  constructed  of  cubical  sub¬ 
volumes  have  been  studied  by  volume  integral  moment  method 
techniques  [3],  [4).  In  these  whole-body  irradiation  studies, 
the  power  absorption  is  determined  as  a  function  of  the  angle 
of  incidence,  polarization,  and  frequency  of  the  incident  EM 
wave.  Usually  the  results  are  shown  as  a  plot  of  specific  ab¬ 
sorption  rate  in  watts/kilogram  versus  frequency,  where  an 
incident  power  density  of  I  inW/cm 7  is  assumed.  Typically,  for 
a  given  angle  of  incidence  and  polarization,  the  power  absorp¬ 
tion  increases  relatively  rapidly  with  frequency  to  a  resonant 
peak  and  then  slowly  decreases  to  an  asymptotic  high-frequency 
limit. 

There  is  another  interaction  which  takes  place  which  has  not 
been  considered  in  recent  whole-body  calculations.  Both 
homogeneous  and  inhomogeneous  models  have  not  accounted 
for  the  EM  interaction  due  to  the  surface  layering  of  biological 
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bodies,  eg.,  the  layers  of  skin.  fat.  and  muscle.  An  investigation 
of  these  effects  is  the  subject  of  this  paper. 

Absorption  L elects  Due  to  Layering 

Early  calculations  to  assess  the  biological  significance  of  EM 
radiation  used  planar  models  consisting  of  skin.  fat.  and  semi¬ 
infinite  muscle  layers  [ 5 1  Later  a  two-lavered (fat  and  muscle) 
spherical  model  was  used  |6] .  These  earlier  investigations  are 
summarized  in  (7|.  Both  studies  showed  that  the  absorption 
characteristics  are  critically  dependeni  on  the  number  and 
thickness  of  the  surface  layers.  More  recent  whole-body  calcu¬ 
lations  have  not  considered  the  effect  of  surface  layers.  One- 
reason  for  this  is  due  to  the  unavailability  of  multilayered  ana¬ 
lytical  solutions  for  any  geometry  other  than  the  plane,  infinite 
circular  or  elliptical  cylinder,  and  sphere.  In  the  case  of  nu¬ 
merical  approaches,  such  as  the  volume  moment  method 
technique,  practical  limitations  on  the  number  of  cells  which 
may  be  used  prohibit  consideration  of  layering  effects. 

A  six-layered  sphere  has  recently  been  used  to  model  the  iso¬ 
lated  head  ( 8 1  - ( 1 0 1  .  Of  particular  interest  are  results  in  (9] 
and  [10]  which  compare  the  frequency  dependent  absorption 
in  stx-layered  spheres  with  that  in  corresponding  homogeneous 
spheres.  One  of  the  power  absorption  results  in  |9|  has  been 
recalculated  and  is  shown  in  Eig.  1  as  a  plot  of  absorption 
efficiency  versus  frequency.  This  and  subsequent  multilayered 
sphere  calculations  were  made  using  a  Mie  theory  computer 
program  based  on  tl  ■  mathematical  development  in  |8|  .  The 
absorption  eltici.  s  the  total  power  absorbed  divided  by 
the  powe--  .Rj  .  ,  the  geometrical  cross  section.  An 
absorpti<-  m  iency  ater  than  unity  indicates  that  the  in¬ 
fluence  of  the  body  on  the  incident  wave  extends  beyond  the 
geometrical  boundary.  The  absorption  efficiency  for  the  homo¬ 
geneous  bram  tissue  sphere  increases  with  frequency  until  the 
maximum  is  reached  and  then  slowly  decreases.  The  amplitude 
and  frequency  of  the  resonance  are  dependent  upon  the  size, 
shape,  and  dielectric  characteristics  of  the  model.  This  max¬ 
imum  absorption  condition  can  be  called  the  geometrical 
resonance.  The  absorption  efficiency  for  the  six-layered  model 
has  an  additional  resonance  which  is  due  to  the  impedance 
matching  effects  of  the  surface  layers  and  at  this  resonance 
the  absorption  efficiency  is  almost  30  percent  greater  than  that 
at  the  geometrical  resonance.  This  second  maximum  absorption 
condition  can  be  called  the  layering  resonance. 
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Fig.  2.  Absorption  characteristics  of  a  six-layer  plana:  lodel.  The  five 
outside  layers  are  the  same  as  those  in  Fig.  1 .  The  innermost  layer  is 
semiinfinite  brain  material. 


FREQUENCY, GHZ 

Fig.  1.  Absorption  characteristics  of  the  six-layered  10-cm  radius  spher¬ 
ical  head  model  of  Joines  and  Spiegel  (9|.  The  five  outside  layers 
(and  thickness)  are  skin  (0.15  cm),  fat  (0.12  cm),  bone  (0.43  cm), 
dura  (0.1  cm),  andCSF  (0.3  cm).  The  radius  of  the  inner  brain  sphere 
is  8.9  cm. 

The  most  interesting  feature  shown  in  Fig.  I  is  that  the  pres¬ 
ence  of  layering  enhances  the  absorption  of  the  spherical  model. 
Another  feature  is  that  the  geometrical  and  layering  resonances 
appear  to  be  independent  of  one  another,  i.e.,  the  geometrical 
resonances  of  the  homogeneous  and  layered  spheres  are  almost 
identical,  while  the  layering  resonance  of  the  six-layered  sphere 
appears  merely  as  an  enhancement  of  the  absorption  over  that 
of  the  homogeneous  sphere. 

Calculations  for  a  homogeneous  and  six-layered  planar  model 
have  been  made  for  the  same  tissue  thicknesses  and  dielectric 
characteristics  as  were  used  to  obtain  the  results  in  Fig.  1 .  The 
calculations  were  made  using  conventional  planar  techniques 
[11].  The  brain  tissue  core  in  the  spherical  model  was  assumed 
to  extend  to  infinity  in  the  planar  model.  These  results  are 
shown  in  Fig.  2.  The  planar  model,  which  does  not  exhibit  a 
geometrical  resonance,  does  however  show  the  same  layering 
resonance  as  is  observed  in  Fig.  1  for  the  layered  spherical 
model.  The  location  of  the  layering  resonance  and  the  en¬ 
hancement  of  the  absorption  efficiency  due  to  the  layering  are 
almost  identical  for  the  planar  and  spherical  models.  Taking 
the  ratio  of  the  layered  to  homogeneous  absorption  efficiencies 
for  the  planar  model  in  Fig.  2  pnd  multiplying  by  the  homo¬ 
geneous  sphere  result  in  Fig.  1  gives  the  dashed  curve  in  Fig.  3, 
a  predicted  absorption  for  the  layered  sphere.  Comparing  this 
predicted  absorption  to  the  actual  absorption  shows  that  the 
layering  enhancement  of  the  planar  model  can  be  applied  to 
the  homogeneous  sphere  solution  to  predict  the  absorption  in 
the  layered  sphere  to  within  5  percent,  both  in  absorption  ef¬ 
ficiency  at  a  given  frequency  and  in  the  location  of  the  reso¬ 
nant  peaks. 
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Fig.  3.  Actual  and  planar  predicted  absorption  characteristics  of  the  six¬ 
layered  spherical  head  model  of  Fig.  1 . 


There  are  two  obvious  questions  that  neeu  to  be  answered: 
Why  are  the  geometrical  and  layering  resonances  essentially 
independent  in  Fig.  1 ,  and  why  does  the  planar  model  predict 
the  layering  resonance  of  the  head  model  with  such  accuracy? 

To  answer  the  first  question,  consider  the  geometrical  reso¬ 
nance  of  the  head  model  in  Fig.  I .  The  resonance  of  the  homo¬ 
geneous  sphere  is  dependent  orr  the  radius  and  constitutive 
makeup  of  the  sphere.  For  high-loss  dielectric  bodies  of  the 
type  considered  here,  the  peak  absorption  occurs  when  the 
ratio  of  circumference  to  free-space  wavelength  (2 na/\  =  ka) 
is  about  unity.  For  the  homogeneous  sphere  in  Fig.  1,  reso¬ 
nance  occurs  when  ka  =s  0.94.  Referring  to  the  scale  drawing 
of  the  sphere  in  Fig.  1 ,  it  is  clear  that  the  surface  layers  repre- 
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sent  a  small  fraction  of  the  total  sphere  radius  (about  10  per¬ 
cent).  One  would  expect  then  that  replacing  the  outer  portion 
of  the  homogeneous  sphere  by  a  layered  segment  which  on  the 
average  has  similar  dielectric  characteristics  and  whose  thick¬ 
ness  is  much  less  than  a  wavelength  in  the  material  will  have 
little  effect  on  the  geometrical  resonant  frequency.  This  is  the 
case  as  noted  in  Fig.  1.  The  layers  have  an  independent  reso¬ 
nant  frequency,  which  in  this  case  is  approximately  four  times 
the  geometrical  resonant  frequency.  The  frequency  at  which 
the  layers  resonate  is  a  function  of  their  thickness  and  constit¬ 
utive  parameters.  Given  the  electrical  characteristics  of  the 
layers,  the  layering  resonant  frequency  is  inversely  proportional 
to  the  thickness  of  the  layers.  The  fact  that  the  layering  reso¬ 
nant  frequency  is  well  removed  from  the  geometrical  resonant 
frequency  of  the  homogeneous  sphere  accounts  for  the  inde¬ 
pendence  of  the  two  resonances.  Numerical  calculations  show 
that  when  the  layers  are  a  larger  fraction  of  the  sphere  radius, 
i.e.,  the  frequency  of  surface  layer  resonance  approaches  the 
geometric  resonance  frequency  of  the  homogeneous  sphere, 
then  the  geometrical  resonance  of  the  inhomogeneous  sphere 
is  no  longer  the  same  as  the  geometrical  resonance  of  the 
homogeneous  sphere. 

Now  consider  the  second  question  as  to  why  the  planar 
models  can  predict  both  the  resonant  frequency  and  enhance¬ 
ment  over  the  homogeneous  case  for  the  layered  sphere. 
Looking  first  at  the  resonant  frequency,  we  note  that  in  the 
planar  case  the  incident  and  transmitted  waves  are  normal  to 
the  surface.  In  the  case  of  the  sphere,  the  only  statement  that 
can  be  made  with  certainty  is  that  the  incident  wave  is  in  gen¬ 
eral  not  normal  to  the  surface.  However,  the  absorption  is  due 
to  an  internal  interaction  and  the  behavior  of  the  fields  trans¬ 
mitted  into  the  sphere  must  be  considered. 

Snell’s  law  can  be  used  to  relate  the  angle  of  the  transmitted 
wave  to  the  angle  of  the  incident  wave  at  local  regions  on  the 
surface.  Snell’s  law  for  a  wave  propagating  from  free  space 
with  permittivity  e0  into  a  region  with  permittivity  e  is  given 
by 


sin  0,  =  (e0/e)l/2  sin  9 , 

where  0,  and  9,  are  the  angles  of  incidence  and  transmission, 
respectively,  both  angles  being  measured  from  the  local  normal 
to  the  surface.  The  dielectric  characteristics  of  biological  tis¬ 
sue  are  characterized  by  a  relatively  large  complex  permittivity. 
For  example,  at  the  layering  resonant  frequency  of  2.06  GHz, 
the  dielectric  constant  of  the  outer  sk'.,  layer  is  47.5-j  1 1 .4. 
Substituting  this  into  Snell's  law  for  0,  =  90°  (grazing  inci¬ 
dence).  it  is  found  that  the  angle  of  transmission  is  only  8°  from 
the  normal.  Since  the  angle  of  transmission  will  be  even  smaller 
for  other  angles  of  incidence,  it  is  clear  that  the  transmitted 
wave  propagates  almost  normally  into  the  tissue  regardless  of 
the  angle  of  incidence,  i.e.,  the  wave  transmitted  into  the 
sphere  interacts  with  the  surface  layers  in  essentially  the  same 
manner  as  it  does  in  the  planar  model,  as  shown  in  Fig.  4. 
Therefore,  the  resonant  frequency  resulting  from  this  inter¬ 
action  should  be  the  same  whether  the  layers  reside  on  the  sur¬ 
face  of  a  sphere  or  on  the  surface  of  a  plane. 


I- i|i.  4.  The  local  behavior  of  the  transmitted  wave  is  almost  identical, 
i.e.,  normal  to  the  surface,  at  (a)  planar  and  (b)  nonplanar  interfaces 
for  materials  with  a  high  dielectric  constant. 


The  fractional  enhancement  of  the  absorption  due  to  the 
layering  can  best  be  understood  by  considering  the  surface 
layers  as  a  frequency-dependent  impedance  matching  device. 
At  frequencies  well  removed  from  the  layering  resonant  fre¬ 
quency  (either  far  above  or  far  below),  the  surface  impedance 
of  the  layered  model  is  the  same  as  for  the  homogeneous 
model,  i.e.,  at  the  lower  frequencies  the  layers  are  so  thin  as 
to  have  a  negligible  effect  and  at  very  high  frequencies  the 
depth  of  penetration  is  so  low  that  the  transmitted  power  is 
all  absorbed  m  the  surface  skin  layer  which  has  electrical  char¬ 
acteristics  almost  identical  to  those  of  the  brain  material  of 
the  homogeneous  model.  In  the  region  of  the  layering  reso¬ 
nant  frequency,  an  enhanced  power  transmission  occurs,  and 
this  enhancement  is  dependent  only  on  the  surface  layering 
configuration,  and  therefore  the  fractional  enhancement  in 
power  transmission  will  be  the  same  for  both  spherical  and 
planar  models.  In  the  layered  planar  model,  the  power  trans¬ 
mitted  into  the  innermost  layer  is  completely  absorbed.  The 
fractional  power  absorption  enhancement  in  the  sphere  can 
only  be  the  same  if  all  the  transmitted  power  is  also  absorbed. 
This  will  occur  if  the  diameter  of  the  innermost  brain  material 
is  many  skin  depths  thick.  A  calculation  at  the  layering  resonant 
frequency  of  2.06  GHz  shows  that  the  diameter  of  the  inner 
brain  sphere  is  greater  than  six  skin  depths,  insuring  almost 
total  absorption. 

The  planar-prediction  procedure  which  has  been  outlined 
here  for  the  spherical  head  model  is  equally  applicable  to  arbi¬ 
trary  nonspherical  bodies,  because  the  mechanism  of  absorp¬ 
tion  enhancement  depends  only  on  the  transmitted  wave  being 
nearly  perpendicular  to  the  surface  layers,  a  phenomenon  which 
is  independent  of  the  body  shape,  hxtensive  calculations  for 
other  spherical  bodies  and  infinite  cylinders  (II  Massoudi,  pri¬ 
vate  communication)  have  verified  the  validity  of  the  planar- 
prediction  method.  Summarizing: 

1 )  The  layering  and  geometrical  absorption  resonances  for  a 
layered  three-dimensional  object  will  be  independent  if  the 
layering  resonant  frequency  is  well  above  the  geometrical 
resonance,  i.e.,  the  surface  layers  responsible  for  the  layering 
resonance  are  a  small  fraction  of  the  overall  size  of  the  object. 

2)  If  the  resonances  are  independent  (a  condition  which  can 
be  tested  by  comparing  the  geometric  resonance  of  the  homo¬ 
geneous  three-dimensional  object  with  the  layering  resonance 
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of  a  semiinfmite  planar  model),  then  a  planar  model  can  be 
used  to  predict  the  absorption  layering  resonance  for  non- 
planar  geometries  provided  that  the  complex  dielectric  constant 
of  the  object  is  large  so  that  the  transmitted  fields  are  normal 
to  the  layers  and  the  skin  depth  in  the  material  is  small  enough 
that  the  power  transmitted  into  the  interior  of  the  three- 
dimensional  object  can  be  assumed  to  be  completely  absorbed. 
Note  that  1)  must  be  satisfied  before  2)  can  be  applied,  i.e., 
the  geometric  resonance  must  be  independent  of  the  layering 
because  the  planar  model  cannot  predict  changes  in  the  geo¬ 
metric  resonance,  only  in  the  layering  resonance.  Conditions 
1 )  and  2)  may  appear  to  be  so  restrictive  as  to  have  little  appli¬ 
cation.  however,  they  are  in  fact  satisfied  by  most  biological 
models. 


Layering  Absorption  Resonance  in  Man 

It  has  been  shown  that  a  semiinfinite  planar  model  can  accu¬ 
rately  predict  the  layering  resonance  in  a  nonplanar  biological 
model.  The  interest  here  is  in  determining  the  effect  of  layer¬ 
ing  on  the  whole-body  absorption  of  EM  waves  by  man.  This 
will  be  done  by  using  a  planar  model  to  determine  the  absorp¬ 
tion  enhancement  due  to  surface  layering  and  then  applying 
the  resulting  enhancement  factor  to  whole-body  absorption 
results  previously  obtained  for  a  homogeneous  prolate  spher¬ 
oidal  model  of  man. 

In  the  numerical  calculations  which  follow,  the  dielectric 
characteristics  of  the  various  tissue  types  are  taken  front  (12) 
for  frequencies  up  to  and  including  10  GHz.  Above  10 GHz, 
the  permittivity  and  conductivity  are  based  on  the  characteris¬ 
tics  of  electrically  polarizable  molecules  [10],  modified  to 
provide  continuity  at  10  GHz  with  the  lower  frequency  dielec¬ 
tric  characteristics. 

The  surface  layers  of  man  in  general  consist  of  skin-fat- 
muscle  or  skin-fat-muscle-bone-muscle  arrangements.  The 
surface  layering  information  required  for  the  multilayered 
planar  model  was  obtained  from  published  anatomical  cross- 
sectional  data  [13] .  [14] .  Tissue  thicknesses  were  examined  in 
79  horizontal  cross  sections  of  man.  Average  thickness  of  the 
surface  layers  of  skin,  fat,  muscle,  and  bone  were  calculated 
over  the  front  half  of  each  cross  section.  It  was  found  that  the 
surface  characteristics  could  be  represented  by  three-layer 
skin-fat-muscle  configurations  in  37  of  the  cross  sections, 
while  the  remaining  42  cross  sections  were  better  represented 
by  a  five-layer  skin-fat-muscle-bone-muscle  configuration. 
The  mean  thickness  (and  standard  deviation)  of  the  skin  and 
fat  layers  for  79  cross  sections  are  2.25  mm  (0.73mm)  and 
2.62  mm  (1.47  mm),  respectively.  Corresponding  values  for 
the  muscle  and  hone  layers  for  42  cross  sections  are  9.74  mm 
(4.32  mm)  and  9.01  mm  (4.81  mm),  respectively.  Calculations 
similar  to  those  shown  in  Fig.  2  were  made  for  each  of  the  79 
cross  sections.  The  results  were  weighted  by  the  fractional 
frontal  surface  area  which  they  represented  and  then  averaged 
to  obtain  the  layered  curve  in  Fig,  5.  Note  that  there  is  an 
absorption  enhancement  over  homogeneous  muscle  tissue  of 
34  percent  at  1.8  GHz  and  a  reduction  of  24  percent  at  5.7 
GHz 
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I-'ip .  5.  Predicted  whole-body  layering  resonance  in  man  utilizing  a 
planar  model. 


f  ig.  6.  Whole-body  absorption  versus  frequency  for  a  model  of  man 
exposed  to  a  broadside  incident  plane  wave  polarized  parallel  to  the 
long  axis  of  the  model.  The  dashed  curve  was  calculated  using  homo¬ 
geneous  prolate  spheroidal  and  circular  cylindrical  models  (o/ft  = 
6.34,  a  -  0.875  m).  The  solid  curve  is  the  planar-predicted  absorp¬ 
tion  in  a  corresponding  multilayered  model  of  man. 

Fig.  6  shows  the  whole-body  absorption  efficiency  for  both 
a  homogeneous  and  a  multilayered  model  of  man.  The  homo¬ 
geneous  result  was  obtained  using  a  combination  of  prolate 
spheroidal  and  cylindrical  models  ]  1 5] .  The  specific  result 
shown  is  for  the  maximum  absorption  case  which  occurs  when 
the  model  is  illuminated  broadside  by  a  wave  polarized  parallel 
to  the  long  axis  of  the  model.  The  homogeneous  results  in 
[15]  were  obtained  using  the  dielectric  characteristics  for  a 
homogenized  mixture  of  fat.  bone,  and  muscle  rather  than 
muscle  tissue  alone.  Homogeneous  calculations  for  muscle  tis¬ 
sue  alone  would  result  in  a  more  accurate  multilayered  predic¬ 
tion;  however,  these  results  are  not  presently  available.  The 
multilayered  result  was  obtained  by  taking  the  ratio  of  the 
layered  to  homogeneous  absorption  efficiencies  for  the  planar 
model  in  Fig.  5  and  multiplying  the  layering  enhancement  factor 
by  the  homogeneous  result  in  Fig.  6.  Inasmuch  as  the  multi¬ 
layered  data  used  to  generate  Fig.  5  were  obtained  over  the 
front  surface  of  the  cross  sections  of  man.  the  multilayered 
curve  in  Fig.  6  is  applicable  to  frontal  illumination  of  a  man 
model.  The  net  effect  of  the  surface  layering  is  generally  to 
increase  the  absorption  in  the  postresonance  region.  The 
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Lisciinji  ciih jilt. oiiKMi t  ji  I  S  (ill/  and  rcductnm  at  5.7  (ill/ 
represent.',  the  eieaiesi  denaiinn  lioin  the  Imnnipeneous  result 
Calculations  hesoiul  10  (ill/  shoes  that  the  two  results  metes 
when  the  depth  of  penetration  becomes  so  small  that  the  total 
absoiption  in  the  multilayered  model  occurs  in  the  muscle-like 
skin  tissue. 

Comparini:  Fie  <1  to  Fig.  1 .  it  Is  clear  that  the  enhancement 
of  the  absorption  due  to  layering  is  not  nearly  as  great  for  the 
whole-body  model  as  it  is  for  the  isolated  head  model.  There 
are  two  reasons  for  this.  First,  for  the  polarization  of  the  inci¬ 
dent  wave  shown  here,  the  long  slender  shape  of  the  prolate 
spheroidal  whole-body  model  results  in  much  greater  absorp¬ 
tion  at  the  geometrical  resonance  (absorption  efficiency  =  4.43) 
as  compared  to  the  spherical  head  model  (absorption  effi¬ 
ciency  =  1.08)  Second,  the  layering  enhancement  for  the 
whole-body  case  is  much  less  than  for  the  spherical  head 
model  because  of  the  averaging  over  many  different  layering 
combinations  in  the  whole-body  calculations  rather  than  using 
a  single  layering  configuration  as  was  done  for  the  spherical 
head  model.  If  an  averaging  scheme  were  used  in  the  head 
model  calculations,  a  reduction  in  the  layering  enchancement 
would  occur  there  also. 


Conclusions 

It  has  been  shown  that  a  planar  model  can  accurately  predict 
the  layering  resonance  for  a  nonplanar  geometry.  Specifically, 
results  for  multilayered  and  homogeneous  semiinfinite  planar 
models  determine  a  layering  enhancement  factor  which  can  be 
applied  to  whole-body  absorption  results  obtained  from  non- 
layered  three-dimensional  geometries  to  predict  the  absorp¬ 
tion  characteristics  of  three-dimensional  layered  geometries. 
The  technique  has  been  used  to  predict  the  power  absorption 
characteristics  of  a  multilayered  model  of  man. 
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ABSTRACT 

The  factors  of  microwave  radiation  frequency  and  orientation  of  an  animal 
in  the  microwave  beam  are  important  parameters  in  determining  dose  rate.  Dur¬ 
ing  microwave  irradiation  these  factors  may,  on  the  one  hand,  account  for  dif¬ 
ferential  temperature  rise  in  the  rat  and,  on  the  other  hand,  produce  differen¬ 
tial  disruption  of  behavioral  performance.  Fifteen  male  Long  Evans  rats,  in 
Experiment  I,  were  exposed  to  continuous  wave  (CW)  microwave  energy  in  the  220- 
500  MHz  region  in  three  body  orientations  to  determine  the  most  absorbant  fre¬ 
quency  and  body  orientation  within  the  microwave  field  as  evidenced  by  colonic 
temperature  rise.  The  most  absorbant  orientation  was  found  with  the  animal 
length  aligned  parallel  to  the  electric  field.  In  Experiment  II,  six  male  Long 
Evans  rats  were  trained  to  lever  press  for  dry  food  pellets.  Each  animal  was 
then  randomly  exposed  while  lever  pressing  to  three  microwave  frequencies  (360, 
480,  and  50 0  MHz)  in  two  body  orientations  (parallel  to  the  electric  field  or 
parallel  to  the  direction  of  microwave  travel).  The  differential  absorption 
rates  due  to  the  systematic  manipulation  of  frequency  and  orientation  variables 
did  selectively  disrupt  behavioral  performance  of  the  rat.  Greatest  disruption 
of  behavioral  performance  occurred  when  exposures  were  made  with  the  animal 
aligned  parallel  to  the  electric  field  at  500  MHz.  Less  disruption  of  perfor¬ 
mance  occurred  at  360  or  480  MHz  in  this  orientation.  Little  disruption  of 
performance  was  observed  when  the  animal  was  aligned  parallel  to  the  direction 
of  the  microwave  travel  at  360,  480,  or  500  MHz. 


INTRODUCTION 

In  recent  years  the  relationship  between  electromagnetic  phenomena  and 
biological  systems  has  generated  considerable  interest.  Of  particular  interest 
have  been  the  variables  which  determine  the  dose  rate  on  the  biological  target. 
A  variety  of  investigators  [1,  2,  3,  4]  have  noted  that  the  dose  rate  for  the 
biological  organism  depends  strongly  on  such  factors  as:  (1)  average  and 
instantaneous  field  intensity,  (2)  duration  of  exposure,  (3)  size  and  shape  of 
the  organism,  (4)  dielectric  properties  of  the  organism,  (5)  wavelength  of 
radiation,  and  (6)  orientation  of  the  organism  in  microwave  fields.  Recently 
Gandhi  [5,  6]  has  noted  a  strong  resonance  phenomena  based  on  animal  size, 
wavelength  of  radiation,  and  animal  orientation  in  microwave  fields.  Since 
these  factors  should  account  for  differential  absorption  and  temperature  rise 

Hiolo'jical  FJffeobB  of  Elccivoma<mct.i  '  .V.'  V;!,  Vo]  .  1,  Selected  Papers 
of  the  1975  USNC/URSI  Meeting,  Boulder,  Colorado  [HEW  Publication  (FDA) 
77-8011] 


in  the  rat  during  irradiation,  these  conditions  should  also  produce  differential 
disruption  of  behavioral  performance  in  the  rat  during  irradiation.  To  investi¬ 
gate  these  factors  further,  two  experiments  were  conducted  to  determine  the 
effect  of  these  factors  on  rat  heating  and  behavioral  performance  during  irradi¬ 
ation. 


EXPERIMENT  I:  MATERIALS  AND  METHOD 
Subjects 

Fifteen  male  Long  Evans  rats  350-3S0  grams  in  body  weight  were  used.  All 
rats  were  deprived  of  food  to  85  percent  of  their  free-feeding  weights. 

Microwave  Irradiation  Apparatus 

The  microwave  radiation  chamber  consisted  of  two  parallel  copper  plates. 

The  lower  plate  was  2  mm  in  thickness  and  measured  120  cm  x  198.2  cm.  The 
upper  plate  was  1  mm  in  thickness  and  measured  63.5  cm  x  198.2  cm.  Both  plates 
were  enclosed  in  a  rectangular  plywood  box  215  cm  x  88  cm  x  150  cm.  The  inte¬ 
rior  of  the  box  was  lined  with  sheets  of  microwave  absorbing  material  (11.8  cm 
thick  Eccosorb).  A  ventilating  fan  was  placed  at  one  end  of  the  chamber  to 
remove  air  from  the  chamber  at  the  rate  of  310  cubic  feet  per  minute.  Addi¬ 
tional  ventilation  openings  were  placed  near  the  position  of  the  animal  to 
allow  entry  of  fresh  air.  Two  small  doors  (63  x  63  cm)  were  constructed  in 
each  side  of  the  chamber. 

A  block  diagram  of  the  exposure  apparatus  is  shown  in  Fig.  1.  Microwave 
power  was  obtained  from  an  Epsco  signal  source  operating  in  the  200-500  MHz 
band  with  a  minimum  output  of  100  watts  power.  A  Philco  model  164B  bidirec¬ 
tional  power  meter  was  used  to  monitor  power  delivered  to  and  reflected  from  the 
parallel  plate  chamber.  Power  leaving  the  parallel  plate  radiation  chamber  was 
monitored  by  a  Hewlett-Packard  No.  434A  calorimetric  power  meter.  A  Heathkit 
cantenna  was  used  for  final  dissipation  of  power  from  the  parallel  plate  cham¬ 
ber.  A  Weinschell  No.  109L  double  stub  tuner  was  used  to  balance  the  irradia¬ 
tion  system  and  reduce  reflected  power  to  a  minimum.  Frequency  of  the  micro- 
wave  energy  delivered  to  the  parallel  plate  chamber  was  monitored  by  a  Hewlett- 
Packard  No.  5253B  frequency  counter.  Interconnection  of  the  various  components 
of  the  irradiation  system  utilized  UG58/U  coaxial  cable. 

Determi nation  of  Absorbed  Microwave  Power 

The  Intensity  of  the  microwave  field  within  the  parallel  place  chamber  at 
the  position  of  the  rat  was  measured  as  mW/cm^  by  a  General  Microwave  Corpora¬ 
tion  Raham  model  1  (300-18,000  MHz)  field  intensity  probe. 

A  calorimetric  method  [2]  was  used  to  determine  how  much  of  the  microwave 
power  was  actually  absorbed  by  the  animal  during  each  exposure.  The  amount  of 
absorbed  power  was  expressed  as 

AT 

mW/gm  =  4180  x  — 


where 


AT  =  rectal  temperature  increase  of  rat  in  degrees  Celsius  during 
microwave  exposure 


t  =  time  of  exposure  in  seconds 
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4180  =  conversion  constant 


Procedure 

Each  rat  was  habituated  for  thirty  minutes  a  day  over  a  four-day  period  to 
a  cylindrical  Plexiglas  rat  holder  (.32  cm  x  6.5  cm  x  24  cm).  The  holder  was 
constructed  with  1.27  cm  slots  on  both  upper  and  lower  surfaces  to  interfere 
little  with  normal  heat  dissipation  by  the  animal  and  to  ensure  release  of  urine 
and  fecal  material.  On  exposure  days  a  microwave  transparent  liquid  crystal 
fiber  optic  temperature  probe  [7]  was  inserted  colonically  (4  cm)  and  held  in 
place  by  masking  tape  applied  to  the  rat's  tail.  Each  rat  was  then  placed  in 
the  rat  holder  and  set  inside  the  radiation  chamber  such  that  the  length  of  the 
rat  body  was  parallel  to  the  electric  component  of  the  microwave  field  (E 1 1  L) . 
Each  rat  was  allowed  approximately  twenty  minutes  before  microwave  exposure  to 
allow  for  colonic  thermal  equilibrium.  Exposure  to  microwave  irradiation  com¬ 
menced  when  colonic  temperature  had  stabilized.  In  each  case  microwave  power 
was  gradually  increased  over  a  5-second  period  to  the  established  field  density 
of  25  raW/cm^.  This  was  done  to  avoid  peak  surges  of  power  above  25  mW/cm2  at 
radiation  onset.  Each  rat  was  then  exposed  to  ten  minutes  of  microwave  radia¬ 
tion  at  one  of  twenty  frequencies  in  the  220-500  MHz  range.  The  microwave  field 
intensity,  raW/cm^,  and  amount  of  power  absorbed  by  the  animal,  mW/gm,  were 
recorded  for  each  exposure.  All  rats  were  exposed  randomly  to  the  twenty  fre¬ 
quencies.  This  procedure  was  repeated  four  times.  Some  of  the  rats  were  addi¬ 
tionally  exposed  to  microwave  radiation  at  360,  480,  and  500  MHz  with  the 
length  of  the  rat  body  parallel  to  the  magnetic  field  (H 1 1  L)  and  the  direction 
of  microwave  propagation  through  the  chamber  (C||  L).  Relative  humidity  during 
the  radiation  exposures  varied  between  12  percent  and  40  percent.  Ambient 
temperature  of  the  radiation  chamber  was  controlled  by  room  temperature  (70-72° 
F).  All  exposures  were  CW  radiation. 

EXPERIMENT  I:  RESULTS 

Figure  2  shows  the  mean  relative  absorption  of  microwave  power  ((mW/gm)/ 
(mW/cm^))  and  a  scale  indicating  the  animal  rise  in  body  temperature  (centigrade) 
per  minute  from  220-500  MHz  for  animals  aligned  with  the  long  axis  of  the  body 
parallel  to  the  electric  field  (e||  L).  Also  shown  are  the  mean  ((mW/gm) /(mW/ 
cm2))  values  of  absorbed  power  for  animals  aligned  parallel  to  either  the  mag¬ 
netic  field  (H ! |  L)  or  with  the  direction  of  wave  propagation  through  the  cham¬ 
ber  (C 1 1  L)  at  360,  480,  and  500  MHz.  The  most  microwave  power-absorbant 
animal  orientation  in  this  frequency  range  clearly  was  with  the  long  axis  of  the 
animal  parallel  to  the  electric  field  (E||  L).  In  this  orientation  500  MHz 
appears  to  be  the  most  absorbant  frequency.  Analysis  of  variance  techniques 
(8]  indicates  that  there  is  a  differential  absorption  of  power  in  the  electric 
orientation  depending  on  frequency  (£  <  .01).  Use  of  Newman  Kuels  multiple 
comparison  techniques  further  indicates  that  exposures  to  500  MHz  irradiation 
in  the  electric  orientation  were  clearly  different  from  all  other  exposure  fre¬ 
quencies  in  this  orientation  (g  <  .05)  except  360  and  440  MHz.  Exposure  of 
rats  to  360  MHz  radiation  in  this  orientation  was  also  different  from  other 
frequency  exposures  (g  <  .05)  except  at  260,  440,  460,  480,  and  500  MHz.  These 
results  indicate  that  for  the  passive  rat,  simple  body  heating  as  measured  by  a 
change  In  colonic  temperature  is' dependent  on  both  frequency  (g  <  .01)  and  body 
orientation  relative  to  the  microwave  fields  within  the  radiation  chamber. 


EXPERIMENT  II:  MATERIALS  AND  METHOD 


Subjects 


Six  male  Long-Evaas  rats  350-380  grams  in  body  weight  were  used  as 
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Figure  2.  Mean  frequency  and  body  orientation  dependent  microwave 
absorption  and  colonic  temperature  rise  in  the  Long  Evans  rat. 
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subjects.  All  animals  were  deprived  of  food  to,  and  maintained  at,  85  percent 
of  free-feeding  body  weights. 


Apparatus 

The  parallel  plate  irradiation  system  and  rat  holder  were  identical  to 
that  of  Experiment  I.  The  rat  holder  was  modified  by  attaching  a  small  Plexi¬ 
glas  response  lever  and  small  1.25  cm  diameter  glass  food  cup.  Since  standard 
metallic  microswitches  cannot  be  used  in  a  microwave  field,  an  alternate 
system  of  monitoring  lever-pressing  behavior  was  devised.  This  was  accomplished 
by  passing  a  concentrated  light  beam  through  47  mil  plastic  Dupont  Crofon  fiber 
optics.  The  fiber  optics  were  arranged  in  such  a  way  that  each  lever  press 
made  by  the  rat  interrupted  the  passage  of  the  light  beam.  A  sensitive  photo¬ 
cell  arrangement  outside  of  the  radiation  chamber  could  then  detect  each  inter¬ 
ruption  of  the  light  beam  and  consequently  each  lever  press.  A  Lehigh  Valley 
Electronics  Corporation  pellet  feeder  was  mounted  outside  the  radiation  cham¬ 
ber  and  delivered  45  mg  Noyes  food  pellets  to  the  glass  cup  of  the  rat  holder 
inside  the  radiation  chamber  via  a  1.27  cm  x  1  m  section  of  Tygon  flexible 
plastic  tubing.  Scheduling  and  recording  of  behavioral  test  sessions  was  ac¬ 
complished  via  a  system  of  relay  circuitry,  digital  printing  counters,  and 
cumulative  response  recorders. 


Procedure 

All  six  rats  were  trained  to  lever  press  for  food  pellets  on  a  random 
interval  thirty-second  schedule  of  reinforcement.  The  clock  time  used  for  the 
scheduling  was  three  seconds,  with  a  probability  of  .10  food-pellet  delivery. 
With  this  arrangement  a  food  pellet  could  be  delivered  from  three  seconds  to 
every  180  seconds,  with  the  average  time  of  each  food-pellet  delivery  being 
every  thirty  seconds.  Each  rat  was  given  daily  thirty-minute  lever-pressing 
sessions  until  a  stable  rate  of  responding  was  achieved.  Lever-pressing 
behavior  was  considered  stable  when  the  number  of  lever  presses  during  the 
first  five  minutes  of  a  session  varied  less  than  10  percent  from  a  similar 
measure  of  the  previous  daily  session.  On  the  average,  approximately  twenty 
daily  lever-pressing  sessions  were  required  to  achieve  stable  response  rates. 

Once  stable  responding  rates  were  achieved,  all  animals  were  repeatedly 
and  randomly  exposed  to  different  conditions  of  microwave  radiation.  For  all 
microwave  exposures,  the  field  density  within  the  chamber  at  the  position  of 
the  rat  was  maintained  at  25  mW/cm^  as  measured  by  the  field  intensity  meter. 

All  rats  were  exposed  randomly  to  360,  480,  and  500  MHz  microwave  radiation, 
with  both  the  length  of  the  rat  body  parallel  to  the  electric  field  (e||  L)  and 
direction  of  wave  propagation  through  the  microwave  chamber  (c||  L) .  In  this 
fashion  each  rat  was  exposed  to  the  six  possible  combinations  of  three  frequen¬ 
cies  and  two  body  orientations  of  microwave  exposure  while  performing  a  stable 
lever-pressing  task.  Radiation  exposures  with  the  rat  body  parallel  to  the 
magnetic  field  were  not  done  due  to  the  difficulty  in  rotating  the  radiation 
chamber.  Retraining  sessions  were  given  each  rat  after  each  microwave  exposure 
to  maintain  stable  rates  of  responding  before  the  next  microwave  exposure  and  to 
check  for  possible  carryover  effects.  Microwave  irradiation  began  at  the  begin¬ 
ning  of  the  sixth  minute  of  each  exposure  session  and  terminated  at  the  end  of 
the  first  one-minute  interval  during  which  responding  fell  below  one-third  of 
the  rat's  normal  response  rate  per  minute  determined  from  the  previous  daily 
session.  If  the  rat  responded  throughout  the  session  at  a  rate  above  one-third 
the  baseline,  irradiation  ceased  at  the  end  of  the  thirty-minute  session, 
resulting  In  a  maximum  exposure  of  25  minutes.  The  primary  measure  of  lever¬ 
pressing  disruption  was  time  to  work  stoppage  defined  as  the  time  from  radiation 
onset  until  the  end  of  the  minute  during  which  responding  fell  below  one-third 
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of  the  baseline  rate.  Each  rat  was  removed  from  tie  radiation  chamber  and 
animal  holder  at  work  stoppage  or  at  the  end  of  the  session. 

Colonic  temperature  measurements  were  made  at  the  beginning  of  the  radia¬ 
tion  sessions  and  again  at  work  stoppage  or  the  end  of  an  exposure  session  to 
determine  the  amount  of  microwave- induced  heating  for  each  of  the  exposure 
conditions.  Water  intake  during  the  24-hour  period  between  each  daily  session 
while  each  animal  was  in  its  home  cage  was  measured  and  recorded.  Relative 
humidity  during  the  radiation  sessions  varied  between  12  percent  and  40  per¬ 
cent.  Ambient  temperature  of  the  radiation  chamber  was  controlled  by  room 
temperature  (70-72°  F). 

Because  of  a  lack  of  closed  circuit  television  equipment,  visual  observa¬ 
tion  of  the  animal's  performance  during  microwave  radiation  was  not  possible. 

EXPERIMENT  II:  RESULTS 

Frequency  of  the  microwave  radiation  and  orientation  in  the  microwave 
field  both  proved  to  be  of  major  importance  in  determining  disruption  of  lever¬ 
pressing  behavior.  Figure  3  shows  the  time  to  work  stoppage  from  radiation  on¬ 
set  for  the  frequency  and  body  orientation  parameters.  For  microwave  exposure 
sessions  with  the  animal's  length  aligned  parallel  to  the  direction  of  wave 
propagation  through  the  chamber  (c||  L) ,  nearly  all  animals  lever  pressed  at 
near  normal  rates  throughout  the  session.  Only  one  animal,  and  this  at  500  MHz 
radiation  in  the  c||  L  orientation,'  reduced  responding  to  the  work-stoppage 
criterion.  With  the  length  of  the  rat  body  aligned  parallel  to  the  electric 
field  (E ||  L),  however,  time  to  work  stoppage  was  much  shorter.  Average  values 
of  time  to  work  stoppage  and  standard  deviations  are  shown  in  Fig.  3.  Animals 
exposed  to  360  and  480  MHz  radiation  in  the  E | [  L  orientation  appeared  nearly 
normal  upon  removal  from  the  microwave  chamber. 

Animals  exposed  to  500  MHz  radiation  in  the  E||  L  orientation  appeared 
flushed  and  remained  in  a  prone  position  with  the  body  extended.  The  condition 
of  the  animal  appeared  similar  to  the  description  of  what  Justesen  and  King  [9] 
termed  "flaccid  paralysis".  Four  of  the  animals  exposed  under  these  conditions 
had  urinated  profusely  in  the  animal  holder.  The  fur  on  the  underside  of  the 
animal's  body  was  wet  in  spite  of  the  slots  made  in  the  holder  to  release  urine 

from  the  holder  and  away  from  the  microwave  fields.  The  rats  were  clearly  heat 

stressed.  Each  animal,  however,  quickly  regained  its  upright  position  within 
minutes  after  removal  from  the  radiation  chamber.  A  randomized  block  factorial 
analysis  of  variance  [8]  was  used  to  evaluate  the  effects  of  the  treatments. 
This  test  indicated  a  significant  difference  between  the  frequency  of  radiation 
during  exposure  (F  =  17.68,  df  =  2/25,  g  <  .001),  the  animal  orientation  in  the 
microwave  field  (F  =  119.24,  df  =  1/25,  g  <  .001),  and  a  frequency  by  orienta¬ 
tion  interaction  (F  =  12.37,  df  =  2/25,  p  <  .001). 

A  Newman  Kuels  test  found  no  difference  in  the  times  to  work  stoppage  for 
the  exposures  with  the  length  of  the  animal  aligned  parallel  to  the  direction 
of  the  microwave  propagation  through  the  chamber  (c||  L)  at  360,  480,  and  500 
MHz.  Each  of  these  exposure  conditions  was,  however,  different  (p  <  .05  or 
greater)  from  the  exposures  with  the  animal  aligned  with  the  length  of  the  body 
parallel  to  the  electric  field  (E ||  L)  at  360,  480,  and  500  MHz.  With  the 
animal  aligned  parallel  to  the  electric  field,  times  to  work  stoppage  at  360 
and  480  MHz  were  not  significantly  different.  Exposure  to  500  MHz  radiation  ir. 
this  orientation  was,  however,  different  (g  <  .01)  from  all  other  exposure 
conditions. 

The  amounts  of  body  heating,  as  measured  by  colonic  temperature  change 
during  irradiation,  paralleled  the  behavioral  effects.  Figure  4  shows  the 
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relative  absorption  of  microwave  energy  as  determined  calorimet r i cal ] y  and 
given  as  mW/gm.  Also  shown  in  Fig.  4  is  a  scale  showing  the  animal  body 
temperature  fisc  (centigrade  per  minute).  For  animals  aligned  parallel  to  the 
wave  propagation  (c||L),  body  heating  was  minimal.  For  animals  aligned  paral¬ 
lel  to  the  electric  field  (H ] |  L) ,  however,  heating  was  significant  and  in¬ 
creased  as  the  frequency  of  radiation  increased.  A  randomized  block  factorial 
analysis  of  variance  [8]  was  used  to  evaluate  the  effect  of  the  treatments  on 
rat  whole  body  heating  as  measured  rectally.  This  test  indicated  a  difference 
between  the  frequency  of  radiation  during  exposure  (F  =  21.70,  df  =  2/25,  p  < 
.01),  the  animal's  orientation  in  the  microwave  field  (F  =  129.81,  df  =  1/25, 
p  <  .01),  and  a  frequency  of  radiation  by  animal  orientation  interaction  (F  = 
8.92,  df  =  2/25,  p  <  .01). 

A  Newman  Fuels  test  indicated  that  rat  body  heating  during  exposure  was 
not  different  at  360,  480,  or  500  MHz  when  the  rat  body  was  parallel  to  direc¬ 
tion  of  wave  propagation  (c||  L)  through  the  radiation  chamber.  With  the  rat 
aligned  parallel  to  the  electric  field  (E j |  L) ,  however,  500  MHz  exposures  were 
different  (p  <  .01)  from  360  and  480  MHz  and  also  from  the  exposures  in  the 
wave  propagation  orientation  at  all  frequencies.  In  the  electric  field  orienta¬ 
tion,  rat  body  heating  at  360  and  480  MHz  were  not  different  from  each  other 
but  were  clearly  different  from  exposures  in  the  wave  of  propagation  orienta¬ 
tion  (g  <  .05). 

The  measures  of  time  to  work  stoppage  and  rat  body  heating  expressed  as 
mW/gm  followed  each  other  very  closely  in  a  reverse  fashion.  As  body  heating 
increased,  time  to  work  stoppage  became  much  shorter.  The.  correlation  coef¬ 
ficient  between  the  time  to  work  stoppage  and  mW/gm  measures  for  all  animal 
radiation  exposures  was  -.79  (p  <  .01). 

Individual  records  of  response  rate  are  shown  in  Figs.  5  and  6.  Respond¬ 
ing  during  radiation  in  the  e||  L  orientation  at  500  MHz  is  shown  in  Fig.  5  with 
comparable  baseline  performance  and  recovery  responding  the  day  following  radia¬ 
tion  exposure  shown.  Responding  for  the  same  animal  during  radiation  in  the 
C||  L  orientation  at  500  MHz  is  shown  in  Fig.  6.  In  both  figures  radiation  on¬ 
set  is  shown  (RFon)  as  well  as  time  to  work  stoppage  and  radiation  termination 
«*off). 

An  examination  of  the  average  number  of  responses  per  minute  during  micro- 
wave  irradiation  revealed  an  interesting  finding.  A  comparison  of  each  animal's 
average  response  rate  during  microwave  irradiation  until  time  to  work  stoppage 
or  the  end  of  the  exposure  session  with  a  similar  measure  from  the  previous 
day's  baseline  performance  showed  no  difference.  Figure  7  shows  this  finding. 
Total  lever-press  rates  for  all  exposure  conditions  were  transformed  into  sup¬ 
pression  of  responding  ratios.  The  suppression  ratios  were  used  as  an  index  of 
the  effect  of  microwave  radiation  exposure  on  responding.  Each  suppression 
ratio  was  indexed  by  the  ratio  of  B/A  +  B  where  A  represents  the  total  number 
of  lever  presses  during  baseline  performance.  The  score  B  represents  the  total 
number  of  lever  presses  during  the  period  in  radiation.  Thus  a  ratio  of  0.50 
would  indicate  no  relative  change  of  lever  presses  during  microwave  irradiation 
compared  to  the  same  measure  from  baseline  responding  with  no  microwave  irradi¬ 
ation  present.  A  ratio  of  0.00  would  indicate  a  complete  cessation  of  lever 
pressing  during  microwave  radiation  exposure.  At  each  of  the  microwave  exposure 
conditions,  responding  remained  fairly  stable  with  an  abrupt  drop  or  cessation 
of  responding  during  the  one-minute  interval  when  the  work  stoppage  criterion 
was  met.  A  randomized  block  factorial  analysis  of  variance  [8]  revealed  that 
suppression  of  responding  up  until  the  work  stoppage  criterion  was  not  different 
for  either  the  frequency  of  microwave  radiation  (g  >  .10)  or  the  orientation  of 
the  animal  in  the  microwave  field  (p  >  .10). 


Similarly  a  measure  of  the  amount  of  water  consumed  by  the  rat  between 


Figure  5.  Cumulative  response  record  of  rat 
MW-20,  500  MHz,  E j |  L,  radiation  exposure. 
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Figure  7.  Mean  suppression  ratios  of  responding  during  microwave  exposure 


daily  lever-pressing  sessions  was  also  not  influenced  by  exposure  to  microwave 
radiation.  The  amount  of  water  consumed  by  each  rat  after  microwave  exposure 
until  the  next  daily  session  was  measured.  This  amount  of  water  intake  (cc) 
was  then  compared  to  a  similar  measure  of  water  intake  before  microwave  expo¬ 
sure.  Figure  8  shows  Lhe  mean  amount  and  standard  deviations  of  water  consumed 
after  microwave  radiation  exposure  as  a  percentage  cf  each  animal's  baseline 
water  intake.  A  randomized  block  factorial  analysis  of  variance  [8]  revealed 
no  significant  difference  in  water  intake  as  a  function  of  frequency  of  micro- 
wave  radiation  (g  >  .10)  and  orientation  of  the  animal  in  the  microwave  field 

(g  >  .10). 

To  determine  whether  the  effects  of  a  particular  radiation  session  had  any 
effect  on  subsequent  responding  in  recovery  test  sessions,  each  animal's  per¬ 
formance  during  the  first  five  minutes  of  the  baseline  session  was  compared  to 
performance  on  the  day  after  or  two  days  after  microwave  radiation  exposure. 
Standard  suppression  ratios  of  B/A  +  B  were  used.  The  animal's  responding  dur¬ 
ing  the  first  five  minutes  of  the  baseline  sessions  is  represented  by  the  score 
A,  whereas  total  responding  during  the  first  five  minutes  of  the  recovery  ses¬ 
sion,  either  one  or  two  days  after  microwave  irradiation,  is  represented  as  the 
score  B.  Thus  a  ratio  of  .50  would  indicate  no  relative  change  in  responding 
during  the  recovery  sessions  compared  to  performance  during  previous  baseline 
measures.  Figure  9  shows  the  result  of  this  comparison,  suggesting  little  car¬ 
ryover  effect  from  radiation  sessions  into  recovery  sessions  one  or  two  days 
later.  A  randomized  block  factorial  analysis  of  variance  [8]  revealed  that 
there  were  in  fact  no  significant  differences  in  the  two  performance  measures 
due  to  either  the  frequency  of  radiation  factor  (g  >  .10)  or  the  orientation  of 
the  animal  in  the  microwave  field  (g  >  .10). 

DISCUSSION 

This  study  was  conducted  to  answer  two  questions.  First,  does  the  fre¬ 
quency  of  microwave  radiation  and  orientation  of  an  animal  in  a  microwave  field 
affect  the  absorption  of  microwave  energy  and  consequent  animal  body  heating? 
Second,  ran  these  factors,  by  producing  differential  heating,  significantly 
alter  the  behavior  of  an  animal  performing  in  a  microwave  field?  The  data 
reported  here  support  an  affirmative  answer  to  both  of  these  questions. 

The  results  of  both  experiments  clearly  indicate  a  strong  effect  of  micro- 
wave  frequency  and  animal  body  orientation  in  determining  microwave  absorption, 
rat  body  Seating,  and  consequent  suppression  of  behavioral  responding.  Animals 
exposed  1>  the  most  absorbant  frequency  and  body  orientation  combination  (50C 
MHz,  E ||  L.  were  clearly  heat  stressed.  Vasodilation  and  assumption  of  a  prone 
body  position  in  these  animals  after  radiation  exposure  were  marked.  Time  to 
work  stoppage  for  animals  exposed  under  these  conditions  was  short.  For  animals 
exposed  to  less  absorbant  frequencies  and  body  orientations,  the  rate  of  body 
heating  was  less  and  consequent  disruption  of  behavior  was  not  evident.  The 
high  correlation  between  rate  of  body  heating  and  time  to  work  stoppage  further 
supports  the  hypothesis  that  differential  microwave  absorption,  body  heating, 
and  behavioral  suppression  can  be  governed  by  the  factors  of  microwave  fre¬ 
quency  and  animal  body  orientation  in  the  microwave  fields. 

The  question  may  be  asked  as  to  how  mild  hyperthermia  may  be  responsible 
for  suppression  of  a  well-trained  behavior  in  the  rat.  It  is  interesting  to 
note  that  the  lever-pressing  animal  maintained  near  normal  average  response 
rates  during  microwave  irradiation.  There  seemed  to  be  an  abrupt  cessation 
of  lever  pressing  during  microwave  exposure,  under  most  absorbant  conditions, 
rather  than  a  gradual  decline  in  behavior  below  the  work-stoppage  criterion. 

This  effect  was  unexpected,  but  in  agreement  with  the  earlier  report  by 
Justescn  and  King  [9]  of  a  similar  finding  in  rats  exposed  to  microwave 


87 


AD-AlOO  029 
UNCLASSIFIED 


UTAH  UNIV  SALT  LAKE  CITY  P/6  w,. 

ELECTROMAGNET I C  ENERGY  ABSORPTION  AND  ITS  DISTRIBUTION  FOR  MAN  .-ETCtUl 
NOV  78  OP  GANDHI,  J  *  D* ANDREA,  M  J  HA6MANN  DAMD17-07-C-Io92 
UTEC-HD-78-109 


o  -j 
o  — 

in  u* 


O  -± 
00  — 


O  _i 

•o  = 
co  «ju 


O  -i 
sO  = 
CO  V 


3  >1  V1NI  831VM 


FREQUENCY  MHz  &  ORIENTATION 


Figure  9.  Carryover  effects  suppression  ratios. 


radiation.  Roberts,  et  at.  ,  [10]  have  shown  that  rats,  heated  by  raising  the 
environmental  temperature,  show  characteristic  changes  in  behavioral  thermo¬ 
regulation.  Rats  heated  in  this  fashion  show  an  initial  increased  level  of 
body  grooming  and  activity  during  the  first  ten  minutes  of  exposure  to  ambient 
temperatures  of  37.7°  C.  Signs  of  body  relaxation  and  extension,  characteris¬ 
tic  of  the  heat-stressed  rat,  occurred  only  after  45  minutes  exposure  to  37.7°C 
ambient  temperature.  The  possibility  then  exists  that  efferent  mechanisms 
responsible  for  grooming  behaviors  as  a  thermoregulatory  response  may  also  be 
elicited  in  the  microwave  heated  rat.  Mild  hyperthermia  induced  by  microwave 
exposure  may  elicit  thermoregulatory  behavior  in  the  rat  such  as  grooming  which 
conflict  with  lever-pressing  behavior.  Thus  time  to  work  stoppage  measured  in 
this  experiment  may  in  fact,  for  some  animals  at  least,  be  measuring  time  to 
elicitation  of  grooming  behavior.  Since  microwave  exposure  may  be  a  more  ef¬ 
ficient  and  rapid  method  of  inducing  mild  hyperthermia,  several  classes  of  the 
classic  thermoregulatory  behaviors  may  be  elicited  simultaneously  or  in  rapid 
succession.  This  latter  possibility  may  be  true  for  the  animals  exposed  under 
the  most  absorbant  radiation  conditions.  Since  direct  observation  of  each 
animal  during  irradiation  could  not  be  done,  further  investigation  of  this 
hypothesis  is  required. 

It  is  possible  that  other  factors  may  be  responsible  for  behavioral  sup¬ 
pression  than  just  mild  hyperthermia  induction  and  elicitation  of  behavioral 
thermoregulation.  One  may  speculate  that  evaporative  water  loss  during  micro- 
wave  exposure  was  enhanced  and  that  consequent  thirst  induction  may  contribute 
to  behavioral  suppression.  A  finding  of  no  change  in  normal  water  Intake  over 
the  24-hour  period  after  microwave  exposure,  however,  seems  contrary  to  this 
hypothesis.  The  lack  of  any  evidence  of  proactive  effects  due  to  microwave 
exposure  suggests  that  the  most  reliable  effect  of  such  exposure  was  body  heat¬ 
ing  sufficient  to  suppress  behavior  but  not  intense  enough  to  produce  actual 
physical  thermogenic  damage  to  the  animal.  This  may  suggest  that  the  primary 
effect  of  moderate  intensity  microwave  exposure  is  to  induce  mild  hyperthermia 
and  consequent  suppression  of  behavior,  both  of  which  are  transient  in  nature. 

It  is  important  to  note  in  Experiment  I  that  even  though  each  animal  was 
allowed  sufficient  time  for  colonic  thermal  equilibration  before  exposure  to 
microwave  radiation,  each  animal  was  alert  during  the  exposure  period.  It  is 
well  known  that  emotional  reactions  produce  psychogenic  effects  in  the  rodent 
and  that  temperature  regulation  can  be  altered  by  such  factors.  It  is  difficult 
to  rule  out  such  reactions  to  the  environment  when  using  colonic  temperature  in 
the  alert  animal  as  the  dependent  variable.  The  use  of  habituation  techniques 
such  as  familiarizing  the  animal  to  the  animal  holder  and  radiation  chamber  may 
help  minimize  such  reactions.  Nevertheless  the  animal  cannot  be  readily 
habituated  to  the  microwave  exposure  without  long-term  adaptation  to  induced 
heating  effects.  Determination  of  absorbed  microwave  power  as  described  in 
this  report  must  be  distinguished  as  a  dynamic  measure,  since  the  rat  is  cap¬ 
able  of  active  thermoregulation  during  irradiation. 
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I.  INTRODUCTION 

In  electromagnetics,  discretization  for  transforming  an  inte¬ 
gral  equation  to  a  matrix  equation  is  commonly  accomplished 
using  pulse  functions  as  a  basis  [  1  ]  — [ 4J .  More  elaborate  bases 
such  as  modal  fields  may  be  justified  in  a  particular  problem 
but  are  avoided  in  the  matrix  formulation  for  the  general 
scattering  body. 

Discretization  with  pulse  functions  requires  that  the  scatter¬ 
ing  body  be  partitioned  into  a  number  N  of  cells,  where  N  is 
large  enough  that  complex  permittivity  and  the  complex  time- 
independent  electric  field  may  be  assumed  constant  within 
each  subvolume.  Procedures  for  improving  convergence  with 
the  general  scattering  body  are  of  interest  since  the  cost  of 
computation  is  a  rapidly  increasing  function  of  N. 

This  communication  describes  two  new  procedures  which 
have  been  found  to  improve  convergence  in  solution  of  the 
two-dimensional  problem  of  transverse  magnetic  (TM)  exci¬ 
tation  of  an  infinite  cylinder  of  arbitrary  cross-section  shapt 
Convergence  is  improved  because  the  calculations  include  part 
of  the  variation  of  the  field  within  each  cell. 


II.  REVIEW  OF  EARLIER  SOLUTIONS  USING 
PULSE  FUNCTIONS 

For  TM  excitation  of  an  infinite  cylinder,  both  incident  and 
scattered  fields  have  E  directed  parallel  to  the  infinite  dimen¬ 
sion  (z  axis).  The  electric  field  is  described  by  a  scalar  integral 
equation: 

y )  =  Ej(x,  y)  -  J  J  [er(x  ,  y  )~  I  ] 

•  Ez(x\  y')Ho(2>(k0p)  dx'  dv' ,  ( 1  I 
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where  (x,  y)  and  (x' ,  y')  are  coordinates  of  the  observation 
point  and  source  point,  respectively ;  p  is  the  distance  between 
the  two  points;  Et‘  and  E,  are  incident  and  total  electric  field 
intensity,  respectively,  with  e‘u>t  time  dependence;  er(x’,  y')  is 
complex  permittivity  at  the  source  point  relative  to  free  space; 
*0  =  wVPo£o  i  anc*  Ho <2)(*oP)  is  the  Hankel  function  of  zero 
order. 

The  discrete  analog  of  (1)  consists  of  the  A^-by-A' system  of 
linear  equations: 

N 

=  m  =  1 .  2,  '**,  A.  (2) 

n  - 1 

Richmond  performed  the  discretization  using  pulse  func¬ 
tions  [  1 J .  Each  square  cell  was  approximated  by  a  circle  of 
radius  a  having  an  equal  area  to  simplify  integration  of  the 
Bessel  functions.  His  expressions  for  the  matrix  elements 
follow: 

(3) 

Amn  =— f-(ern  -  1  Vl(MVA)<!?)(*oPmn).  m  *  n. 

(4) 


=  1  +  (Cn 


1) 


ink  na 

i  + - «t<2>(M) 


001 8-926X/78/0900-0743S00.75  ©  1978  IEEE 


-  92 


744 


IhlK  TH  ANSACTIONS  ON  ANTI  NNAS  AND  HkOt'ACi  ANON .  VOI..  AH  26,  NO  5,  ShHTKMBKK  1978 


Fig.  l.  Notation  used  in  evaluation  of  off-diagonal  mainx  elements 


In  this  .method,  the  electric  field  intensity  is  considered  to 
be  constant  over  the  area  of  a  cell,  an  approximation  that 
results  in  significant  error  if  the  cells  are  not  very  small.  Larger 
cells  can  be  used  if  the  variation  ot  the  electric  field  is  included 
in  the  calculations.  A  method  of  approximating  the  variation 
of  the  electric  field,  called  “plane-wave  correction,”  is  described 
in  the  next  section. 

Ill  PLANE-WAVE  CORRECTION 

Variation  of  ZTz  is  initially  unknown,  but  we  may  make  the 
approximation  of  representing  the  fields  within  each  cell  by  a 

superposition  of  plane  waves.  All  member  plane  waves  are  p  1S  al  the  center  of  the  mth  cen  as  shown  ln  Flg  ]  Let  Pmn 

required  to  have  propagation  vectors  of  magnitude  k  corre-  an(j  /  be  distances  from  the  center  of  the  nth  cell  to  the 

sponding  to  the  frequency  of  the  incident  wave  and  local  observation  point  and  source  point,  respectively.  If  the  angle 

complex  permittivity  of  the  dielectric  body.  If  a,  and  fi,  repre-  between  pmn  and  r  is  0  ,  then  the  distance  between  the  source 

sent  orientation  of  propagation  and  amplitude  of  the  ith  com-  point  anJ  observation  point  is 

ponent,  respectively,  the  superposition  of  place  waves  may  be 
written  as 


Pmn2  -  >'Pmn  COS  O’. 


CM 
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(5) 


If  er  is  constant  within  the  nth  cell  and  the  approximation  of 
circular  cells  is  used  in  (1).  then  the  off-diagonal  matrix 
If  complex  permittivity  is  constant  within  a  cell  and  the  elements  are 
approximation  of  circular  cells  is  us- .'  in  tl ),  then  the  diagonal 
matrix  elements  are 
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Using  (5)  for  /■*( r  .  0  V  neglevUng  terms  higher  than  second 
order  in  kr'  in  the  expansion  of  the  exponential,  and  perform¬ 
ing  the  O'  integrations, 
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Performing  the  integrations  over  the  Hankel  functions  gives, 
for  the  diagonal  elements. 


Using  (5)  for  £■;(>'.  O')  and  neglecting  terms  higher  than  second 
eder  in  kr'  in  the  expansion  of  the  exponential. 
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In  evaluating  the  off-diagonal  matrix  elements,  suppose  that 
the  source  point  /’’  is  in  the  nth  cell  and  the  observation  point 
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Performing  the  d'  integrations  in  (12),  and  neglecting  terms  If  quadratic  and  higher  order  terms  in  16)  are  neglected,  ihen 
with  I  l  |  >  1 , 
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Since  only  the  last  member  in  il'i  leper  ter  ■ 
unknown  coefficients  in  the  plane  wave  expanu-m  *>  ..  . 
delete  the  term  and  use  the  first  three  mrmt-ei-  t  •  ■>.< 
diagonal  matrix  elements. 
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The  summations  in  all  three  integrands  of  (13)  contain  , 

factors  independent  of  the  summation  index  that  may  be 

removed  from  the  summations.  Also  changing  to  a  dummy  The  use  of  (18)  is  justified  only  if  the  Ian  member  of  (  1  7|  is 
variable,  small.  The  fractional  error  in  Amn  due  to  neglecting  the  last 

_ member  of  ( 1 7)  is 
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Since  the  magnitude  of  the  ratio  of  the  Hankel  functions  is 
greatest  for  a  small  argument,  we  set  pmn  =  \A ra  for  adjacent 
cells.  Expanding  the  Bessel  functions  for  small  arguments  and 
using  ( 1 6)  to  estimate  y,  we  obtain 
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where  T  *=  1 .78107. 

In  the  plane  wave  correction  method,  (8)  and  (18)  are  used 
for  the  matrix  elements,  and  the  solution  corresponds  to  the 
(14)  value  of  Ez  at  each  cell  center.  Numerical  considerations 
suggest  that  deletion  of  the  last  member  in  (17)  is  the  principal 
source  of  error.  Equation  (20)  requires  that  if  4\Jrr 
In  (2/(rVffk0a))  exceeds  the  condition  number  of  the  matrix, 
then  the  error  in  Et  calculated  for  the  cell  centers  will  be  less 
than  the  variation  of  Et  within  each  cell  relative  to  values  at 
the  cell  centers.  Hence  the  plane  wave  correction  will  account 
for  much  of  the  variation  of  Et  within  each  cell.  We  must  have 
a  <  0,0876\o  or  the  criterion  will  not  be  satisfied  even  for  a 
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condition  number  of  unity.  Similar  convergence  criteria  are 
not  available  in  the  method  using  pulse  functions  that  is 
described  in  Section  11.  A  second  method  of  approximating 
the  variation  of  the  electric  field,  called  "cylindrical-cell 
correction."  is  described  in  the  next  section. 


IV.  CYLINDRICAL-CELL  CORRECTION 


In  the  methods  described  in  Sections  II  and  III,  the  scatter¬ 
ing  body  is  divided  into  cells  that  are  approximated  by  circular 
cylinders  in  calculation  of  the  matrix  elements.  The  method  of 
cylindrical-cell  correction,  which  is  presented  for  the  first  time 
in  this  communication,  emphasizes  the  properties  of  such  a 
model  of  the  original  body. 

Let  the  cross  section  of  the  scatterer  be  divided  into  square 
cells,  each  having  side  s.  For  the  model,  replace  each  cell  with 
a  circular  cylinder  having  radius  a  and  assume  properties  of 
free  space  between  the  cylinders.  For  the  chosen  geometry  [61 , 
the  effective  relative  complex  permittivity  creff  of  the  model 
and  the  complex  permittivity  er’  of  each  cylinder  are  related  by 


crell  -  I  +  —  Ifc'r 
S i 


I  ). 


(21  1 


It  we  require  that  the  effective  relative  complex  permittivity 
of  the  model  equal  that  of  the  original  scatterer,  then 


rra 


used  to  calculate  the  diagonal  matrix  elements: 

in  , 

4  m  m  ~  ~~  [  k  uJ  i  tk  a  *  2  T  ) 

~  koaJ0(k'a)H .  <  ;<, i 

From  (10)  and  (24).  the  off-diagonal  matrix  elements  are  given 
by 


Mo2(Cr 


1) 


A’ 

J  fl 


dr 


J0(k  r  )//(/  2  *(  k0sm  „ 


I  JO' 


Using  Graf’s  addition  theorem  [5  ]  and  the  notation  of  Fig.  1. 
jk  2  * 

-4  mn  =  (£>n  ~  O  ///*  2  ’( k  oPm  n  I 

4 

•  f  Jt{k0r')JoV<'r')r'  Jr' 

■>0 

Performing  the  9'  integration,  which  collapses  the  summation, 
and  using  a  relationship  from  [7]  for  the  >■'  integration. 
obtain  the  following  expression  used  to  calculate  the  oil- 
diagonal  matrix  elements: 


/ 


cos  10  JO 


i  2S  i 


If  the  electrical  size  of  the  cylinders  is  sufficiently  small, 
the  inhomogeneities  of  the  model  will  have  no  significant 
effect,  so  that  scattering  from  the  model  satisfying  (22)  will 
duplicate  that  of  the  physical  scatterer.  The  model  is  useful 
since  we  may  approximate  the  variation  of  the  electric  field 
within  each  cylinder.  A  circular  cylinder  with  TM  excitation 
will  have  internal  fields  given  by 

~  2  k’nJnikr)  cos  (nd  +  cn),  (23) 

n  =  0 


-4  m  n  =  ~  th(2)(koPmn)[k'aJi(k'j)J0(k0a) 

-  k0uJ0(k'a)Jl(k0a)\ .  m  =£  n 

In  the  cylindrical-cell  correction  method,  (26)  and  (24)  are 
used  for  the  matrix  elements,  and  the  solution  corresponds  to 
the  value  of  at  each  cell  center.  A  solution  is  defined  for 
any  ratio  of  a/s,  but  only  two  cases  have  been  studied:  1)  tan¬ 
gent  cylindrical  cells 


where  the  bn  and  c„  are  determined  by  the  excitation.  If  the 
cylinder  is  sufficiently  small,  the  zero-order  term  will  domi¬ 
nate,  so  we  may  use  the  approximation 


=  b0J0(kr).  (24) 

From  (6)  and  (24),  the  diagonal  matrix  elements  are  given 


=  1  + 


/*  o2(er 


1) 


F 

Jo 


J0(k‘r')H0<2>k0r'r‘  dr' 


(25) 


j  ,  4 

a  ~  ~  ,  er  =14 —  (cr  —  I ) 
2  7T 


and  2)  overlapping  cylindrical  cells 


r 


For  very  large  or  very  small  values  of  a/s,  the  structure  of  the 
model  differs  significantly  from  that  of  the  scatterer,  so  that 
the  solution,  which  actually  corresponds  to  scattering  from  the 
composite  of  circular  cylinders,  will  differ  from  that  for  the 
desired  scatterer.  Advantages  of  the  method  will  be  seen  in  the 
example  in  Section  VI. 

V.  CALCULATION  OF  POWER  ABSORPTION 


where  k'  =y/erm  k0. 

Performing  the  O'  integration  and  using  a  relationship  from 
[7]  for  the  r  integration,  we  obtain  the  following  expression 


With  both  plane-wave  correction  and  cylindrical-cell 
correction,  the  solution  corresponds  to  the  value  of  t/{  at  each 
cell  center.  The  plane-wave  correction  leaves  the  variation  of 
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/_.  within  each  cell  undefined,  hut  with  cylindrical-cell 
correction,  (24)  describes  the  approximate  variation  of  /•_, 
within  each  circular  cylinder  of  the  model. 

If  the  value  of  power  absorption  is  desired,  then  1/2  o * 
may  he  computed  within  each  cell  using  the  values  of  A*  found 
using  pulse  functions  or  plane-wave  corrections.  In  this  section 
the  approximate  variation  of  A.'*  within  each  cell  is  used  for  an 
improved  expression  for  power  absorption  if  cylindrical-cell 
correction  is  used. 

Total  power  absorption  in  one  circular  cylinder  is  given  by 


/T'' 


l:,h2*r  JrJO. 


(30) 


From  t  221  the  conductivity  of  the  cylinder  is  related  to  that  of 
the  scatterer  by 


(31) 


Substituting  ( 24)  and  ( 3  1 )  in  1 30)  gives 


s'-n  fa 

I  ”  hohO*  I  Joik  rVo*(k  r)r  dr. 

u-  J0 


(32) 


Using  a  relationship  from  [  7 1  for  the  r  integration  and  simpli¬ 
fying.  the  average  absorbed  power  density  within  the  cell  is 
given  by 


Im  [/.  ’uy(,*(A  u)i/j(A'a)| 

Re  (k  a)  Im  (ka) 


(33) 


where  /  .  is  the  calculated  value  of  / 
using  cylindrical-cell  correction. 


found  for  the  cell  center 


Fig.  2.  Variation  of  average  specific  absorption  rate  with  frequency 
for  infinite  cylinder  model  of  man.  Incident  power  density  is 
I  mW/cm2. 


VI.  NUMERICAL.  F XAMPLF:  TIIF.  INFINITE  CYLINDER 
MODEL  OF  MAN 

The  infinite  cylinder  has  been  suggested  for  use  as  a  model 
tor  body  extremities  or  the  chest  cavity  in  the  evaluation  of 
biological  ha/.ards  from  electromagnetic  radiation  [81.  Approx¬ 
imation  of  an  extremity  or  the  whole  body  by  an  infinite 
cylinder  may  be  made  using  the  ratio  of  volume  to  length  for 
the  cross-sectional  area. 

Ihe  model  considered  in  this  example  is  an  infinitely  long 
cylinder  with  square  cross  section  of  20  by  20  cm,  having  TM 
excitation  by  a  plane  wave  incident  perpendicular  to  one  of 
the  lour  congruent  faces.  The  body  is  homogeneous  with 
values  of  permittivity  and  conductivity  that  are  typical  for 
muscle,  skin,  and  other  tissues  having  high  water  content  [9] . 

The  methods  described  in  this  communication  have  been 
used  to  calculate  average  specific  absorbed  power  for  the 
model  Suitability  of  the  model  can  only  be  justified  at  high 
trequencies  where  end  effects  may  be  neglected,  but  calcu¬ 
lations  have  been  made  over  a  wide  range  of  frequencies  to 
allow  comparison  of  convergence  with  the  different  methods. 

fig  2  illustrates  the  frequency  dependence  of  average 
specific  absorption  rate  (SAR)  for  the  model.  The  values 
obtained  using  the  accepted  method  of  pulse  functions  with 
100  cells  are  the  standard  for  comparison.  Values  obtained 
using  pulse  functions  with  81  cells  were  found  to  differ  from 


those  with  100  cells  by  less  than  1  percent  for  frequencies  up 
to  200  MHz.  Values  found  using  9  cells  with  pulse  functions, 
plane  wave  correction,  and  cylindrical-cell  correction  with 
tangent  and  overlapping  cells  are  also  given  in  the  figure.  For  9 
cells,  all  four  calculations  are  in  good  agreement  with  the 
standard  at  low  frequencies,  where  the  cells  have  small  elec¬ 
trical  size  Calculations  with  the  new  methods  appear  tc>  have 
comparable  accuracy  and  appear  to  offer  a  significant  improve¬ 
ment  in  convergence  when  compared  to  pulse  functions. 

The  criterion  a  <  0.0876X0  suggests  that  the  plane  wave 
correction  with  9  cells  will  cease  to  be  useful  for  frequencies 
exceeding  700  MHz,  but  the  exact  limit  depends  on  the  con¬ 
dition  number  as  mentioned  in  Section  111.  Increasing  error  is 
observed  for  frequencies  exceeding  200  MHz  which  is  on- 
sistent  with  the  convergence  criteria. 

VII.  CONCLUSIONS 

Two  methods  have  been  found  which  improve  convergence 
m  solution  of  the  two-dimensional  problem  of  TM  excitation 
of  an  infinite  cylinder  of  arbitrary  cross-section  shape.  The 
extension  to  the  three-dimensional  problem  has  not  yet  been 
completed. 

Convergence  is  improved  because  the  calculations  include 
approximation  for  the  variation  of  the  field  within  each  cell. 
The  new  methods  give  expressions  for  the  matrix  elements 
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that  are  different  from  those  found  using  pulse  functions,  but 
values  of  the  matrix  elements  are  in  good  agreement  for  elec¬ 
trically  small  cells. 

The  matrix  elements  found  with  the  new  methods  do  not 
require  a  substantial  increase  in  computational  effort.  For  a 
homogeneous  scatterer.  if  equal  size  cells  are  used,  the  new 
methods  require  one-time  computation  of  several  additional 
Bessel  functions,  but  time  for  calculation  of  the  matrix 
elements  is  dominated  by  recalculation  of  the  zero-order 
Hankel  function  which  is  needed  with  or  without  the 
corrections  made  in  the  new  methods.  If  the  scatterer  is  not 
homogeneous  or  if  different  cell  sizes  are  used,  then  time  for 
calculation  of  each  matrix  element  is  approximately  doubled 
in  the  new  methods.  Time  spent  in  calculation  of  the  matrix 
elements  is  proportional  to  A'2,  whereas  time  spent  in  solving 
the  matrix  equation  is  proportional  to  A'3 .  so  for  large  numbers 
of  cells  there  is  still  no  significant  increase  in  computation 
time  with  the  new  methods. 
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Abstract — Two  interpolants  are  described  that  may  be  used  to  correct 
the  results  of  a  moment-method  solution  using  pulse  functions  as  a  basis 
and  delta  functions  for  testing.  The  interpolants  allow  for  some  of  the 
variation  of  the  fields  within  each  cell  and  thereby  increase  accuracy  and 
improve  convergence.  The  interpolants  are  usable  for  a  general  scatterer 
and  typically  require  about  1  percent  of  the  cost  of  the  initial  numerical 
solution. 

I.  Introduction 

OMENT-METHOD  solutions  for  electromagnetic 
scattering  with  a  general  dielectric  object  often  use 
pulse  functions  as  a  basis  and  delta  functions  for  testing 
[ !]-(-]-  The  object  of  this  paper  is  to  show  that  an  inter- 
polant  may  be  used  to  allow  for  some  of  the  variation  of 
the  fields  within  each  cell  and  thereby  increase  accuracy 
and  improve  convergence.  The  computational  expense  of 
using  the  interpolant  is  typically  about  1  percent  of  the 
cost  of  the  initial  numerical  solution. 

A  numerical  solution  using  a  pulse  function  basis  re¬ 
sults  in  a  single  value  representing  £  within  each  cell.  The 
delta  functions  used  for  testing  enforce  the  integral  equa¬ 
tion  at  the  center  of  each  cell  so  that  the  calculated  E 
values  are  most  representative  of  the  cell  centers.  Experi¬ 
mental  tests  have  shown  that  the  error  in  E  calculated  for 
the  cell  centers  is  relatively  small  even  when  adjacent  cells 
have  values  which  may  differ  by  an  order  of  magnitude 
[t>|.  Since  the  values  of  E  may  be  assigned  to  points  in 
space,  interpolation  is  possible. 

The  present  study  has  been  restricted  to  piece-wise 
interpolation  in  which  the  value  of  the  corrected  specific 
absorption  rate  (SAR)  is  calculated  for  one  cell  at  a  time. 
In  this  paper,  two  interpolants  that  have  been  developed 
by  the  authors  are  described.  The  triquadratic  interpolant 
is  usable  when  cell  centroids  are  on  a  Cartesian  product 
mesh.  The  NEWSUD  interpolant  is  useful  for  problems  in 
which  the  cubical  cells  may  have  different  sizes  and 
general  locations. 
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II.  Test  oe  a  One-Dimensional  Interpolant 

Moment-method  solutions  with  any  subsectional  basis, 
such  as  pulse  functions,  require  that  the  scatterer  be 
approximated  by  a  composite  of  cells  [l]-(5].  Since  ana¬ 
lytical  solutions  for  such  composites  are  unknown  in 
three-dimensional  problems,  a  test  of  convergence  has 
been  made  for  a  one-dimensional  problem. 

Consider  the  one-dimensional  problem  of  a  plane  wave 
polarized  in  the  x  direction  incident  upon  a  dielectric  slab 
which  extends  from  ;  =  0  to  z  =  a.  The  electric  field  is 
described  by  a  scalar  integial  equation: 

EA:)-E[(z)-J-^  f\<r(z-)~i)Ex(z)e  ,k'*'dz  (1) 

where  z  and  z  are  coordinates  of  the  observation  point 
and  source  point,  respectively,  p  is  the  distance  between 
the  two  points,  E[  and  Ex  are  the  incident  and  total 
electric-field  intensity,  respectively,  with  e'“'  time  depen¬ 
dence,  er(z’)  is  the  complex  permittivity  at  the  source  point 
relative  to  free  space,  and  k0  =  wV  p,,t0. 

The  discrete  analogue  of  (1)  consists  of  the  N-by-N 
system  of  linear  equations: 

A 

2  £«.*£«»,.  rn=1.2 . N  (2) 

n  «  1 

where  the  dielectric  is  partitioned  into  ,V  cells  which  are 
thin  slabs,  and  n  and  m  are  indices  for  the  cells. 

The  matrix  elements  are  readily  evaluated  using  ( I )  with 
the  procedure  in  [2]: 

Amn=  1  +U„„- |)(1  -e  :)  (3) 

/‘"',-sin  j.  m*n  (4) 

where  A  is  the  thickness  of  each  cell  and  pm„  is  the 

distance  between  the  centers  of  the  ruth  and  nth  cells. 

The  analytical  solution  for  a  homogeneous  dielectric 
slab  is  given  by 

£,U) 

£'(())[  2  \  t,  cos  \  tr  a  -  z)  +  2j  sin  \  t,  A0(  a  -  z )  ] 

2  \  t,  cos  \  t,  k„a  + j(  I  +  «,)  sin  \  t,  k()a 

(5) 
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Fig.  I.  Calculated  average  specific  absorption  rate  for  a  2-cm  thick 
layer  of  muscle  at  2450  MHz.  Incident  power  density  is  1  MW/cm2, 
«r,  “50,  <r  =  2.2  mho/m. 


Fig.  I  gives  the  average  SAR  for  a  2-cm  thick  layer  of 
muscle  at  2450  MHz  for  three  different  computational 
procedures.  The  analytical  value  was  found  by  integration 
of  1/2  oEE *  using  (5).  Numerical  solutions  were  ob¬ 
tained  using  (2)-(4).  Numerical  values  without  the  inter- 
polant  were  found  by  averaging  1  /2  a  E  E*  for  each  cell. 
Numerical  values  were  also  obtained  by  using  a  piece-wise 
quadratic  interpolant  (7)  to  obtain  an  expression  for  £ 
within  each  cell  and  then  integrating  the  expressions  for  a 
volume  average  of  1/2  o  E  E*. 

It  is  readily  seen  in  Fig.  I  that  the  interpolant  causes  a 
significant  improvement  in  convergence  to  the  known 
analytical  solution.  We  may  relate  the  improvement  in 
convergence  to  the  well-known  preference  of  Simpson’s 
rule  to  the  trapezoidal  rule  for  numerical  quadrature. 
Simpson's  rule  requires  fitting  a  piecewise  quadratic 
through  the  data  points  and  generally  gives  a  better  ap¬ 
proximation  of  the  integrand  function  for  greater  ac¬ 
curacy  than  the  trapezoidal  rule. 

Ill  Triqi  adratic  Interpolant 

We  have  developed  a  triquadratic  interpolant  which 
may  be  used  when  the  cell  centroids  are  on  a  Cartesian 
product  mesh.  As  the  name  suggests,  quadratic  interpola¬ 
tion  is  used  parallel  to  each  of  the  three  Cartesian  axes. 

A  ,otal  of  27  cells  is  used  in  the  stencil  for  each 
calculation.  Fig.  2  shows  the  location  of  the  26  cells 
surrounding  the  cell  in  which  the  correction  is  made.  The 
Cartesian  product  is  shown  for  X  =  -  1,  0.  I.  Y  —  -1,0,  1, 
Z=  -  I.  0,  I,  but  scaling  is  readily  used  to  adjust  to  a 
specified  cel!  size. 

Quadratic  interpolation  may  be  used  with  univariate 
data  for  the  function  H X  )  at  the  three  points  .V  =  1 , 0,  i 

by  the  rule 

m>=  H  -  l)|(.V  I )  +  £'(0)(  I  -  ,V  ) 

(l  +  T)+£(l)y(A  +  I).  (6) 
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Fig.  2.  Location  of  cells  in  stencil  for  the  tnquadranc  inierpolanl 

Triquadratic  interpolation  may  be  performed  on  the  sten¬ 
cil  shown  in  Fig.  2  by  using  products  formed  from  three 
univariate  interpolants.  For  example,  the  contribution  of 
the  function  value  at  .V  =  -  I.  >'  =  -  1.  7.  =0  (cell  number 
2)  is 

F{  -  1,  -  1.0)  y  (.V  -  I )  y  (  >  I )(  I  /)(!+/). 

Triquadratic  interpolation  as  described  in  the  last  para¬ 
graph  is  only  usable  for  corrections  in  a  cell  surrounded 
by  26  other  cells,  as  shown  in  Fig.  2  It  is  possible  to  use  a 
single  interpolation  rule  based  on  the  full  27-cell  stencil  if 
E  values  are  estimated  for  unoccupied  positions  by  means 
of  a  series  of  fill-in  rules.  We  have  used  the  following 
series  of  rules  in  which  the  E  values  are  filled  in  the  order 
of  increasing  distance  from  the  central  cell  by  averaging  E 
values  known  in  adjacent  cells. 

1)  Fill  in  the  six  closest  cells  (numbers  5.  II.  13.  15.  17. 
23).  If  any  are  unoccupied,  use  the  value  of  /.’  in  the 
central  cell. 

2)  Fill  in  the  twelve  next  closest  cells  (numbers  2.  4.  6.  8. 
10,  12.  16.  18.  20.  22.  24,  26).  If  anv  are  unoccupied,  use 
one-half  the  sum  of  E  for  the  pair  from  the  six  closest  cells 
of  part  1  which  share  a  face  with  the  unoccupied  cell. 

3)  Fill  in  the  eight  corner  cells  (numbers  I.  3.  7.  9,  |9. 
21.  25,  27).  If  any  are  unoccupied,  use  one  third  the  sum 
of  E  for  the  three  cells  sharing  a  face  with  the  unoccupied 
cell 

Application  of  triquadratic  interpolation  to  the  general 
scatterer  requires  the  following  calculations  be  made  for 
each  cell.  First,  a  sieve  is  used  to  find  which  of  the  26 
surrounding  cells  are  occupies.  Next,  the  fill-in  rules  are 
used  as  required  to  complete  the  stencil.  Finally,  the 
triquadratic  interpolant  is  used  to  evaluate  the  integral  of 
1/2  o  E  E*  over  the  central  cell.  If  the  cells  have  different 
complex  permittivities,  both  the  fill-in  and  interpolation 
rules  must  be  modified  by  suitably  multiplying  b\  the 
ratio  of  permittivities  so  that  the  interpolant  has  both  D 
normal  and  E  tangential  to  each  boundary  continuous. 
We  have  made  such  modifications  and  have  used  inter¬ 
polation  in  solutions  for  inhomogeneous  dielectric 
scatterers. 

Since  analytical  solutions  are  not  available  for  arrays  of 
cubes,  a  standard  for  evaluating  results  obtained  by 
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1(8  3.  C  alculated  average  specific  absorption  rate  fo:  a  12-cm  cube  of 
muscle  at  I  MHz  without  interpolation.  Incident  power  densits  is  1 
MW  cm\  er,  =  200C.  n  =  0  336  mho  m. 


TABU  I 

Error  is  Cam  i  t  Atm  At iRvur  Sint  it  k  Absorimios  Raii 
kir  A  12-cm  CT  bf  of  Ml  st  i  i  ai  I  MHz.  t,,  =  2000,  n  =  0.556 
mho/m 
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numerical  methods  is  not  easily  obtained,  but  one  conver¬ 
gence  test  has  been  made  for  the  triquadratic  interpolant. 
hg.  3  gives  the  average  specific  absorption  rate  calculated 
without  interpolation  for  a  12-cm  cube  of  muscle  with 
plane  wave  at  I  MHz.  normally  incident  on  one  face 
I, inear  convergence  is  demonstrated  in  the  figure.  Ex¬ 
trapolation  using  the  values  for  216  and  512  cells  gives  an 
SAR  of  0.3275  x  10  6  W /kg  for  an  incident  power  density 
of  I  MW/cm2.  Table  I  gives  the  estimated  errors  in  SAR 
without  interpolation  and  with  the  triquadratic  interpolant 
found  by  comparison  with  the  extrapolated  value  of  SAR. 
Results  in  Table  I  suggest  that  the  triquadratic  interpolant 
causes  a  significant  improvement  in  convergence. 

IV.  NEWSUD  Interpoi  am 

We  have  developed  the  NEWSUD  interpolant  which 
may  be  used  when  the  cells  have  different  sizes  and/or 
arbitrary  placement.  The  NEWSUD  interpolant  is  more 
general  than  the  triquadratic  interpolant  and  does  not 
require  fill-in  rules  but  is  limited  to  linear  rather  than 
quadratic  correction.  The  name  NEWSUD  is  an  acronym 
of  the  words  “north."  “east,"  “west."  “south,"  “up,"  and 


Tig  4  Configuration  for  linear  interpolation  m  the  NhWSl'D 
method 


“down"  since  the  stencil  contains  only  one  cell  near  each 
of  the  six  faces  of  a  central  cube  for  calculations  within 
that  cube. 

When  the  NEWSUD  interpolant  is  used  for  calcula¬ 
tions  within  a  cell,  first,  one  imaginary  plane  is  passed 
through  each  edge  of  the  cube  and  the  veil  centroid.  The 
planes  d'vide  the  cube  into  six  congruem  pyramids  as  well 
as  dividing  the  exterioi  into  six  corresponding  solid  an¬ 
gles.  Next,  each  solid  angle  is  searched  to  find  a  cell 
centroid  as  close  as  possible  to  the  interpolated  cube. 
Interpolation  is  performed  separately  in  each  of  the  six 
pyramids  using  only  the  centroids  of  the  interpolated  cube 
and  corresponding  external  cell.  If  a  solid  angle  does  not 
contain  an  external  cell  centroid  within  a  reasonable  dis¬ 
tance.  say  one  side  of  the  central  cell  or  less,  the  uninter¬ 
polated  value  of  E  is  used  within  the  corresponding 
pyramid. 

I  he  configuration  for  linear  interpolation  within  a 
py  ramid  in  the  NEW  St  D  method  is  shown  in  lag.  4.  The 
point  at  which  interpolation  is  desired  (  V.  )  .  / 1.  the 
centroid  of  the  external  cell  (  V,.  /,).  and  the  centroid 

of  the  interpolated  cell  which  is  used  as  the  origin  de¬ 
termine  a  plane  which  is  used  for  the  figure. 

The  equation  of  the  plane  perpendicular  to  the  line 
connecting  (.V,.  Yt.  7,)  to  the  origin  and  passing  through 
the  point  (.V.  Y.  7.)  is 

V.V,  +  If,  +  //,  =  try/Vf-E  T,;+Zf  (7) 

where  a  is  the  distance  from  the  origin  to  the  plane.  The 
linear  interpolation  used  in  the  NEWSUD  method  when 
an  external  cell  is  found  corresponding  to  a  pyramid  of 
the  interpolated  cell  is  accomplished  by  the  approxima¬ 
tion  of  identifying  E( A.  Y ,  Z)  with  E(.\".  V",  Z  )  so  that 


FA  A  .  I'./)  - 

V'  i;+  *  i  +  s.  i 


Vf  +  >  f  +  7. 1 


/;  (0. 0.0).  (8) 


Linear  interpolation  defined  by  (7)  and  (8)  is  used  to 
evaluate  the  integral  of  1/2  a  E  E*  in  a  pyramid  for 
which  a  corresponding  external  cell  is  found.  If  the  cells 
have  different  complex  permittivities,  (8)  must  be  mod¬ 
ified  bv  suitably  multiplying  by  the  ratio  of  permittivities 
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Fig  5.  An  impro1  ed  model  of  man  for  numerical  calculations. 


Fig.  6.  Whole-body  SAR  for  a  homogeneous  model  of  man.  F|U—  k 
front-to-back.  incident  intensity  =  I  MW/cm*’. 


so  that  both  D  normal  and  E  tangential  to  a  boundary  are 
continuous. 

Fig.  5  illustrates  a  model  of  man  which  has  been  used 
for  numerical  calculations  [8|.  Moment-method  solutions 
have  been  made  for  the  model  using  pulse  functions  as  a 
basis  and  delta  functions  for  testing.  The  NEWSUD  inter- 
polant  has  been  used  in  calculations  since  the  model  has 
cells  which  are  cubes  of  different  sizes.  Fig.  6  shows  the 
frequency  dependence  of  the  average  SAR  calculated  for 
the  model  of  man  with  a'  vertically  polarized  plane  wave 
having  frontal  incidence.  Values  found  both  with  and 
without  interpolation  are  given. 

If  the  local  E  values  are  exact  samples  and  if  there  is 
appreciable  variation  between  adjacent  cells,  then  statisti¬ 
cal  procedures  may  be  used  to  calculate  confidence  limits 
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Fig.  7.  Correction  to  whole-body  SAR  for  homogeneous  model  of 

man 


for  values  of  the  average  SAR  found  without  aid  of 
interpolants  [9).  The  80-percent  confidence  half  widths  for 
average  SAR  values  have  been  calculated  and  plotted  in 
Fig.  7.  The  percent  correction  which  the  NEWSUD  inter- 
polant  makes  is  also  shown  in  Fig.  7. 

The  frequency-dependent  variation  of  the  true  E  must 
cause  increasing  error  in  calculations  of  local  E  at 
frequencies  greater  than  200  MHz  for  the  cell  sizes  used 
with  the  model  of  man  [10].  The  80-percent  confidence 
half  widths  shown  in  Fig.  7  are  not  valid  at  frequencies 
above  200  MHz  since  their  calculation  requires  the 
assumption  that  the  local  E  values  are  exact  samples.  Fig. 
7  shows  that  the  NEWSUD  interpolant  makes  relatively 
little  correction  to  the  SAR  at  low  frequencies  where 
relatively  little  error  is  expected,  but  the  correction  in¬ 
creases  sharply  at  frequencies  for  which  the  error  is  ex¬ 
pected  to  increase. 


V.  Contusions 

The  interpolative  correction  to  calculated  values  of  the 
average  SAR  has  been  tested  in  a  one-dimensional  prob¬ 
lem  in  which  the  analytical  solution  is  available  for  com¬ 
parison  with  numerical  solutions.  A  test  was  also  made  in 
a  three-dimensional  problem  in  which  the  true  solution 
was  estimated  by  extrapolation.  In  both  tests  the  inter¬ 
polant  makes  small  corrections  when  few  cells  are  used 
since  sampling  is  insufficient  to  allow  proper  estimation  of 
the  variation  of  the  fields.  The  fractional  correction  is  also 
observed  to  decrease  as  many  cells  are  used  since  the 
interpolated  and  uninterpolated  solutions  must  converge 
to  the  same  answer.  In  both  tests  the  interpolant  was 
found  to  provide  a  significant  reduction  in  error  for  a 
negligible  increase  in  computational  expense. 

It  is  hoped  that  the  two  interpolants  presented  in  this 
paper  arc  only  the  first  steps  in  the  development  of  a 
procedure  which  will  find  general  usage  as  a  follow-up  to 
moment-method  solutions.  Use  of  the  interpolant  for 
calculation  of  the  external  scattered  fields  and  develop¬ 
ment  of  a  general  interpolant  consistent  with  the  wave 
equation  appear  to  be  the  next  necessary  steps. 
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